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Idealplastische isotrope und orthotrope Platten bei Vollausnutzung 
aller Elemente. 


Von H. Craemer. 


1. Problemstellung. In der Plastostatik der Stabwerke, insbesondere der durchlaufenden 
Balken ist es sehr einfach, die Bemessung so vorzunehmen, da fiir einen bestimmten Be- 
lastungsfall alle Querschnitte auf Biegung voll ausgenutzt sind. Man kann z. B. von einem 
Balken mit tiberall gleichem plastischem Widerstandsmoment ausgehen und fiir diesen die 
FlieBgelenke bestimmen; damit ist dann der Momentenverlauf im ganzen Tragwerk aus 
Gleichgewichtsbetrachtungen allein bestimmbar, d. h. statisch bestimmt. Zwischen den 
FlieBgelenken sind dabei die Querschnitte nicht plastiziert, und es steht nichts im Wege, die 
Widerstandsmomente so lange zu verringern, bis alle Schnitte voll ausgenutzt sind, also der 
kontinuierliche Spannungsausgleich? vorliegt. DaB diese Vollausnutzung unter der Hichst- 
last tatsachlich eintritt, folgt aus den Extremalsatzen von Prager und Greenberg ?; dies gilt 
auch fir die in den spateren Abschnitten gebrachten Lésungen fiir Platten. 

Bei einer Platte iiberall gleicher Starke treten an die Stelle der FlieBgelenke die Bruchlinien?; 
beim gleichmafig belasteten Quadrat z. B. sind dies im allgemeinen die Diagonalen. Der 
Spannungszustand in den dazwischen liegenden Bruchstiicken, die nicht plastiziert sind, ist 
aber, im Gegensatz zum Stabwerk, nicht statisch bestimmt. Um den Momentenverlauf in ihnen 
zu bestimmen, mite also die Elastostatik herangezogen werden, was wegen ihrer unregel- 
maBigen Gestalt meist auf praktisch unlésbare Probleme fihrt. 

Die Aufgabe liegt also insoweit anders als bei Stabwerken. Sie ist besonders fiir den Stahl- 
beton von Wichtigkeit und besteht darin, die Platte bei gegebener Belastung und Stiitzung so 
zu bemessen, daf sie in ihren samtlichen Elementen voll ausgenutzt ist. Ihre Auflager setzen 
wir dabei so bemessen voraus, da8 sie die an sie abgegebenen Belastungen und Momonte auf- 
nehmen kénnen, ohne ihrerseits plastische Formanderungen zu erleiden. 


2. Gleichgewichtsbedingungen. Die Gleichgewichtshedingungen eines beliebig gerichteten 
Elements dx, dy in kartesischen Koordinaten (x, y), (Abb. la) sind bekanntlich 


—s 0m, Omxy aa omy Omxy 
ear a dy ’ Ere Fagen (2.1) 

und. 

Oqx ogy _ = 

Tan ae [rare P> (2.2) 
also . 

Omx O?maxy 0? my 

Se ee ve 
eae aay: P (23) 


Hierbei sind gq und m die Querkrafte und Momente je Langeneinheit und p die Flachen- 
belastung. 
Bei Polarkoordinaten und polarer Symmetrie (Abb. 1b) gilt statt (2.3) 


rqg=—m,+—(rm). (2.4) 


3. FlieBzustand isotroper Platten. Fiir jedes Plattenelement gibt es bekanntlich* zwei auf- 
einander senkrechte Richtungen 1 und 2, die Hauptrichtungen, in denen die Biegemomente m, 


1 H. Craemer, Ing.-Arch. 21 (1953), 5S. 187. 

2H. J. Greenberg u. W. Prager, Proc. Am. Soc. Civ. Egrs., Sept. 59, 1952. 
3 K, W. Johansen, Brudlinieteorier, Kopenhagen 1943. 

4 4, Nadai, Elastische Platten, Berlin 1925. 
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und m,, die Hauptmomente, algebraische GréBt- und Kleinstwerte sind und die Schub- | 


momente verschwinden; m sei der absolut gréBte der beiden Werte. 


Als FlieSbedingung itbernehmen wir die bei Stabwerken gebrauchliche; bei voller Plasti- | 


zierung wird also 
h2 
Bee gn oo ayo 
m= Os A 9 ( ) 


wobei h die Plattendicke und o, die Flie8spannung des Baustoffs ist. Die Dehnungen und | 


infolgedessen die Krimmung in Richtung dieses Hauptmoments werden hierbei theoretisch | 


unendlich grof. 
Fiir gerissenen Stahlbeton gilt statt dessen angenahert 


mm afi(n— 28), 3.2 


20> 


wenn f, die Bewehrung je Langeneinheit, o, ihre FlieBspannung und o, der die FlieBspannung | 
des Betons vertretende Wert ist. Hierbei wurde angenommen, dafi die Bewehrung den Schnitt | 
senkrecht kreuzt, doch ist eine Erweiterung fiir beliebige Richtungen méglich +. Entsprechend | 


unserer Aufgabe sind h, f, und mim allgemeinen ortsveranderlich. 


a) Gleichsinnige Vollausnutzung. Als solche bezeichnen wir nun einen Zu- | 


stand, bei dem in allen Elementen in beiden Richtungen die Hauptmomente gleich den | 


FleBmomenten werden: 


ti, ==, == mm. (3.1.1) 


Der Momentenkreis wird dann zum Punkt, sodaB keine ausgezeichneten Hauptrichtungen | 


bestehen und fir ein belie big gerichtetes kartesisches System 


mM, = My =M, May = 0 (3.1.2) 


wird. Damit gehen (2.1) und (2.3) tiber in 


am om ; 
= ee Vy = oy ry (3.1.3) 
0? Og 
(sz ! aa)” Bere (3.1.4) 


Man sieht, dafs diese Gleichung formal mit derjenigen fiir die Ausbiegung einer gespannten 
Haut unter Querbelastung iibereinstimmt; bei unregelmafigen Umrissen der Platte kann also 
die Lésung experimentell auf Grund des Membrangleichnisses gefunden werden. 


Bei polarer Symmetrie wird 


nN, =m, =, (3.1.5) 
also zufolge (2.4) 
dm 
Bored (3.1.6) 


Statt der Bruchlinien tritt bei Vollausnutzung ein differentiales Bruchnetz auf; die — 
Kriimmung der Biegeflache hat im vorliegenden Falle in allen Richtungen gleiches Vorzeichen. _ 


b) Ungleichsinnige Vollausnutzung. Vollausnutzung kann nun auch bestehen, 


wenn m, = — m, = mist; der Momentenkreis hat dann den Halbmesser m und das Maximum | 
++ m ist nurin den Hauptrichtungen erreicht. Die Trajektorien dieser Hauptrichtungen kénnen | 


dabei, je nach dem vorliegenden Fall, noch verschiedenen Verlauf haben; wir beschranken uns 
auf drei wichtige Sonderfalle. 


Ist die Trajektorienschar kartesisch, d.h. besteht sie aus zwei Geradenscharen und 
sind diese den Richtungen x, y parallel, so wird 


mM, = —My =m, My = 0 (3.2.1) 
und wegen (2.1) und (2.3) 
he Om ed om 
liane q=- aa (3.2.2) 


1S. FuBnote 3 von Seite 151. 
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und 


0 0? 
(aaa = a3) = ie (3.2.3) 
Kennt man eine Lisung von (3.1.4), so gewinnt man eine Lésung von (3.2.3) durch Ubergang 
von y auf y//—1. 


Kreuzt eine kartesische Trajektorienschar die Richtungen x, y unter 45°, so ist 


mM, =m, =0, Mey =m, (3.2.4) 
also 
Ee om __ Om 
| Ix ~ Oy ? a ae (3.2.5) 
und 
o2m 1 
Ge by 7 P (3.2.6) 
| Besteht also, dritter Sonderfall, polare Symmetrie, so ist 
m,=—m =m, (3.2.7) 
also 
tel d 
: rqg=m-—+ ar (rm). (3.2.8) 


Die Kriimmungen der Biegefliche haben hier in den beiden Hauptrichtungen entgegen- 
| gesetztes Vorzeichen. 


c)Superpositionsprinzip. Die Uberlagerung von Lisungen zwei verschiedener 
Lastfalle ist in der Plastostatik stets dann und nur dann streng zulassig 1, wenn zu beiden das 
gleiche Bruchbild gehért. Hier trifft dies dann 
zu, wenn die ganze Platte entweder gleichsinnig ye ip 
oder ungleichsinnig ausgenutzt ist oder wenn die, ie 
Grenze zwischen den diesbeziiglichen Gebieten fiir 
beide Lastfalle die gleiche ist. 


d) Vergleich mitdemelastischen 

Zustand. Wird ein Stabwerk so bemessen, 
_da8 unter der Héchstlast alle Querschnitte mit 
dem vollen FlieBmoment ausgenutzt sind, so tritt 
durch die Plastizierung keine Neuverteilung der 
Momente, sondern nur ein Spannungsausgleich 
innerhalb der einzelnen Querschnitte? ein; wenn 
dies auch fiir Platten zutrifft, so mite auch im Gr ty 
elastischen Zustand die Randspannung go, aller 
Elemente in Richtung x, y absolut gleich sein a 7 b 
| und (3. 1. 4) und (3. a. 3) miiBten auf geometrisch Abb. 1. Plattenelemente mit Spannungsresultanten, 
| vertragliche elastische Formanderungen a mit kartesischen Koordinaten, b bei polarer Symmetric. 
fihren. Sind o,, y,7 die Spannungen des Platten- 
'elementes Abb. la in beliebigem Abstand z von seiner Mittelebene und y die Querdehnungs- 
-zahl, so muBte also 3 
2 92 Orr 
Fra (Gn — 7.05) + 5 (8y —¥ 04) = 2 (1 CY) 35 dy 


) sein, wobei, je nachdem gleichsinniger oder ungleichsinniger Ausgleich vorliegt, 


22 


1 H. Craemer, Abh. Int. Ver. Briicken u. Hochb, 12 (1952), 5S. 79. 
2 S. FuBnote 1 von Seite 151. 
3S. Timoshenko, Theory of Elasticity, 5. 24, New York 1951. 
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wire. Dies in die vorige Gleichung eingesetzt, liefert aber 


0° GON Hi 
(zs ay)(z) = 
eine Bedingung, die nur in Ausnahmefallen mit (3.1.4) oder (3.2.3) zugleich bestehen kann. Im 


Gegensatz zum Stabwerk wird daher im allgemeinen nicht nur eine Spannungsumlagerung 
innerhalb der Querschnitte, sondern auch eine solche der Momente eintreten. 


4. Isotrope Platte, Anwendungen. a) Kreis- un d Ringplatten. Wir stellen einige 
Falle polarer Symmetrie voran; da hier q aus dem Gleichgewicht sofort bestimmbar ist, sind | 
die Lésungen besonders einfach. Eine Kreisplatte vom Halbmesser ry sei mit p = konst. be-} 
lastet, so daf | 


1 
a ene (4.1.1) | 
ist. Die Platte sei am Rande frei gestiitzt, demnach ist 


Th = OMe fure tear (4.1.2) | 


Abb, 2. Ringplatte am AuBenrande eingespannt. Abb. 3. Volle Kreisplatte am Au®enrande eingespannt, 

Die Senkungsflaiche im Bruch hat Kriimmungen tberall gleichen Vorzeichens, so daB gleich-| 
sinnige Vollausnutzung besteht. Aus (3.1.6) folgt dann 

m= —pr/44+C; 


die Konstante C bestimmt sich aus (4.1.2), und es wird 


m= f. (re — r?). (4.1.3) 

Fiir eine in Plattenmitte wirkende Einzellast P erhalt man auf ganz ahnliche Weise 
x 12 1 Dp . | 
m= >In. (4.1.4) 


Ist auBerdem eine Flichenlast p vorhanden, so folgt aus Abschnitt 3 c), daB die rechten Seiten; 
von (4.1.3) und (4.1.4) zu addieren sind; dasselbe erhalt man, wenn man q in (3.1.6) als Funktion) 
von p und P ausdriickt. 

Eine Kreisringplatte nach Abb. 2 sei am auferen Rande eingespannt und mit p gleich-| 
maBig belastet. Da die radiale Kriimmung negativ, die tangentiale positiv ist, kann nur} 
ungleichsinnige Vollausnutzung eintreten. Es ist | 


rq=5 P(e pce. 


also nach (3.2.8) 


C 
m == + + (2% — ry. 


Die Konstante C bestimmt sich daraus, da8 am inneren Rande m = 0 sein mu, so daB 


P Cn 
m= — me —2ce4 =) (4.1.5) 
wird. 


Ks sei nun eine volle Kreisplatte am Rande eingespannt (Abb. 3). Aus dem Bruchbild folgt 
daB fiir 0 <r Sc die Ausnutzung gleichsinnig, fiir ec < r < ry ungleichsinnig ist. Da fir r = ¢ 
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ee eS eee 


-raschen Abklingen der Exponentialfunktion Gebrauch machen kann; 


die radialen Momente m, nicht gleichzeitig positiv und negativ sein kinnen, muB 
Mia fur 7 = ¢ (4.1.6) 


sein; (4.1.1) gilt fiir beide Bereiche. 
Fir den inneren Bezirk ist somit (4.1.3) anwendbar, wenn man c statt ro einfiihrt: 


m= (2 — 7), UPS Se (4.1.7) 


wi 
Fur den auBeren Teil erhalt man aus (4.1.1) und (3.2.8) 


ee 7 
ie 8 
oder unter Beachtung von (4.1.6) 
4 
Vea (r 7) (<r<r,). (4.1.8) 


Hierbei ist die GréBe c frei wahlbar; diese Tatsache ist das Gegenstick zu der frei wahlbaren 
SchluBlinie der durchlaufenden Balken. 


b) Rechteckplatten. Eine rechteckige Platte nach Abb. 4 sei gleichmaBig mit p 
belastet und frei gestiitzt, so daB die Randbedingungen 
1 
ne = © fir 0 eary und im 2.0" fur y=ist (4.2.1) 


bestehen. Die Ausnutzung ist offenbar gleichsinnig. 
Fir die Lésung empfiehlt sich die Einfithrung der in der Ab- 
bildung dargestellten Hilfsordinaten y, und y,, damit man von dem 


sie lautet dann 


2 aS a7 sy 
Po {Ps Bae #4 Inf, _ exp ( nw y¥,/a + exp ( cen 


=e = a n3 1 + exp (— na b/a) 


J 
(4.2.2) 
wobei (— 1), anzeigt, daB die Reihe alterniert. Daf die hier maB- 
gebende Gleichung (3.1.4) erfillt ist, laBt sich nachweisen, wenn man 
dort p als Fourierreihe mit Periode 2 a entwickelt und dann (4.2.2) anwendet. Ebenso lassen 
sich die Randbedingungen (4.2.1) leicht bestatigen, wenn man fiir x, y, und y, die ent- 
sprechenden Werte in (4.2.2) einsetzt. Es laBt sich ferner zeigen, da fir a/b = 0, wie zu 
erwarten, der Verlauf von m in x parabolisch wird. 

Statt dessen kann man auch die Doppelreihe 


Abb. 4. Rechteckplatte. 


eet (—1)i(— 1), ab UO &% kay 


: cos cos 
ut Lk (b? i? a? i? a b 
i ere Ce ( i. ) 


(4.2.3) 


verwenden. Hier kann der Nachweis der Richtigkeit dadurch erbracht werden, daB man in 
(3.1.4) p durch eine Doppelreihe mit Perioden 2 a und 2 b ausdriickt und (4.2.3) einsetzt. Die 


_ Erfillung von (4.2.1) ist sofort ersichtlich. Fir x = y = 0 z. B. liefert die Auswertung 


m = 0,076 pab. (4.2.4) 


Steht eine Einzellast P in der Mitte der Platte, so ist die Lésung bei y= 0 unstetig; sie kann 


fir y >0 geschrieben werden 


has ie 1 ag Tut exe Bimiyle) = ep (= mm Bie 6) sexp. (= xt 9/4) (4.2.5) 
— n a , 1 + exp (—n a b/a) 
Sie befriedigt ersichtlich die Randbedingungen (4.2.1) an den Randern x = + a/2 und y, = 0. 


Sie befriedigt ferner, wie es sein soll, die Gleichung (3.1.4) mit p = 0. Bildet man schlieBlich 
nach (3.1.3) die Querkraft q, fiir einen Schnitt y = 0 dicht neben der Last, so muf diese 
Querkraft nach Gréf8e und Verteilung mit der halben Belastung P iibereinstimmen; daf dies 
tatsdchlich der Fall ist, kann durch Darstellung von P als Fourierreihe nachgewiesen werden, 

Wir behandeln weiter ein Beispiel fir ungleichsinnige Vollausnutzung. Eine quadratische 


3% 1 é as : 
Platte mit Seitenlange a und Koordinaten nach Abb. 4 sei in x = +4 frei gestiitzt; die 
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Rander y = + = a seien vollig frei und je durch eine Linienlast p = konst. belastet. Offen- 


sichtlich ist die Hauptkriimmung in Richtung x positiv, in Richtung y negativ. Es gelten die 
Randbedingungen (4.2.1) und auBerdem mubB. 


pe 1 
Gy ee Pe AOS SY ae (4.2.6) 
sein. Die Lisung 
m= oe cos “= * cos ase (4.2.7) 
1, 3,5 


geniigt, wie leicht nachzuweisen, der hier maBgebenden Gleichung (3.2.3) mit p = 0, sowie | 
den Bedingungen (4.2.1). Da8 sie auch (4.2.6) befriedigt, folgt, wenn man qy nach (3.1.3) bildet — 
und p als Reihe darstellt. 

SchlieBlich untersuchen wir noch eine Rechteckplatte a- b mit gleichmafiger Last p, die 
an zwei aneinanderstoBenden Seiten eingespannt, an den beiden anderen vollig frei ist; den 
Ursprung des kartesischen Koordinatensystems wahlen wir aber in der freien Ecke; es muB also 


m=0,q,.=0 bzw. qy=0 fir «=0 baw. y=—0 (4.2.8) 


sein. Die Betrachtung der Bruchverformung zeigt, daB hier der Fall von (3.2.6) vorliegt. Die 
Loésung 


m= — Spry (4.2.9) 


geniigt (3.2.6) und den Randbedingungen (4.2.8). Die Randtriger werden nur durch Quer- 
krafte 


ee bzw. gy = ——=ph 
und Schubmomente 
May = — > pay bzw May = — pha 


gema (3.2.5) und (3.2.4) beansprucht. 


Auch bei den vorangegangenen Beispielen lassen sich die Auflagerdriicke aus den Quer- 
kraften fir den betreffenden Rand ermitteln. Bei der nicht voll ausgenutzten Platte hingegen 
kénnen sie, da die Bruchstiicke elastisch bleiben, nur aus elastostatischen Betrachtungen ge- 
wonnen werden. 


5. FlieBzustand orthotroper Platten. Im Falle orthogonaler Anisotropie bestehen fiir jedes 
Element zwei aufeinander senkrechte Richtungen a und ), fiir die die Biegemomente m, oder 
m,, unter denen Plastizierung eintritt, absolute GréBt- und Kleinstwerte sind; wir bezeichnen 
sie als Orthotropiehauptrichtungen und HauptflieBmomente. Die 
GréBe von m, und m, ergibt sich wie in (3.1) und (3.2), jedoch mit verschiedenen Werten o, 
und f, fiir die beiden Richtungen. 

Die eben eingefithrten Begriffe haben mit den in Abschnitt 2 definierten Hauptrichtungen 1 
und 2 und Hauptmomenten m, und m, zunachst nichts zu tun, da diese von der Belastung des’ 
Tragwerks abhangen, jene aber reine Baustoffeigenschaften darstellen. Bei einer parallel zum 
Stamm geschnittenen Holzplatte z. B. verlaufen die Richtungen a und b parallel und senkrecht 
zu den Fasern, wahrend sie bei orthogonal bewehrtem Stahlbeton mit den Bewehrungsrich- 
tungen ubereinstimmen. Zur Vereinfachung werden wir aber im folgenden voraussetzen, da 
die Richtungen 1 und 2 mit jenen a und b zusammenfallen; Plastizierung tritt also ein, wenn 

entweder m, =m, oder m, = m, (5.1) 
ist. 

Die Orthotropiehauptrichtungen sind nun im allgemeinen Falle nicht fiir jedes Element 
die gleichen, sondern bilden ein rechtwinkliges Netz von Orthotropietrajektorien, 
wie man z. B. am Faserverlauf in der Nahe eingewachsener Aste erkennt. Wir setzen hier 
entweder kartesische oder polare Orthotropie voraus, so da die Trajektorien entweder zwei 
Geradenscharen oder eine Schar von Radien mit den zugehérigen Kreisen bilden. 
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a) Gleichsinnige Vollausnutzung. In Erweiterung von Abschnitt 3a) ver- 
stehen wir unter gleichsinniger Vollausnutzung, daB simtliche Elemente mit 


m =mMm=m, m=—=m=um (shui) 


} ausgenutzt sind; das ,Orthotropieverhdltnis“ u kann hierbei beliebige positive 
| Werte annehmen und sei vorlaufig als Konstante angenommen. Man sieht, daB, abweichend 
/ vom Bisherigen, hier kein Ausgleich der Spannungen im ganzen Tragwerk vorliegt, sondern 
nur ein solcher der Ausnutzungen. 

Ist die Orthotropie kartesisch, so lassen wir die Koordinaten mit ihren Hauptrichtungen 
zusammenfallen, so daB 


m, =m, M=p"Um, my = 0 (Saka) 
wird; (2.1) und (2.3) gehen dann iiber in 
am ye ane (5.1.3) 


und 
e 0? 
(a +- Lu a8) nas-— P (5.1.4) 


als Erweiterung von (3.1.3) und (3.1.4). 

Das Membrangleichnis kann auch hier herangezogen werden, wenn man eine in den beiden 
Richtungen verschieden stark gespannte Haut verwendet. Lésungen von (3.1.4) fihren auf 
Lésungen von (5.1.4), wenn man von y auf y/Vu ubergeht 1. 

Bei polarer Symmetrie wird 

m =m, m,=—"wm (5:15) 
und wegen (2.4) 


rq=(l—p)m4+r@ (5.1.6) 


als Erweiterung von (3.1.6). 


b) Ungleichsinnige Vollausnutzung. Man sieht, da®B die Angaben des 
vorigen Abschnitts giiltig bleiben, wenn man dort fu mit umgekehrtem Vorzeichen einsetzt. 


6. Orthotrope Platten, Anwendungen. a) Kreisplatten. Es sei wieder die in Ab- 
schnitt 4 a) untersuchte frei gestiitzte Kreisplatte mit Belastung nach (4.1.1) betrachtet; die 
Ausnutzung ist gleichsinnig, und es wird 


m = pt ae (zy (6.1.1) 


Die Richtigkeit laBt sich durch Einsetzen in (5.1.6) beweisen, die Erfiillung von (4.1.2) ist sofort 
_ ersichtlich. In der Plattenmitte, r = 0, erhalt man m = 0 bzw. m = oo, je nachdem y z List; 
fir 1 = 1 fihrt die Lésung auf (4.1.3) zuriick. 

Wir behandeln weiter die an Hand von Abb. 3 untersuchte eingespannte Platte; fiir den 
_inneren Teil 0 <r <c gilt (6.1.1) mit c an Stelle von ry; auBen (¢c Sr X 1) ist der Ausgleich 
ungleichsinnig, und es ist Absatz 5 b) zu beachten. Die Lésung lautet 


as ike (ean earls (6.1.2) 


m= 2 (3 +4) \\r 


Setzt man dies und (4.1.1) in (5.1.6) mit umgekehrtem Vorzeichen von y ein, so iberzeugt man 
sich leicht von der Richtigkeit; fiir r = c wird m = 0, wie es (4.1.6) verlangt. 


b) Rechteckplatten. Wir beschranken uns auf die Erweiterung von (4.2.2) und 
(4.2.3) auf orthotropen Baustoff. Es wird 


m cos = 
me eam a 1 + exp (— na b/a)u) 


_ Apa (—1)n ala (—nayi/a Vu) + exp (—nza y,/a |) (6.2.1) 
13,5 


1 Vel. auch die Affinitatssatze in der Arbeit von Johansen, FuBnote 3 von S. 151. 
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oder 


16 pab (—1)i(—1), ab tm % kay 2.2 
DFU Ph eee Bes neg A (6.2.2) 


ge. 1y 3,5 
Durch Vergleich mit (4.2.1) und (5.1.4) iberzeugt man sich leicht von der Richtigkeit. 


7. Wechselndes Orthotrepieverhiltnis. In Abschnitt 5 hatten wir vorausgesetzt, daf} das 
Verhiltnis w eine Konstante sei. Im Stahlbeton nun entsteht die Orthotropie gemaB (3.2) bei 
im allgemeinen gleichbleibendem h durch verschiedene Dichte f, der Bewehrungen in den 
beiden Richtungen; j kann also an verschiedenen Stellen der Platte verschieden sein. Damit 
reduziert sich aber die Bemessungsaufgabe auf die Erfillung der Gleichgewichtsbedingungen 
(2.3) mit m,, = 0, d.h. von 


O?mx | O?my has i il 
Oa? 1 Oye ee P Se 
oder von (2.4). 


In der Praxis teilt man eine gleichmaBige Belastung p meist in zwei unveranderliche Teil- 
betrage p, + py auf, so daB 


O?mx 02m 
Oe eet ay? Sie 5S (2) 
also 
a 4 x? pb? 4 y? 
es, =F (1 sf z) my= (1 2 (7.3) 


wird. Das Verhiltnis p,/py wird dabei so gewahlt, daB man dem elastischen Verhalten még- 
lichst nahe kommt. Es ]4Bt sich aber leicht zeigen, da® dies mit einem solchen konstanten Ver- 


: ; : a 1 
haltnis nur unvollkommen erreicht werden kann, da8B vielmehr z. B. an den Randern « = ++ 5 


der Anteil py verschwinden miiBte. Dieses weithin geiibte Verfahren, dessen Begrindung 
immer mit Uberlegungen nur aus der Elastostatik versucht wird, findet somit erst auf Grund 
der Plastostatik — und nur innerhalb ihrer Giltigkeitsgrenzen — seine volle Rechtfertigung. 


8. Zusammenfassung. Aus der Bedingung, daf alle Elemente der Platte in zwei aufeinander 
senkrechten Richtungen mit dem vollplastischen Moment gleichen oder entgegengesetzten 
Vorzeichens beansprucht werden, werden die Differentialgleichungen fiir die Bemessung iso- 
troper und orthotroper Platten in den wichtigsten Sonderfallen entwickelt. Die Anwendung 
wird an Beispielen erlautert. Ferner werden Angaben iiber die Bemessung von Stahlbeton- 
platten gemacht. 


(Eingegangen am 19. April 1954.) 


Anschrift des Verfassers: Prof. Dr. H. Craemer, Lahore, Pakistan. 
Moghalpura College of Engineering & Technology. 
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Beitrag zum Problem der Wirmespannungen in Scheiben. 


Von E. Tremmel. 


1. Einleitung. Unter dem Einflu8 stationarer ebener Temperaturfelder T(x, y) bleiben die 
Querschnitte zylinder- oder prismenférmiger Kérper oder diinne Scheiben, da die Vertriglich- 
keitsbedingungen hier erfillt sind, spannungsfrei, wenn sie keine Kurven enthalten, die eine 
Warmequelle umschlieSen!. Denn nur dann kénnen die aus den Temperaturverzerrungen allein 


| folgenden Verschiebungen eindeutige Ortsfunktionen sein. Ist das nicht der Fall, so werden die 


mehrdeutigen FormanderungsgréBen durch Uberlagerung eines zusitzlichen, einem Selbst- 
spannungszustand zugeordneten Verschiebungsfeldes getilgt. Zu seiner Festlegung ist daher 
eine den Randbedingungen angepaBte Airysche Spannungsfunktion so aufzubauen, dab die aus 
ihr abgeleiteten Verschiebungsgréfen bestimmte, die mehrdeutigen Temperaturverschiebungen 
kompensierende Funktionen enthalten. 

In der vorliegenden Arbeit soll das Problem zunachst allgemein behandelt werden; die Auf- 
gabe, eine Spannungsfunktion aus den berechneten Temperaturverschiebungen zu bestimmen, 
fiihrt hier auf eine Poissonsche Differentialgleichung, deren Partikularintegral aus einem in 
komplexer Form ausgedriickten Ansatz berechnet wird. Da die Vertraglichkeitsbedingungen 
im Falle eines stationaéren Warmeflusses, wie schon erwahnt, von vorneherein erfillt sind, bleibt 
die auf Temperaturwirkungen erweiterte Differentialgleichung der Airyschen Spannungs- 
funktion homogen. Die gesuchte, die Randbedingungen befriedigende endgiiltige Lésung wird 
daher durch Uberlagerung jenes Partikularintegrals mit anderen passend gewahlten Integralen 
der Bipontentialgleichung als Linearkombination biharmonischer Funktionen erhalten. In 
einem Beispiel wird die praktische Durchfithrung der Berechnung an einem von exzentrisch 
liegenden Kreiszylindern begrenzten Rohr gezeigt werden. 


2. Berechnung der Verschiebungsgréfen. Fiir einen prismatischen oder zylindrischen 
Korper, dessen Endflachen (z = -+ a) unverschieblich festgehalten seien, ist ein in der Lings- 
richtung z konstantes stationares Temperaturfeld T(x, y), das also der Laplaceschen Differen- 
tialgleichung 

See ae 
geniigt, durch die Randbedingungen festgelegt. Wegen (1) kann unter Verwendung der 
Cauchy-Riemannschen Differentialgleichungen aus T die Funktion S(x, y) so bestimmt werden, 
dah 

(= tees (2) 
eine analytische Funktion der komplexen Variablen t= x -+ 1 y wird. Bedeuteny die Poisson- 
sche Zahl und ap den Warmeausdehnungskoeffizienten, dann sind die aus den ‘Temperatur- 
dehnungen allein folgenden Verschiebungen u, v in x- und y-Richtung in komplexer Darstellung 
durch 

w(t) =utiv=(14+)orf{Qd (3) 

gegeben. Denn unter der Voraussetzung spannungsfreier Querschnitte gilt fiir die Verzerrungs- 
komponenten ! 


Csg yy Or Exy= 0; 
werden sie durch die Ableitungen der Verschiebungen ausgedriickt: 


Ou Ov 1 & = 


= = ——e | ex 
bo Fie Ng 2 kd dy’ aes 2\dy | ax 


so ergeben sich die Cauchy-Riemannschen Differentialgleichungen 
du Ov ou dv 


Ox dy’ Oy On 


1 E. Melan, Temperaturverteilungen ohne Warmespannungen, Osterr. Ing.-Arch. 6 (1951), 5. ihe 
E. Melan, Osterr. Ing.-Arch, 4 (1950), $.153; E. Melan u. H. Parkus, Warmespannungen, 5.9 ff., Wien 
1953. 
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d. h. auch u und v sind analytische Funktionen der komplexen Variablen t: 


utiv=a(t). 
Mit 
, du, . dv - Ou 
WAU eae e ‘On dy * dy 
folgt 
, Ov Ou 
Dyas ets ofa igs =) 


wobei der Realteil dieser Funktion durch 
Ou Ov 
eer =2(1+»)arT, 
der Imaginarteil daher unter Verwendung der oben berechneten Funktion S durch 


Ov Ou 


Seep mee +) or S 


gegeben ist. Aus der ee ee 
= fw'(t) dt=(1 +) ar {(T +78) de 


ergibt sich schlieBlich die a ie ene Gleichung (3). 
Der Korper wird bis auf die Spannungen in der 2-Richtung , fiir die 


Op eadG (epiant (4) 


gilt, spannungsfrei sein, wenn die Verschiebungen (2) eindeutig sind. Das ist nur dann der Fall, 


wenn f 0 ree 5) 


wird, wenn also, wie aus dem Gaufschen Integralsatz hervorgeht, in dem Querschnitt keine 
Kurven existieren, die eine Warmequelle umschlieBSen. 

Wir teilen das gegebene Temperaturfeld in einem Teil T,(x, y), fur den das Integral f Qy dt 
verschwindet, der also eindeutige Verschiebungen bewirkt, und in einem Teil T,(x, y), fir den 


f Q, dt + 0 
wird, so dafs die von diesem Teil gelieferten Verschiebungen mehrdeutig werden. Es sei 
aH yl bn (6) 
die Verschiebungen zufolge dieses Temperaturfeldes sind nach (3) durch 
u+iv= (1+ )ar{(Q) + Q,) dt (7) 
gegeben und setzen sich aus dem eindeutigen Anteil 
Uy + ivp= (1 +) ar f Qy dt (8a) 
und dem mehrdeutigen Anteil 
uy +iv, = (1 +7) or | Qy dt (8b) 


zusammen. Die mehrdeutigen Verschiebungen werden nun durch Uberlagerung eines zusatz- 
lichen Verschiebungsfeldes getilgt, wobei der diesem Feld entsprechende Spannungszustand 
den Randbedingungen geniigen muB. 


3. Einfithrung der Airyschen Spannungsfunktion: Wir denken uns den Querschnitt so auf- 
geschlitzt, daB samtliche, die Warmequelle umschlieBenden Kurven durchtrennt sind und 
haben nun eine Airysche Spannungsfunktion y zu bestimmen, welche die im geschlossenen 
Querschnitt mehrdeutigen Temperaturverschiebungen kompensiert und die Randbedingungen 
erfillt. Es handelt sich also um ein gemischtes Randwertproblem; die Ausdriicke fiir die Ver- 
schiebungen miissen bestimmte vorgegebene Funktionen enthalten, wahrend die Randspan- 
nungen verschwinden sollen. 


Die Spannungsfunktion wird nun in der Form 


Y= + Ye (9) 
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angesetzt, wobei das zu y, gehérige Verschiebungsfeld die mehrdeutigen Temperaturver- 
schiebungen aufheben soll, wahrend die mit y, verkniipften Verschiebungen eindeutig sein 
miissen. 

Fir die aus einer Spannungsfunktion py abgeleiteten Verschiebungen gilt im Falle eines 
ebenen Verzerrungszustandes bekanntlich !: 


P l—y : 1 /oyp -0 
ut+iv=r— [¢ | i m) dt — 55 (5e +854), (10) 


wobei 
l= Ayp=o,.+0, (11) 


und m mit / durch die Cauchy-Riemannschen Differentialgleichungen verknipft ist, so daB 
| + im eine analytische Funktion der komplexen Vaiablen t = x + 7 y ist. 

Da das Auftreten mehrdeutiger Verschiebungen bei eindeutiger Spannungsfunktion an den 
Integralterm in (10) gebunden ist, mu8, um das Verschwinden der mehrdeutigen Temperatur- 
verschiebungen zu erzwingen, wie aus dem Vergleich von (8b) mit (10) hervorgeht, gelten 


(1 +9)ar [ Qdt—= — ATH [tim a (12) 
oder nach Differentation und Aufspaltung in Real- und Imaginarteil 
(1 + 9) arT, = SA. (13) 
woraus fiir y, die Poissonsche Gleichung 
Ay, = — 26 * oof, (14) 
folgt. Setzen wir zur Vereinfachung noch 
= 265 aH (15) 
dann geht (14) tiber in 
JA Speed Ie (14a) 


und es besteht nun die Aufgabe, ein Partikularintegral von (14a) zu bestimmen. 
Da es sich voraussetzungsgemaB um einen stationdren Temperaturzustand mit AT = 0 
handelt, geniigt y, der homogenen Bipotentialgleichung 


AAy, = 0; (A) 
Die Lésung von (14a) 148t sich daher durch Integrale von (A) darstellen. Das gleiche gilt in 
diesem Falle natiirlich auch fiir die vollstandige Lésung y des Randwertproblems. 


4, Partikularintegral der Differentialgleichung Ay, = x T,. Wird in Anlehnung an bekannte 
Verfahren 2 die reelle Funktion y, in Abhangigkeit von den zueinander konjugiert komplexen 
Minimalkoordinaten 

t=x+1y t= x — Ly 
dargestellt: y F, 
vi=tf) +ef(t), CO) 
wobei bekanntlich “ponies: 
f)=fO cD 
die konjugiert komplexe Funktion von f(t) ist, und beachtet man, daB fir den in Minimal- 
koordinaten angeschriebenen Laplaceschen Operator gilt 


so folgt aus (16) 
ty, = [F') +FQ)]s as) 


1 C. B. Biezeno u. R. Grammel, Technische Dynamik Bd. 1, Seite 108 ff., 2. Aufl., Berlin 1953. 
2.N, MuSelisvili, Z. angew. Math. Mech. 13 (1933), S. 264. 
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wegen (17) aber bestehen die Beziehungen 


4f(i)}=ptiq, 4f(t)=p—igq, 
so daB 


Alf’) — fie f= 2iq (20) 


als Imaginarteil der komplexen Funktion 


als Realteil und 


2(p +14) = 8f'(t) (21) 
betrachtet werden kénnen. Mit den Bezeichnungen von (10) bzw. (11) gilt ferner 
8f'(t)=l+im. (21a) 
Vergleicht man nun (14a) unter Beriicksichtigung von (19) mit (21), so folgt mit 
Q= 7, +715, 


die Gleichung 
xQ, = 8f"(t), (22) 


aus der sich durch Integration die Funktion 


fo=F [a dt (23) 


ergibt, die voraussetzungsgemaB mehrdeutig ist. 
Formen wir die rechte Seite durch Teilintegration um: 


f() = 3 (9 gen | ae i) (23a) 


und bilden aus dem ersten Klammerterm mit f, ,o(t) = - Q,t die Spannungsfunktion entsprechend 


(16): 
Vi0= = til i(t) + Q,(¢)| ’ 


die mit 


tt =r? 
und 
Q,(t) + Q(t) = 2 Re Q,(t) = 2 T, 


iibergeht in 
Hw 
ets) eas r’T,, 


so ist y,,) nur vom Quadrat des Radiusvektors und der gegebenen Temperaturverteilung ab- 
hangig und daher offenbar eindeutig. 

Der aus dem zweiten Glied folgende Verschiebungsanteil kann je nach der Querschnittsform 
und den gegebenen Randtemperaturen ein- oder mehrdeutig sein. Im ersten Fall, der z. B. beim 
Kreisring mit rotationssymmetrischer Verteilung der Randtemperaturen auftritt, ist er ein- 
deutig und daher ohne Einflu® auf den Spannungszustand. Im zweiten Fall ist eine eindeutige 


Spannungsfunktion y,, zu bestimmen, welche die nunmehr an die Stelle von (22) tretende Be- 
dingung 


pf dt +3! dQ, = 0 (22a) 


erfillt. Diese Forderung, die fiir jede Warmequelle, d. h. bei mehrfach zusammenhangenden Be- 
reichen fiir jeden Lochrand aufzustellen ist, besagt: der Zuwachs, den die mehrdeutige Ver- 
schiebungsfunktion u,, +7v,, nach einem Umlauf entlang einer das Loch umschlieBenden Kurve 
erfahrt, muB den entsprechenden, von den Temperaturverzerrungen abhangigen Zuwachs, die 
»,Restverschiebung”, aufheben. 

In praktischen Fallen kann aus der Beschaffenheit von f t dQ, unmittelbar auf die zur Kom- 
pensation der Restverschiebung erforderliche Funktion u,; + i v;, und damit auf f,,(t) geschlossen 
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werden; die Aufgabe besteht demnach nur darin, die den vorgegebenen Verschiebungsanteilen 
zugeordnete eindeutige Spannungsfunktion zu finden, die im Verein mit Y1,9 und wy, den Rand- 
bedingungen geniigt. Unter Verwendung von (21a) erhalten wir aus (10) 


G ; 
as) ee WE) (uy + if). (23b) 
Fihrt man (23b) in (16) ein, so folgt 
G = : : 
PA | A) [i (ur, + ivf) +t (us, —7 vs) (16a) 


und schlieBlich in reeller Darstellung 


G 
Yu = ia Why EY V4) « (16b) 


Die so gewonnene Spannungsfunktion wird aber die Eindeutigkeitsbedingung nicht erfiillen, da 
eine der beiden Verschiebungskomponenten mehrdeutig ist. Die gleichen Verschiebungsanteile 
Uj,» Vj, aber werden, wie sich mit (16a) bzw. (16b) leicht nachpriifen lat, auch von den Spannungs- 
funktionen 
G 
Yuu) = + =r u* (16c) 


und 


G 
Pi?) te (16d) 


geliefert, unter denen nun die im speziellen Fall eindeutige Funktion auszuwahlen ist. Aus der 
Uberlagerung von y,,) und y,, folgt schlieBlich die vollstandige Spannungsfunktion y,. 

Die Spannungsfunktion y, ist nun aus Linearkombinationen von bekannten, dem speziellen 
Problem angepaBten Partikularlésungen der Bipotentialgleichung so aufzubauen, daB mit der 
Funktion 

Pa ae) 
die Randbedingungen erfillt werden. Damit ist der Spannungszustand des Kérpers zufolge 
des stationaren Warmeflusses bestimmt. 


Die zur (x, y)-Ebene parallelen Spannungen sind mit der Airyschen Spannungsfunktion 
bekanntlich durch die Beziehung 


ay 
Oik = =e Onn — ik Ay), Xin Xk = 


verkniipft, wobei das Kroneckersymbol wie tiblich mit 6;; bezeichnet ist. 


Fiir die im ebenen Verzerrungszustand in der z-Richtung auftretende Spannung gilt 
G,.=V Ay. (a) 


Aus der Uberlagerung dieser Spannung mit der nach (4) berechneten, dem geschlitzten Quer- 
schnitt zugeordneten Spannung o, erhalt man schlieBlich die endgiltige Spannung in der 
Langsrichtung zu 


O6,= —2G6(1+»)o7T +rAy. (24) 


5. Von exzentrisch liegenden Kreiszylindern begrenztes Rohr; Temperaturspannungen. Als 
Beispiel soll hier der durch einen stationéaren Warmeflu8 hervorgerufene Spannungszustand 
sowie das zugehérige Verschiebungsfeld in einem zylindrischen Rohr untersucht werden, dessen 
Querschnitt von zwei exzentrisch zueinander liegenden Kreisen begrenzt sei. Unter der Vor- 
aussetzung, da die Temperaturverteilung in der Laingsrichtung, der 2-Achse unverdnderlich 
sei, und da8 die Endflachen unverschieblich festgehalten sind, werden die Querschnitte einem 
ebenen Verzerrungszustand unterworfen sein. Fiir die Berechnung sind im vorliegenden Falle 
zweckmabig die durch die komplexe Transformation 1 


tox +iy=ibGtg Sih (Ging — isiny) (25) 


1 Siehe etwa W. Miiller, Ing.-Arch. 13 (1942), S. 37ff. 
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baw. plied gah (26) 


t—ib 


pete a Zot (27) 
Co}  — cos 


bestimmten Bipolarkoordinaten zu verwenden. 
Den Koordinaten £ = konst. enstprechen bekanntlich Kreise mit den (kartesischen) 


Mittelpunktskoordinaten (0, z)-) und den Radien rz, wobei 2; und r; durch 
y %ye= 1 big €, 
b (28a) 


' Ciné 
gegeben sind, und den Koordi- 
naten 7= konst. die durch die 
Pole Q,(0, +6), Q,(0, — b) gehen- 
den Kreise mit den (kartesischen) 
Mittelpunktskoordinaten (z,,, 0) 
und den Radien r,, wobei hier 


n? 


fy, — octgy, | 


ob (28b) 
See sin 4 (| 


eA 


jee gilt. 

ost N= & Mit den Beziehungen (28a) 
Be a kénnen bei gegebenen Radien R,, 
> R, der Randkreise und der Ex- 


zentrizitat e die diesen Radien 
entsprechenden Koordinaten a, f 
und die Konstante b berechnet 
werden (Abb. 1). 

Sei nun die Anderung gegen- 
iiber der zum spannungslosen 
Ausgangszustand gehérigen 
Temperatur am duBeren Rande 
« mit T,, die am inneren Rande fp 
mit T, bezeichnet, dann wird 
die Laplacesche Differential- 
gleichung AT = 0 durch denin & 


linearen Ansatz 


dStg P 


7 


Abb. 1. Bipolarkoordinaten, geometrische Daten. Ih a Co + Cy é (29) 


befriedigt 1 und die Konstanten errechnen sich aus den Randbedingungen zu 


Tx B —T T3 —T 
cy = EE Beret (30) 


Mehrdeutige Verschiebungen werden hier offensichtlich nur im Zusammenhang mit dem zweiten 
Term von (29) auftreten. Wir setzen demnach T, = cy und T, =c,& Mit T, =c,& ergibt 
sich nach (14) 


il 
Ay, = —27 "Coreg é (31) 
oder 
Ay, =xo€, (31a) 
wenn zur Vereinfachung der Koeffizient von c, € entsprechend (15) 
1 
—2 i = ; Corp=x 


1H. S. Carslaw, Conduction of heat in solids, S. 98, New York 1945. 
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gesetzt wird. Mit 
Ty =¢,é 


folgt aus den Cauchy-Riemannschen Differentialgleichungen bis auf eine unwesentliche Kon- 
stante 


Sy] = 7). 
und wir erhalten 
Qi = ¢ (§ +77) 
oder nach Einsetzen von (26) 
Q=%¢ 


_ und damit nach (23) 


f) = op oat. (32) 
Wird di durch Bildung der Ableitung von (25) in der Verdnderlichen¢ ausgedriickt: 


aoe re 
so geht (32) tiber in 
fi) = — mal Teo dt, (32a) 
und wir erhalten aus dem Integral 
fi) = +2 Sib ( Gtg = —2In Gin a (33) 


und unter Beriicksichtigung von (25) sowie nach Umformen und Aufspaltung des zweiten 
in der Klammer stehenden Terms: 


: a8 Street oj E cos —1 F Co}  — cos 7 
2ibIn Gin = b arc cos beens, -iblIn 5 ; 
fie (ce Bie gil (33a) 
wenn 
CofE cos 7 —1 
arc cos 


oj § — cos 4 aay 


gesetzt wird, wobei m das auf den Kreisen = konst. von der positiven Y-Achse (7 = 0) aus 
nach rechts gezahlte Azimut bedeutet, wie unter Verwendung von (25) und (28a, b) leicht nach- 
zuweisen ist. 

Fir die gesuchte Spannungsfunktion yp}. ergibt sich nun nach (16) in komplexer Darstellung 


Yio= 4 “l(c | Att —210(0m Sin S — tin gin 4)) (34) 


Setzt man hier C, gee § —in,tund t nach (25), (26) unter Einfihrung von h nach (27) ein, so 
folgt in reeller Form 


Yro= to le £ (Gin’ € + sin’ y) + bh (p sin 7 — Gin é In Snare | (34a) 


und nach weiterer Umformung 


vio= +O AE (Cos + cosn) +p sinn — in En 
Mit f(t) erhalt man ferner unter Beriicksichtigung von (8b) und (23a) 
u+iv = 2£0 + 9) or 
und damit fiir die von Temperaturverzerrungen allein bewirkten Verschiebungen 
m= + (1+ rare (x§—yy +b), (35a) 


v= + (1 + a)are (#7 ee |: (35b) 
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Mit den beiden aus (€ -+-C) tt folgenden Gliedern der Spannungsfunktion *,, kénnen die ¢ t ent- 
sprechenden mehrdeutigen Verschiebungsanteile in (35a) und (35b) identisch zum Verschwinden 
gebracht werden. 

Dagegen ist der mit g verkniipfte Ausdruck in (34a), (34b) offensichtlich mehrdeutig; damit 
wirde auch der vom zweiten Klammerterm in (10) herrithrende Verschiebungsanteil 


Leave | 2 0v Van. 1 op 
x6 (3s | Pie ce (i) 


mehrdeutig werden, was unseren durch (12) ausgedriickten Voraussetzungen widersprache. Zur 
Kompensation der aus (35a) folgenden mehrdeutigen Verschiebung 


i + tl + vr) are, by 
kommt daher das Glied sin 7, der Spannungsfunktion pf, ) nicht in Betracht. 
Nach Weglassen dieses Gliedes sowie des mit ihm gemafs (33) verbundenen Terms 


ey i ee 


5 erhalten wir fiir die Spannungsfunktion y,,, den Ausdruck 


Yio = + “Gb hE Gof € + cosy). (34) 


Die durch das explizite Auftreten der Funktion gy verursachte Mehrdeutigkeit der Verschie- 
bungen kann aber in folgender Weise eliminiert werden: 

Nach einem Umlauf haben die Verschiebungen u* um den Betrag (1 + »)arpce,b2z7 zu- 
genommen, einen hierzu proportionalen Zuwachs aber erfahrt eine Verschiebungsfunktion, deren 
Komponente uf mit den aus (10) folgenden Koeffizienten durch 


| — 1 
ut — D = Cy OE (é, 7) (36) 


gegeben ist, wobei C eine reelle Konstante und g(&, 7) eine eindeutige Funktion, die also keinen 
Beitrag zu dem Zuwachs liefert, bedeuten. Bestimmt man nun in (36) die Konstante C so, daB 
bei der Uberlagerung u* + u* 


(1 +») ore 620 +=" Con =0 


und damit 


I+» 
| C= —27— ore, ©b (37) 
wird, der Zuwachs also verschwindet, dann ist u* +- uf eine eindeutige Funktion, da die Haupt- 
werte von 

Cof € cos 7 —1 


Leer a end oj § — cos n 


mit den Werten von 7 in eindeutiger Weise verkniipft sind, demnach 


1l—y 
u; x6 Cy 


der eindeutigen Funktion ~ — 7 proportional ist. 

Um die zur Verschiebungskomponente uf gehérige eindeutige Spannungsfunktion y* zu er- 
halten, fiihren wir in (36) den oben bestimmten Wert von C ein; damit folgt fiir den maBgebenden 
ersten Term von (36), der mit uj, bezeichnet werden mége: 


ui, =—( 4+ varegbn; (38) | 


und aus den Cauchy-Riemannschen Differentialgleichungen bis auf eine unwesentliche Kon- 
stante fiir die entsprechende, zur Y-Achse parallele Komponente 


vi, =+(l+rarqbé, . (39) | 
die offensichtlich eindeutig ist. Die aus ihr gemaf (16d) gebildete Spannungsfunktion lautet nun | 
il 
Yu= +722 Corgbyé (40) 


| 
| 


XXIII. Band 1955 Tremmel: Beitrag zum Problem der Warmespannungen in Scheiben. 167 


oder unter Einfiihrung von (15) und (25) 
1 A 
Fir die Spannungsfunktion y, = yo + yy; erhalt man nun nach Einsetzen von (34c) und (41) 


yr = + Pb hs Coj £ + sin») 


% Cy 


=~, bbe ems (42) 


y 


| und mit 

mec; Lt py - 

.° 2 (1—») Gare, 6 = K 

(43) fe 
schlieBlich ne 


yi = h K [§ (Coj & + cos 7) 
—2&OGiné]. (42a) 
Die Spannungsfunktion yw, wird 
nun in der Form 


Y,=h{(A + C@oj 2é + D Gin 2 €) 1% Tyé 
x cosy + BE (Coj € — cos n) Pe 
-|- H Go} ra (44) Tn 


-angeschrieben, wobei die einzelnen 
Terme den fiir Bipolarkoordinaten 
entwickelten Partikularlésungen 1 
der Bipotentialgleichung entspre- 
chen. Mehrdeutige Verschiebungen 
treten zufolge yw, nicht auf, da mit 
Ay, = 1, offenbar f (I, +im,) dt=0 
wird. Aus der Uberlagerung von 
(42) mit (44) erhalten wir die voll- 
standige Spannungsfunktion: 


a h {[4 te (kK Fai B) g ae CC} 2 3 Abb. 2. Positive Spannungen und Verschiebungen 
+ Dein 2&] cos7 —- [H+ (K +B) &] im Bipolarkoordinatensystem. 
x Cof £—2 KE Ging}. (45) 


Fiir die Spannungen (Abb. 2) folgen aus (45) nachstehende Ausdriicke: 
op=— {A+ H+C+2KE + Gin [2 CSinE +2 DGols —(K + B)Cofé + 2K Ging] 
—[(K — B) + 2C Gin 2£ + 2 DOof 2 £)] Gin é cos 7} , 
lo, + (44H+ C+2KE+4+ Giné [2 C Giné+- 2 DCoj E+ (K+ B)Coj §—2 K(1 + Goj?§)] 
— 4 (CGoj 2& + D Sin 2 €) cos? n + [2 C (2 Cof 2E Cof & — Sin 2E Sin £) 
+2D (2 Gin 2 &Coj é — Coj 2 Gin) — (3 K + B) Gin —2 K Goj £] cos 7} , 
ian = (Coj é —cos7)[K —B+2(C Gin2& + D@oj 2 &)] sinn. 


(46) 


) Werden die Randbedingungen 
I 


= {5| Of = 1, = 9 


‘in (46) eingefiihrt, so erhalt man, da sowohl der von 7 freie, wie auch der mit cos 7 behaftete 
_Ausdruck fiir « und f verschwinden miissen, die erforderlichen vier Gleichungen fiir die unbe- 
kanten Konstanten 


Am. B,, Ci: 


| 1 W. Miiller, a.a.O. 5S. 39. 
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A und H sind hierbei nur als eine Konstante zu werten, da sie in den Ausdriicken fiir die Span- 
nungen, wie aus (46) ersichtlich, nur in der Kombination A + H auftreten; ferner in den 
durch (10) gegebenen VerschiebungsgréBen fir y = Ah cos 1) bzw. y= H h Coj &€ die Terme 
dp/dx baw. dp/dy verschwinden, so da die Verschiebungen nur von den Spannungssummen 
og +o, abhangen und daher ebenfalls den gemeinsamen Faktor A + H haben missen. Di 
Lésungen des Gleichungssystems lauten mit den neuen Konstanten K—B=B' und K+ B=F 


A+ Ha — yy PPP O* [ein 2 B — Eg (6 —a)] | 
+ (p +o) (N—D) + F(Gin26 + Gina) +(V—1) Cof2p—y], (Ata) 

BY = + of 26 —Gojf 2a +20), (47b) | 

F = — ~* of 28 —Gof 20 +21), (47c) 

C = + pS 5 Gof 28 — Gof 20 +2 N) ogre : (474) 
jas ny Gof 28 —Cof 20 +2.N) gehen (47e) 


wobei zur Abkirzung gesetzt ist 


1 : : 
L=f$—« —%g(B—a), N=BP—«1 —> (Gin 28 — Sin 2a). (48) | 
Setzt man fiir die Konstanten unter Einfiihrung der weiteren Abkiirzungen 
A+H Boe Fr 


in (46) ein, so erhalt man fiir die Spannungen 


m= Z(2etD +N Cutt (Po ay ea eel eee 


Coj (8 —4) 
—_2N, f Se = 4 Gin £ cos ih (50a) 
Ss 5 (28 fey aN Ga *) 4.2, Gin 2& — 2 (Goj2é +1) 
iy os ee 
a 2 Nelzen@—ay 8 Gin (6 + « —£) + Gin (f +a—38))—Giné) 
Sh (pan | cos 7 | , (50b) | 
tq = 2 Ny i = ead aioe | (Cof € — cos») in 7. (50c)_ 


Mit der Erfillung der auf die Normalspannungen g; beziiglichen Randbedingungen ver- | 
schwinden zugleich auch, wie aus dem Vergleich von (50a) mit (50c) hervorgeht, die Schub- | 
spannungen an den Randern a und f. 

Damit sind die in der (x, y) Ebene liegenden Spannungen bestimmt; sie sind, wie aus (30b) 
zu ersehen ist, dem Temperaturgefalle T;— T, proportional, waihrend die durch (24) gegebene 
Spannung o, auch von der gesamten Temperaturanderung abhiangige Glieder enthalt. 

Das hier abgeleitete Spannungsfeld kann man sich auch durch ein senkrecht zur Symmetrie- 
achse wirkendes Zug- oder Druckkraftepaar 1 hervorgerufen denken; der Rohrquerschnitt wird 
in 7 = 0 oder 7 =z geschlitzt angenommen. GréBe und Abstand a dieser Doppelkraft von der 
X-Achse sind so zu bestimmen, da die von den Temperaturverzerrungen bewirkten gegen- 


1 E. Tremmel, Osterr. Ing.-Arch. 8 (1954), S. 20 ff. 
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‘seitigen Verschiebungen der freien Endquerschnitte verschwinden. Fiir die auf dic Langen- 
-einheit der z-Richtung bezogene GréBe der Doppelkraft erhalt man den Wert (Abb. 2) 


KL 
N—L 


| = 
= Cof 28 —Cof 2a + 22%) (51) 


/ und fiir den Hebelsarm 
ey N + 1/2 (oj 2 6 —Coj 2«) 


he 2(N,—LN ee) 
Coj 28 —Cof2a + ch Name 
mit den zur Abkiirzung eingefiihrten Konstanten 
N= @ a)? — Gur (p'—a) (53) 
und 
N= (6—a) + + (Gin 28 — Gin 2.) (54) 


‘sowie L und N gemaB (48a) und (48b). 

Die Bestimmung des resultierenden Verschiebungsfeldes bietet weiter keine Schwierig- 
keiten. Die aus den Temperaturverzerrungen allein, also im offenen Querschnitt folgenden 
Verschiebungen sind allgemein durch (3) gegeben. Mit der durch (29) festgelegten Temperatur- 
verteilung erhalten wir fiir die mit cy verkniipften Glieder gema8 (8a) nach Durchfiihrung der 
Integration 


W(t) = Up + ivy =(1 +r) aregt 
und nach Aufspaltung in Real- und Ima- 
ginarteile die Verschiebungen u, v in den 


Richtungen x, y des kartesischen Ko- 
ordinatensystems 


U =(1+»)oreghsinyn, 


1 65 
% =(1lt+r)aregh Giné. 


-wihrend die mit c, verbundenen Ver- 
schiebungsglieder bereits in (35a) und 
(35b) angegeben wurden. 

Durch Zusammenfassung von (55) 
bzw. (35) ergeben sich schlieBlich die Ver- 
schiebungen fiir dasin 7 =a geschlitzte 
Rohr. 

Uberlagert man diesen Verschiebungs- 
' gréBen die nach (10) berechneten, aus der 
Spannungsfunktion y folgenden Aus- 
driicke, dann verbleibt das dem geschlosse- 
‘men Rohr zugeordnete eindeutige Ver- 

schiebungsfeld. Abb. 3. 
| 


6. Zahlenbeispiel. Die Anwendung des hier geschilderten Rechnungsganges soll nur an 
‘einem praktischen Beispiel gezeigt werden. Unter dem Flur der Maschinenhalle eines Kraft- 
hauses sei eine Rohrleitung im Fundamentbeton gemaB Abb. 3 gefiihrt. Die Hallentemperatur 
-werde konstant auf dem Wert T,, gehalten, wahrend die Temperatur an der RohrauBenwand 
lingere Zeitraume hindurch den Wert T; haben mége. Bleibt der EinfluB der AuSentempe- 
‘ratur unberiicksichtigt — bei den verhaltnismaBig groBen Abmessungen derartiger Hallen ist 
diese Vernachlassigung zulassig —, dann lassen sich Temperatur und Spannungsfeld, zu- 
‘mindest in der unmittelbaren Umgebung des Rohres, nach der hier abgeleiteten Berechnungs- 
methode mit hinreichender Genauigkeit bestimmen. 

Wir kénnen die Fundamentoberkante als Kreis mit unendlich groBem Halbmesser und 
naherungsweise konstanter Umfangstemperatur T,, auffassen und der Untersuchung das durch 
diesen ,,Kreis‘‘ und den Rohrquerschnitt gegebene Bipolarkoordinatensystem zugrundelegen. 
‘Mit r, = oo folgt aus (28) Gina = 0 und damit a = 0. Fiir den Rohrumfang gilt (Abb. 3) 


foal = ra +4 


12% 
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wenn der Abstand des Scheitels von der Fundamentoberkante mit d, bezeichnet wird, ferner 


nach (28 
(28) ‘ é 
p= b@tgP, = GaRe 


aus der Division der beiden letzten Gleichungen folgt schlieBlich 
ds 
Co) B=—+1. 
TB 


Damit ist das den hier vorliegenden geometrischen Verhaltnissen angepafte Bipolarkoordi- 
natensystem festgelegt, da mit dem bekannten Wert von f der Mafistabfaktor 6 aus (28) zu 


b=r,z Gin 


berechnet werden kann. Bemerkt sei noch, daB die zweite Koordinate 7 von der Rohrsohle 
aus, wo sie den Wert Null hat, gegen den Scheitel, fiir den 7 = z wird, zu zahlen ist. 

Mit dem hier angenommenen Verhiltnis d,/rg = 1,5 wird Co) 6 = 2,5, und daraus ergibt 
sich weiter § = 1,5668. Fir die zur Bestimmung der Spannungen allein erforderliche Kon- 
stante c, in 


Tee ht ee 
erhalt man nach (30) 
Tp —T» 
C— ene oe 
oder, wenn die Temperaturdifferenz T,; —T, = AT gesetzt wird, 


a = SF = 0.638247. 


Nehmen wir den Elastizitatsmodul des Betons mit E = 2,5 105 kg/em?, die Poissonsche Zahl 


mit y = 1/6 an, dann folgt der Schubmodul 
E 
= = - 105 
C= Fey 1,071 + 105 kg/cm? , 


und wir erhalten mit a7 = 10° fiir die von der Temperaturdifferenz abhangige Konstante K 
nach (43) 


K= —0AT87b AT. 


Zur Berechnung der Spannungen werden ferner noch die oben angegebenen Konstanten L und 


N bzw. J, L, und Ny bendtigt; es ist 
L= + 0,6503 , J = + 0,2396, 
N= —4,1614, 2L,= + 2,8362, 2 NN, = 0.5248. 


Fiihrt man die so bestimmten Konstanten in die fiir die Spannungen berechneten allgemeinen 
Ausdricke ein, dann fallt der MaBstabfaktor b heraus, und wir erhalten z. B. fiir die vor allem ~ 
interessierenden Spannungen o, in den einzelnen ausgezeichneten Punkten der Symmetrie- 
achse die nachstehenden Werte: 


a) an der Fundamentoberkante fiir § = 0,7 =z 
Of STAT 
b) am Rohrscheitel fiir § = 6, 7 =a 


= — 0,705 AT, 
c) an der Rohrsohle fiir = 6,7 = 0 
o, = — 1,498 AT. 
Um Aufschluf8 ibe: die GréBenordnung der gréBten Schubspannung zu bekommen, wurde 


0 
aus (38) 5: = = (0), - = 0 gebildet und die beiden Koordinaten £&, max Yrwas' Gureh: Iteration 


bestimmt. Mit €, mg, = 0,987 und "zr max == 1,135 ergibt sich 7,,,. = -- 0,131 AP, 
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'¥ Die Ergebnisse zeigen, das schon geringe Temperaturunterschiede Spannungen hervorrufen, 
die die Ri®sicherheit des Fundamentbetons betrachtlich verringern kénnen. Mit der, ent- 
_sprechend der Temperatur des Stollenwassers, angenommenen Rohrtemperatur von etwa 8° und 
der Raumtemperatur in der Halle, die 15° betragen moge, folgt eine Temperaturdifferenz 

AT=T,—T,= — 7°; 


° 


damit ergibt sich am Rohrscheitel eine Zugspannung von —0,705 -—7, also rund 5 kg/cm2, 
| wahrend sie ihren Gré®twert mit 10,5 kg/cm? an der Sohle erreicht. Die Druckspannung an der 
| Fundamentoberkante halt sich mit rund 16 kg/cm? ebenso wie die gréBte Schubspannung mit 
rund 1 kg/cm? in maBigen Grenzen. In praktischen Fallen werden neben den Spannungen auch 
| die Verschiebungen zu bestimmen sein, da durch sie gegebenenfalls unerwiinschte Kraftwirkungen 
_ auf das Rohr verursacht werden kénnen. 

Wenn die vorstehende Untersuchung auch wegen der verschiedenen Idealisierungen kein 
_ genaues Bild des Spannungs- und noch weniger des Verschiebungsfeldes liefert, so erhalt man 
doch einigermafen wirklichkeitstreue Aufschliisse tiher die Beanspruchung des Fundament- 
| blockes. 

Demgegentiber mi®Bten hier im Rahmen einer baustatischen Berechnung, bei welcher der 
, zwischen Maschinenhausflur und Rohr liegende Betonkérper als Gewélbe aufzufassen wire, 
) willkiirliche Annahmen iiber die Lage der Kampferquerschnitte getroffen werden; schlieflich 
, aber ware der Arbeitsaufwand bei halbwegs genauer Erfassung des veranderlichen Tragheits- 
moments keineswegs geringer als nach dem hier erérterten scheibentheoretischen Verfahren. 


7. Ausblick und SchluBfolgerung. Ahnlich wie in dem hier behandelten Zahlenbeispiel 
kénnten auch die Spannungen in einem Motorblock mit kreisférmiger, von einem erhitzten 
» Medium durchstrémten Bohrung naherungsweise berechnet werden. SchlieBlich lassen sich 
» auch die Verschiebungsfelder eingespannter Gewélbemauerlamellen in einfacher Weise gemaB 
» (40) bestimmen und mit ihnen annahernd die aus der Einspannung folgenden SchnittgréBen !. 
Die Untersuchung des von exzentrischen Kreisen begrenzten Ringes bringt weiter keine 
' Schwierigkeiten; fiir den ebenen Spannungszustand wird o,, = 0 und ¢,, = ar T, so daf in (3) 
) der Faktor (1 + yw) durch 1 zu ersetzen ist; da ferner der Beiwert des Integralterms in (10) in 


1 1 


OC (lini) 
| iibergeht *, erhalten wir hier fiir den Koeffizienten des Stérungsgliedes in (14) baw. (14a) 
x= —2G6(1+>»)ar. (15a) 


| Ebenso kann auch der Fall einer exzentrisch liegenden Warmequelle in der vollen Kreisscheibe 
) behandelt werden. Diese Berechnungen aber beziehen sich strenggenommen nur auf Scheiben 
mit warmeabgeschirmten Oberflaichen, die fiir die Praxis kaum Bedeutung haben. Das Ver- 
. fahren ware also, wie es fiir zentrische Kreisringe und andere Scheiben bereits durchgefithrt 
' wurde 3, auf die Beriicksichtigung von Warmeverlusten an den Oberflachen zu erweitern. 

Die hier speziell fiir den von zwei exzentrisch liegenden Kreisen begrenzten Querschnitt 
| erhaltenen Ergebnisse gelten in ahnlicher Form, wie schon erwahnt, im Falle stationarer 
| Warmefliisse fiir alle geschlossenen Ringquerschnitte; bei symmetrischen Temperaturver- 
» teilungen miissen sich, da die Vertraglichkeitsbedingungen erfillt und die Rander spannungsfrei 
sind, die Spannungen zu einer im allgemeinen nicht verschwindenden Doppelkraft und einem 
| Momentenpaar zusammenfassen lassen. Bei unsymmetrischen Temperaturverteilungen oder 
( Querschnitten tritt zu diesen SchnittgréBen noch ein Schubkraftepaar * hinzu. Die Span- 
' nungsfunktion y wird daher immer aus den, diesen Belastungsfallen im offenen Ring zuge- 
ordneten Spannungsfunktionen aufgebaut sein. Sie kann also, wenn diese letzteren bekannt 
} sind, worauf noch, ohne auf Einzelheiten einzugehen, hingewiesen sei, unter Verwendung von 
| (14a) unmittelbar berechnet werden. 


(Eingegangen am 28. Mai 1954.) 
Anschrift des Verfassers: Dr. Erwin Tremmel, Zell am See (Osterreich), Haus Lebzelter. 


1 £, Tremmel, a. a. O. S. 34. 

2 C, B. Biezeno u. R. Grammel, a. a. O. S. 112. 
3 E. Melan, Osterr. Bauz. 8 (1953), S. 89 ff. 

4 KE, Tremmel, a. a. O. S. 26. 
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Eine neue hydrodynamische Analogie zur Torsion prismatischer Stabe. | 
Von E. Pestel. 


1. Einleitung. Es ist eine bekannte Tatsache, daB haufig die wesentlichen Eigenschaften | 
von zwei oder mehreren verschiedenen physikalischen Erscheinungen in identischer mathe- | 
matischer Form zum Ausdruck kommen. Wenn eine vollstandige mathematische Identitat | 
besteht, geniigt es, nur eines dieser physikalischen Systeme zu studieren, um gleichzeitig eine 
Aussage auch iiber das Verhalten der mathematisch analogen zu gewinnen. 

Von dieser Méglichkeit konnte man besonders in der Elastizitatstheorie vielfach Gebrauch | 
machen, da hier zahlreiche mathematische Analogien entdeckt wurden. Unter diesen nehmen || 
die Gleichnisse fiir die Aufgabe, den Spannungs- und Formanderungszustand eines auf Torsion | 
beanspruchten prismatischen Stabes zu untersuchen, insbesondere das fiir dieses Problem vor | 
50 Jahren von Prandtl gefundene Membran-Gleichnis, wohl die bedeutsamste Stellung ein, | 
Neben dem Prandtlschen Gleichnis ist fiir eine anschauliche Erérterung des Torsions-Problems 
gelegentlich auch das Thomsonsche Strémungsgleichnis herangezogen worden, das die mathe- || 
matische Analogie zwischen der elastizitatstheoretischen Aufgabe und dem ebenen Strémungs- | 
vorgang einer mit konstanter Wirbelstarke strémenden reibungslosen Flissigkeit benutzt. Im_ 

y Gegensatz dazu liegt dem neuen Strémungsgleichnis || 

das Verhalten einer viskosen Flissigkeit zugrunde. | 

Dabei betrachten wir einen Strémungsvorgang, der | 

wie folgt entsteht. 

Eine Scheibe beliebigen Umrisses sei gelagert 

a auf einer Flissigkeitsschicht, deren Trager eine un- 

| endlich groBe Ebene ist. Scheibe und Ebene sind 

| | parallel zueinander (Abb. 1). Sie werden relativ 

| | gegeneinander lings ihrer Normalen mit endlicher 
| Geschwindigkeit bewegt. Wird dabei die Spalthéhe 


| Scheibe L ; ; 
| / | h zwischen Scheibe und Ebene verringert, so wird 


, PD eo Ebene Flissigkeit aus der Schicht zwischen Scheibe und | 

VTL TTT TTL, Ebene herausgedriickt. Dadurch entsteht in dem 

RE Re me aN Spalt ein i. a. von Ort zu Ort verschiedener Uberdruck 

a ; p gegeniitber dem Druck py auBerhalb der Scheibe. 

Es 1a6t sich nun u. a. zeigen, daB die Kraft, welche die beschriebene Bewegung der Scheibe 

hervorbringt, dem Drillmoment zugeordnet werden kann, das einen prismatischen Stab von 
einem Querschnitt gleich dem Umri®B der Scheibe verdrillt. 


2. Mathematische Formulierung des Torsionsproblems und der hydrodynamischen Analogie. 
Wenn wir unter den bekannten Voraussetzungen der St. Venantschen Theorie fir die Ermitt- 
lung der Schubspannungen und der Formanderungen in einem auf Torsion beanspruchten 
prismatischen Stabe, dessen Achse mit der z-Achse zusammenfallen mige, eine Spannungs- | 
funktion m derart definieren, daB 

a 
ty =— 5 ud ty =o (1) 
mie fiihrt uns die Beachtung der Vertraglichkeitsbedingungen auf die bekannte Differential- 
gleichung 


ep , yp 
gn t 3 = EO. (2) 


Darin ist } der Drillwinkel und G der Schubmodul. Als Randbedingung erhalten wir bei einfach | 
zusammenhangenden Querschnitten 


pr = hy. (3) 
Das den Drillwinkel # hervorrufende Drillmoment ergibt sich dann zu 
M2) | Oar (4) 


wenn I, die Querschnittsflache des Stabes ist. 
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Diesen Gleichungen lassen sich wie folgt die entsprechenden des Strémungsgleichnisses 
gegeniberstellen. Unter der Annahme, daB die Spaltdicke h Klein gegeniiber den Abmessungen 
der Scheibe ist, und ferner auch die zeitliche Anderung von h genugend klein! ist, sind die 
Voraussetzungen fiir die sogenannte schleichende Strémung gegeben, die z. B. auch der mathe- 
matischen Behandlung der hydrodynamischen Gleitlagertheorie zugrunde liegt. Es 1aBt sich 
ferner ohne Schwierigkeit der Nachweis fiihren, daB dann die Geschwindigkeitskomponente 
der Strémung in z Richtung vernachlassigt werden kann. Wenn wir ferner unter Beriick- 
sichtigung der Haftbedingung an Scheibe und Ebene an Stelle der Geschwindigkeitskompo- 
nenten in der x- baw. y-Richtung durch Mittelbildung die an jedem Ort (x, y) im Spalt stré- 
mende Flissigkeitsmenge q, und qy einfiihren, erhalten wir folgende einfache Beziehungen 
zwischen dem Druck p und der strémenden Fliissigkeitsmenge, wenn mit fu die absolute Zahig- 
keit der Flissigkeit bezeichnet wird: 


Op 12 uw dp 12 u y 
ee aa ee und en mead Tek (ft) 


Die Erfillung der Kontinuitatsbedingung fihrt auf folgende Gleichung: 


Oqx Oqy bs 7 
Ox oe dy Be 
Die Elimination von q, und q, liefert dann 
ot Re an Og 9! 
Ox? + dy? be (2') 
Am Rande herrscht der Druck 
Pr= Po: (3’) 


Ferner ergibt sich, wie man leicht sieht, als Kraft P, welche die Scheibe mit der Geschwindig- 
keit h relativ zur Ebene bewegt, 


Pp pak ip, Fy: (4’) 
Fy 
Wir erkennen somit folgende Zuordnung zwischen den wesentlichen GréBen des Strémungs- 
und des Elastizitatsproblems: 
ch damage 
6u + 
—269—> is ne 


: 
/ 


Wenn wir nun zur Behandlung mehrfach zusammenhangender Querschnitte schreiten (Abb. 2), 
haben wir die Gleichungen fiir das Torsionsproblem noch wie folgt zu erweitern. Zu der Rand- 
hbedingung (3) treten die entsprechenden fiir die Rander der inneren Bereiche, so daB wir 
an Stelle von (3) setzen 


2 ky > Po> 
M,—> P, (5) 
6 uw 
Tzy — no Tx >» 
| 6 u Abb. 2. Mehrfach zusammenhiangender 
a he ae Querschnitt. 
. 


| gr, =k; (Ga ere ie (3a) 
ie Serdem hat die Spannungsfunktion y, damit die Eindeutigkeit der Querschnittsverwélbung 
gewahrleistet ist, an den inneren Randern noch folgende Bedingung zu erfillen: 
A da G0.b.) (1, 227 en)s (6) 
R; ov 
‘Darin bezeichnet F; die Querschnittsflache des i-ten inneren Bereiches. Fir das Drillmoment 
erhalten wir 

M,=2fpdF —2hF, +22 kh F;. (4a) 

i=1 


1 Es 1aBt sich zeigen, daB hals genigend klein angesehen werden kann, solange h hv < 0,list; darin 
ist vy die kinematische Zahigkeit der Flissigkeit. 
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In (4a) ist das Integral ither die schraffierte Flache von Abb. 2 zu erstrecken. Entsprechend, 
erhalten wir nun unter Benutzung der Zuordnungen (5) als Randbedingungen fir das Stré- 


mungsproblem 
sity Pr; = Pi (ii Oly Qe eam) (3’a) | 
und damit fiir die Kraft P 
P= f{ pdF —pyFy + > pik. (4’a) 
i=1 
Ebenso erhalten wir als Zuordnung zu (6) 
Op ie jy © (s. 
$3, 4 = — = hk; (i =1,-2, » n) 
Re 
Im Hinblick darauf, daf8 
0 ATH ey 
ee a") 


ist, ergibt sich folgende Beziehung: 
q,ds =hF; (tg Local is (6’) 
R; 


daraus folgt, da® der Druck p;in den inneren Bereichen i so zu halten ist, daB aus jedem Lochz 
soviel Fliissigkeit herausstrémt, wie ein Stempel vom Querschnitt I’; verdrangen wirde, wenn 
er mit der Geschwindigkeit h abwarts bewegt wiirde. Die experimentelle Erfiillung dieser Be- 
dingung gestaltet sich im Gegensatz zu der entsprechenden beim Seifenhautgleichnis sehr ein- 
fach, wenn man iiber den inneren Bereichen Zylindert vom Querschnitt F; entsprechend 
Abb. 3 anordnet. 


Abb. 3. Modell fiir mehrfach zusammenhangenden Querschnitt. Abb. 4. Schmaler Rechteckquerschnitt. 


3. Anwendungsbeispiele fiir das Strémungsgleichnis. Es lat sich an Hand einfacher Bei- 
spiele zeigen, daB dieses Strémungsgleichnis in seiner Anwendung als Gedankenexperiment 
vielseitige Méglichkeiten bietet. 

a) Schmaler Rechteckquerschnitt. Man erkennt hier ohne weiteres, daB 
bei der in der Einleitung beschriebenen Bewegung einer Scheibe von der Umrifform eines 
schmalen Rechteckquerschnittes (Abb. 4) der Druckanstieg, der der Schubspannung gemadB 
(5) proportional ist, in der Mitte der Langseite des Rechtecks am gréften ist. Wir erhalten 
einen Naherungswert fiir den hier herrschenden Druckanstieg durch die Betrachtung eines 
unendlich langen schmalen Rechtecks, bei dem ja nur noch ein eindimensionaler Strémungs- 
vorgang quer zur Langseite stattfindet; d. h. q, wiirde hier identisch verschwinden, und wir 
wirden nach zweimaliger Integration der iibriggebliebenen Gleichung (1’) fiir einen Platten- 
streifen der Lange 1 (J > d) angenahert folgende Druckkraft erhalten, wenn p, = 0 gesetzt 
wird: 

(eee (7) 
h 
Mit Benutzung der Zuordnungen (5) P—> M, und Ey) > GO ergibt sich der Torsionswider- 
stand des Querschnittes sofort zu 
M:z 1d 
J= GB =s a . (8) 
Dieser Wert ist je nach dem Verhiltnis //d etwas zu hoch, wie man im Hinblick auf das Stré- 
mungsgleichnis leicht einsieht, wenn man bedenkt, daB bei endlichem Wert von 1 qy nicht ver- 
schwindet, und damit der Druck p mit wachsendem |y| abnimmt. 


; 1 Der Verfasser verdankt einer Diskussionsbemerkung von Herrn Dr. Berger, Wien, den Hinweis auf 
diese experimentelle Méglichkeit. 
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Die Schubspannungen lassen sich am Ort y = 0, wie folgt ermitteln. Hier ist aus Symmetrie- 
grinden q, = 0; also liefert die Kontinuitatsbedingung 


Petoihin 
Wenn wir wieder von der Voraussetzung Gebrauch machen, daB 1S d, erhalten wir wegen 
h = — Ph'/uld° [vgl. (7)] und im Hinblick auf (5): 
he 
Ge Uae und P—>M,, 
M: 
Tzy => 6 id Xs 

insbesondere fiir x = d/2 

3 

te sy mea << l@ M,. (9) 


Der so erhaltene Naherungswert ist fiir endliche Werte I/d (J > d) etwas zu klein, da bei ge- 
gebenem M,, also vorgeschriebenem P, der Druckanstieg dp/dx an der Stelle (x = + d/2; 
y = 0) gréSer sein muB, wenn q, fiir y + 0 nicht verschwindet und damit dp/dx mit wachsen- 
dem |y| abnimmt — fiir (x = + d/2; y = + 1/2) wird ja dp/dx = 0 —, als wenn, wie hier vor- 
ausgesetzt, gyidentisch verschwindet und somit dp/dx 
bei festgehaltenem x konstant ist. 

Man erkennt ferner, da8B es fiir die gréBte Schub- 
spannung und den Torsionswiderstand unwesentlich 
ist, ob die Langseite des schmalen Rechteckquer- 
schnittes gerade oder gekriimmt ist (Abb. 5). 


Ss 


Abb. 5. Gekriimmter schmaler Rechteckquerschnitt. 
b) Aus schmalen Rechtecken zusam- 
mengesetzte Querschnitte. Wir betrachten 


‘ rm 
gewalzte, geschweiBte oder genietete einfach zu- 
sammenhangende Querschnitte entsprechend Abb. 6, 5 Fee 
deren wesentliches Kennzeichen darin zu suchen 


ist, daB der Stabquerschnitt aus Querschnittsbestand- Abb. 6. Zusammengesetzte Querschnitte- 
teilen ahnlich den unter a) behandelten aufgebaut ist. 
Es leuchtet hier sofort ein, daB die Kraft P, die fiir die in der Einleitung beschriebene Be- 


wegung der Gesamtquerschnittsflache F mit der Geschwindigkeit h aufgebracht werden muf, 
nur wenig gréBer ist als die Summe der Krafte P;, die jeweils fiir die Bewegung der Quer- 


schnittsbestandteile F; mit der gleichen Geschwindigkeit h erforderlich sind. Folglich gilt, 
wenn F — 5) Fist; 
j 


Pe Py 
df 


Da unter der Voraussetzung (I; > d;) P; nur wenig kleiner als die rechte Seite von (7) ist, kann 
naherungsweise 


h 
Px oF > 1d} (10) 
j 
gesetzt werden. Im Hinblick auf die Zuordnungen (5) folgt dann 


1 
See Syd, (11) 
J 


wenn die Lange jedes Querschnitthestandteiles I, und seine Dicke d; betragt. 
Im Hinblick auf (10) ibernimmt wegen 


Pine i 
Pp? M, 
jeder Querschnittsbestandteil j folgenden Bruchteil des gesamten Torsionsmomentes 17 5 
Ijdj 
Je ae ys 8 
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Wegen (9) gilt 


3 
gk Sere q; = 
also 
See dj 
Tj max © qt — 4, M, =3 —+, z 
jG jd; hile] 
j 


Die gréBte Schubspannung ist also in dem dicksten Querschnittsbestandteil zu erwarten; sie 
betragt naherungsweise 


dma r 

Tmax aaa * VM, (12) 
21; dj 
j 


c) Ableitung der Formel von Bredt. Bei prismatischen Hohlstaben geringer Wand- 
starke gelingt die Bestimmung der interessierenden GréBen mit Hilfe des Strémungsgleichnisses 
in besonders anschaulicher Weise. Nach den Ausfiihrungen am SchluB des Abschnittes 2 (vgl. 
auch Abb. 3) wiirde das einem diinnwandigen Hohlquerschnitt entsprechende Analogiemodell 
etwa, wie in Abb. 7 gezeigt, auszubilden sein. Die Breite 6(s) des ,,Hutrandes des unten 
offenen und oben geschlossenen Zylinders ist jeweils gleich der Wandstairke des Hohlquer- 
schnittes. Aus der Abb. 8 ist zu ersehen, wie man sich die Durchfiihrung eines solchen Modell- 


Ss 
p 
4 
Vv 
Abb. 7. Modell fir diinnwandigen Hohlquerschnitt. Abb. 8. Modellversuch und Druckverteilung. 


versuches vorzustellen hat, und welche Druckverhaltnisse sich dabei ergeben. Unter der Vor- 
aussetzung, daB der ,,Hutrand“ iiberall sehr schmal ist, kann der Druckabfall von p; auf py 
hinreichend genau als geradlinig aufgefaBt werden. Daraus folgt, daB P = (p; — py) Fin ist, 
worin F’,, die von der Mittellinie des Hutrandes umschlossene Flache ist. Wegen 


Op). Pi Pe ae 
dv 0(s) ns dy 
[vgl. (1’’)]ergibt sich 
12 
P = —* q, 8(s) Fn- (13) 
Im Hinblick auf die Zuordnung (5) : 
6 UL. 
“+ q —> Ts 
und P > M, folgt 
VM, ee OSes (14) 


Aus p; — Po = konst. folgt ferner r, 5(s) = konst.; d.h. der sogenannte ,,Schubflug* ist 
konstant. Da nach (6’) und (13) ; 
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ist, liefert die Zuordnung (5): ot —+>— 2G und P- M, die von Bredt auf anderem Wege 


gefundene Formel 


M: ds 
~ 46F% a as) So 
und damit den Torsionswiderstand zu 


a 


F rs : (16) 
Rn os) 


d) Hohlquerschnitte mit Zwischenstegen. Zunachst betrachten wir den einfachen 
Fall eines symmetrischen, diinnwandigen Kastenquerschnittes mit Zwischensteg (vgl. Abb. 9). 
Offenbar strémt dann bei dem hydrodynamischen Modellversuch keine Fliissigkeit von Raum I 
nach Raum II, so daB p,; = p,,; ist; d.h. der ,,Schubflu®“ im Sinne des Abschnittes c) ist Null. 
Daraus folgt aber nicht, daB im Querschnitt des Zwischensteges tberhaupt keine Schubspan- 
nungen wirken, wie in der einschligigen Literatur zuweilen behauptet wird. Die Schub- 


a 


Abb. 9. Symmetrischer Kastenquerschnitt (a) Abb. 10. Unsymmetrischer Hohlquerschnitt (a) 
und Modell (b). und Modell (b). 


_ spannung ist dort lediglich sehr klein. Wir erkennen mit Hilfe unseres Strémungsgleichnisses 
das Vorhandensein von Schubspannungen sofort, da bei Abwartsbewegung des hydrodyna- 
mischen Modells der Zwischensteg Flissigkeit symmetrisch in die Raume I und II verdrangt. 
: Die gréBte im Zwischensteg auftretende Schubspannung wiirde dann durch (9) gegeben sein, 
- worin jedoch fir M, nur der Bruchteil des gesamten Torsionsmomentes einzusetzen ist, der 
- von dem Steg ibernommen wird. Dieser ist jedoch so klein, da wir den Anteil des Zwischen- 
steges bei der Ermittlung des Torsionswiderstandes vernachlassigen kénnen, der somit durch 
(16) hinreichend genau angegeben ist. 
| Beim unsymmetrischen Fall (z. B. Abb. 10) liegen die Verhaltnisse verwickelter. Wir 
wollen uns auch hier auf die Betrachtung des diinnwandigen Falles beschrinken, und daher 
_ lediglich die Schubspannungen, die am ,,SchubfluB« beteiligt sind, beriicksichtigen. Die Kraft 
| P, die zur Abwartsbewegung des Modelles mit der Geschwindigkeit h aufgebracht werden muB, 
| ist dann praktisch hinreichend genau durch 


P = (p; — Po) Fy + (Py-— Po) Fitmn 

_ gegeben. Da j 

) Py — Po = 11 und pir — Po = on Om (17) 
ist, ergibt sich 
P= ens (qv Or Fy, + Gy 11 Orr F'tr,,) 3 (18) 


- ferner folgt 


12 
(Py — Po) — (Prx — Po) = Px — Prx = “po (4v  t — Go 12 O11) - (19) 
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Aus Kontinuitatsgriinden ist 


Pyke = [ards + [armas = arr [5 +oinérn | p> | (20) 
: ds ds 
Fy h a Jarnas a Jorn ds =andn | 52 UES SAS in | 5 3 | 
fiir q, 111 erhalten wir wegen (1) die Beziehung 
eH ain Oru = Pr — Pu: (21) 
Im Hinblick auf (19) und (21) ist also 
Qiu Orn = 191 — Wu On: (22) 


Die Gleichungen (17) und (21) zeigen, daB wegen der Zuordnung (5): Ay >t, der. Schub- 
fluB in den einzelnen Bestandteilen des Querschnittes jeweils konstant ist, also 
7 Ore ONsi a, 1; 0y; = Konst. und. 1, 0,7, = Konst. 


Ferner gilt wegen (22) 

Tr Oras = Tr Oy — Ty On - (23) 
Als Drillmoment liefert (18) 

M, = 2 (110191, +1 Ou Fn,); (24) 
und schlieBlich erhalten wir aus (20) im Hinblick auf die Zuordnung (5) : 


Sti, > —268 und ae 


il ds ds 
Dies 0; (% : la peau iu fae)? | 


1 ds ds 
ee cee, leo 


Die Gleichungen (23) und (25) reichen aus, bei gegebenem Drillmoment den SchubfluB im 
Querschnitt und den Drillwinkel # zu bestimmen. 


(25) 


und 


4, Zusammenfassung. Die Ausfiihrungen im Abschnitt 3 haben gezeigt, daB das neue 
Strémungsgleichnis auf anschauliche Weise quantitative Aussagen tiber das Torsionsverhalten 
prismatischer Stabe liefert, die zwar auch durch die Anwendung der Membrananalogie und des 
Thomsonschen Strémungsgleichnisses, jedoch in ihrer Gesamtheit nicht durch nur eines der 
letzteren allein durch einfache anschauliche Analogiebetrachtungen gefunden werden kénnen. 
So bietet das Seifenhautgleichnis z. B. keine Méglichkeit, die unter 3c) und 3d) erérterten 
Probleme anschaulich zu behandeln. 

Im itbrigen besteht berechtigte Aussicht, mit Hilfe der vorgelegten hydrodynamischen 
Analogie das Problem dickwandiger Hohlquerschnitte von beliebiger Querschnittsform ex- 
perimentell verhaltnismafSig einfach zu lésen, eine Aufgabe, deren Bewaltigung bei Benutzung 
der bisher bekannten Analogien erhebliche technische Schwierigkeiten bereitet. Der Verfasser 
hofft, in Zusammenarbeit mit dem Festigkeitslaboratorium der Technischen Hochschule 
Hannover (Leiter: Prof. Dr.-Ing. 0. Flachsbart) noch in diesem Jahre eine geeignete Versuchs- 
apparatur zu erstellen. . 


(Eingegangen am 26. Juni 1954.) 
Anschrift des Verfassers: Prof. Dr.-Ing. E. Pestel, Hannover, Nienburger Str. 2. 
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Uber Tragfliigelsysteme in ebener Strémung 
bei beliebigen instationaéren Bewegungen. 


Von K. Nickel. 


1. Einleitung. Das Problem des instationar bewegten Tragfliigels ist schon sehr haufig 
behandelt worden. Der einfachste Fall der ebenen, inkompressiblen Strémung um den unend- 
lich langen Tragfliigel wurde zuerst von W. Birnbaum‘ und H. Wagner ? aufgegriffen. Unter 

| gewissen Vernachlassigungen (diinnes, wenig gewolbtes Profil, usw.) fiihrt diese Aufgabe auf 
die Loésung von zwei Integralgleichungen. Die eine ist eine Fredholmsche Integralgleichung 
erster Art und als Integralgleichung der ebenen Wirbelschicht aus dem stationdren Falle wohl- 
) bekannt. Die zweite, eine Volterrasche Integralgleichung erster Art vom Sonineschen Typ, 
beschreibt den instationéren Anteil der Strémung und ist zuerst von Wagner? aufgestellt 
worden. Wahrend die Integralgleichung der Wirbelschicht schon 1919 von A. Betz * aufgelist 
werden konnte, gelang die Lésung der Wagnerschen Integralgleichung in geschlossener Form 
durch Hankelsche Funktionen erst 1936 H. G. Kiissner 4. 
Damit lieBen sich geschlossene Formeln fiir Auftriebsverteilung, Auftrieb und Moment 
_gewinnen in vollkommener Analogie zum stationdren Fall (man vgl. H.G. Kiissner® und 
H. Séhngen °). Als Nebenergebnis stellte sich dabei noch heraus, da beim Eintreten des Fliigels 
| 


in eine beliebige, raumfeste Bée die Auftriebsverteilung tiber die Fliigeltiefe proportional zur 
Auftriebsverteilung der stationar bewegten ebenen Platte wird. Es existiert also zu jedem Zeit- 
punkt ein Anstellwinkel, fiir den die stationaér angestrémte Platte dieselbe Auftriebsverteilung 
_aufweist, wie der Fliigel in der Boe. Damit greift weiterhin der gesamte ,,Boéenauftrieb“ im 
vorderen Neutralpunkt des Fliigels an, was auch durch Versuche ” bestatigt wurde. Dieses 
Resultat gilt ibrigens nicht nur fiir Béen, sondern ganz allgemein fiir beliebige ,,strémungs- 
feste‘“* Bewegungen, so z. B. (vgl. H. G.Kiissner > S. 356) fiir die aus der Ruhe heraus mit kon- 
stanter Geschwindigkeit und konstantem Anstellwinkel bewegte Platte. Das Wort ,,strémungs- 
fest“‘ bezieht sich dabei nicht auf die Bahn der Fligelpunkte, sondern auf die értlichen Normal- 
geschwindigkeiten am Profil. 

Dem Ubergang von der ebenen, inkompressiblen Strémung zum dreidimensionalen Fall 
des endlich langen Tragfliigels in kompressibler, instationarer Strémung wurden in letzter Zeit 
mehrere Arbeiten gewidmet, man beachte insbesondere H. G. Kiissner ®. 

Die folgenden Untersuchungen erweitern die oben geschilderte Theorie der instationaren 
Bewegungen eines Einzelfliigels auf Systeme aus beliebig vielen und beliebig angeordneten 
Fliigeln®, jedoch unter den beiden Einschrankungen der ebenen und inkompressiblen Strémung. 
-AuBerdem werde noch, wie in allen zitierten Bearbeitungen, von der Zahigkeit abgesehen. Die 
Fliigelprofile seien diinn, schwach gewolbt und ihre Sehnen annahernd parallel und wenig gegen 
jdie Strémung angestellt; weiter verandere sich die Strémungsgeschwindigkeit wenig. Nach 
der Birnbaumschen Theorie der diinnen Profile! wird dieses Fliigelsystem dann in Ziff.2 er- 
'setzt durch eine geeignete Verteilung von Wirbeln der Dichte I’ auf den Profilsehnen und von 
jabschwimmenden freien Wirbeln der Dichte ¢ in ihrer Verlangerung (vgl. Abb. 3). Fihrt man 
‘in Analogie zur Theorie der stationdren Bewegung noch gebundene Wirbel der Dichte y pro- 
portional zu den Auftriebsverteilungen iiber die Profiltiefen ein und ersetzt ° durch y, so findet 
/man (Ziff. 3) entsprechend zur Theorie des Einzelfliigels aus der physikalischen die folgende 
t mathematische Problemstellung: die y geniigen einem System von F'redholmschen Integral- 


1 W. Birnbaum, Z. angew. Math, Mech. 4 (1924), S. 277. 
H, Wagner, Z. angew. Math. Mech. 5 (1925), S. 17. 
Geschichtliche Angaben dazufindet man bei L. Schwarz, Luftfahrtforschung 17 (1940), S. 381. 
H. G. Kiissner, Luftfahrtforschung 13 (1936), 5. 410. 
H. G. Kiissner, Luftfahrtforschung 17 (1940), S. 355. 
H. Séhngen, Luftfabrtforschung 17 (1940), S. 410. ; 
H.G. Kiissner, Luftfahrtforschung 13 (1936), S.425. Weiter beachte man die Arbeit von 
\H. G. Kiissner, Z. Flugwissenschaften 3 (1955) S. 1. 
| 8 H. G. Kiissner, Journal of the Aeronautical Sciences 21 (1954), S. 17. 
® Zwei hintereinanderliegende Fliigel betrachtet J. R. M. Radok, The Aeronautical Quarterly, III 
(1952), S$. 297, jedoch ohne Beriicksichtigung der Riickwirkung des Hinterfliigels auf den vorderen. 


} 
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gleichungen erster Art, dem Gleichungssystem der ebenen Wirbelschichten, und die ¢ werden 
durch ein System von Sonineschen Integralgleichungen bestimmt. Schon zur Aufstellung des 
Sonineschen Integralgleichungssystems ist allerdings die Kenntnis spezieller Lésungen des 
Fredholmschen Integralgleichungssystems erforderlich. Ohne diese Integralgleichungssysteme 
aufzulésen, kann man aus dem ersten sofort einen allgemeinen Sachverhalt ablesen, der den 
oben beschriebenen Béensatz auf Fliigelsysteme verallgemeinert und der schon hier ausge- 
sprochen werden soll: 

Satz: Fir alle strémungsfesten Bewegungen eines beliebigen Fliigelsystems der beschrie- 
benen Art gibt es zu jedem Zeitpunkt einen ,,Satz‘* von Anstellwinkeln, fiir den das ent- 
sprechende System aus ebenen Platten bei stationirer Bewegung dieselben Auftriebsvertei- 
lungen besitzt. Dabei stimmen fiir jeden Zeitpunkt die Anstellwinkel der einzelnen Platten 
i, a. nicht iberein. 


Die exakte Lisung des allgemeinen Integralgleichungssystems der ebenen Wirbelschichten 
steht noch aus, doch ist sie wenigstens fiir Sonderfalle bekannt (vgl. K. Nickel1*). Zwei 
dieser Sonderfalle werden in den Ziffern 4 und 5 herausgegriffen, das System aus endlich vielen 

Tragfligeln hintereinander und das senk- 
ormn__ 6 6-0 0-0-0 0)0 00 00 ) Oro - =—rechte Fligelgitter. In beiden Fallen 
lassen sich geschlossene Ausdricke fir 
Auftriebsverteilung, Gesamtauftrieb und 
Gesamtmoment angeben, ohne daB das zugehérige Wagnersche Integralgleichungssystem ge- 
lést zu werden braucht. Fir das Fligelgitter sind diese Formeln besonders einfach gebaut 
und entsprechen vollstandig denen fiir den Einzelfliigel, in die sie fiir unendlich werdende 
Gitterteilung auch iibergehen. Trotzdem das Integralgleichungssystem der Wirbelschichten 
im allgemeinsten Fall noch nicht gelést ist, 1aBt sich dafiir doch wenigstens der Weg zur Auf- 
stellung des Wagnerschen Integralgleichungssystems und seiner Lisung skizzieren (Ziff. 4). 

SchlieBlich werde der oben ausgesprochene allgemeine Satz noch etwas naher betrachtet. 
Fir ein einzelnes Fligelprofil ist er offensichtlich mit dem Kiissnerschen Béensatz identisch. 
Bei mehreren Fligeln iiberrascht zunachst, 
da8 darin nur von einem stationaér bewegten 
Ersatzsystem aus ebenen Platten, nicht aber 
wafborncc vom Angriffspunkt der Luftkrafte gesprochen 

wird. Da namlich nach Kiissner der Einzelfligel 
: in einer Bée durch einen einzelnen tragenden 
Wirbel im Neutralpunkt ersetzt werden kann, 
erwartet man zunachst von einem Fligelsystem 
dasselbe, da®B also auch bei ihm fir beliebige 
; ‘ 4 ae strémungsfeste Bewegungen der Zusatzauftrieb 
Ae Sa ale eal peer eee oe ee en 
die Anfahr- und se mae a Ss einen einzigen Wirbel nicht mehr zutrifft, kann ohne Rechnung an 
einem Gegenbeispiel nachgewiesen werden: 

Es seien zwei gleiche Fliigel unendlicher Spannweite in Tandemanordnung gegeben, deren 
waagrechter Abstand gegeniiber den Fligeltiefen sehr gro sein soll. Damit ist die Riick- 
wirkung der beiden Fligel aufeinander auch sehr klein geworden, wenigstens in stationarer 
Strémung. Betrachtet man jedoch z. B. den Anfahrvorgang, so gerat der Hinterfliigel nach 
einiger Zeit in die Anfahrwirbel des Vorderfliigels (vgl. Abb. 1). Wahrend bis dahin beide 
Fligel annahernd dieselben Anstrémverhaltnisse und damit in jedem Augenblick auch nihe- 
rungsweise denselben Auftrieb besaBen, kommt der Hinterfligel nun zunachst in ein Aufwind- 
und danach Abwindfeld. Da die abschwimmenden Anfahrwirbel in sich die umgekehrte Zir- 
kulation des Vorderfliigels vereinigen, ist ihre Wirkung auf den Hinterfliigel nicht mehr ver- 
nachlassigbar klein. Wahrend also der Gesamtauftrieb im ersten Augenblick des Anfahrens — 
genau so wie nach sehr langer Zeit — in der Mitte zwischen den beiden Fliigeln, also in ihrem 
Neutralpunkt angreift, kann dies fiir den in Abb. 1 skizzierten Zeitpunkt nicht der Fall sein. 
Fir allgemeinere Fliigelanordnungen gilt offensichtlich dasselbe. 

Es gibt manche praktisch vorkommenden Fille, fiir die die nachstehend entwickelte 
Theorie von Interesse ist, so z. B. fiir instationére Bewegungen der Fliigelsysteme Tragfligel — 


Abb. 1. Anfahrvorgang bei einem Doppelfliigel. 


CoC Kandwirbe/ 
== Anfahrwirbe/ 


1 K, Nickel, Ing.-Arch, 20 (1952), S. 363. 
2 K. Nickel, Ing.-Arch, 22 (1954), S. 108 und Ing.-Arch, 23 (1955), S. 102. 


XXIII. Band 1955 Nickel: Uber Tragfliigelsysteme in ebener Strémung usw. 181 


Hohenleitwerk, Vorfliigel—Hauptfliigel—Querruder oder z. B. fiir die Gittersysteme aus 
Leit- und Laufradern von Turbomaschinen. Nur ist bei der Anwendung auf diese Falle der 
Einflu8 der endlichen Spannweite zu beachten, durch den hinter den Fligeln ein System aus 
Anfahr- und Randwirbeln zuriickbleibt (vgl. Abb. 2 fiir die Kombination Hauptfligel — 
HéhenJeitwerk). Diese Raumlichkeit, die hier nicht beriicksichtigt wird, mu selbstverstand- 
lich noch eigens erfaBt werden. 

Die folgenden Ausfithrungen schlieBen sich in Methode und Bezeichnungsweise eng an die 
Arbeit von Séhngen ? fiir den Einzelfliigel an. Wahrend jedoch bei Séhngen den auftretenden 
mathematischen Fragen ein sehr breiter Raum gewidmet wird, soll hier der Kiirze halber 
darauf nicht eingegangen werden. Beziiglich der Zulassigkeit der gebrachten Umformungen 
werde also angenommen, daf die betrachteten Funktionen geniigende Stetigkeits- und Diffe- 
renzierbarkeitseigenschaften besitzen mégen. Ubrigens bietet hierfiir der Ubergang von 
Séhngens Betrachtung eines Fliigels auf Fliigelsysteme keine besonderen Probleme, weil die 


_ Vertauschung von Integration und Differentiation, partielle Integration, Vertauschung der 


Integrationsreihenfolge usw., nur bei der Eigeninduktion eines Fligels auf sich selbst Schwie- 
rigkeiten bereitet, nicht aber bei der Wirkung eines Fliigels auf einen anderen, und gerade diese 


_ Fragen sind ja von Séhngen behandelt worden. Beziiglich der eindeutigen Lisbarkeit der aufzu- 
 stellenden Integralgleichungssysteme jedoch lieBen sich, wie schon oben angedeutet wurde, 


Séhngens Existenz- und Einzigkeitshbeweise bisher leider nicht allgemein erweitern. 


2. Aufstellung der Grundgleichungen. Die betrachteten diimnen und wenig gewélbten 


_ Tragfliigelprofile seien in der komplexen z-Ebene annahernd parallel zur reellen Achse ange- 


i 


' Profilsehnen liegen und in den Ablise- 


t 
: 


ordnet, und werden von links her mit der annaihernd konstanten Geschwindigkeit u, ange- 
strémt (mitbewegtes Koordinatensystem). Die instationire Bewegung mige zur Zeit t = 0 
plétzlich beginnen, zur Zeit t > 0 habe sich die z-Ebene um die Strecke s(t) seit t = 0 ver- 
schoben *. Fliigeldeformationen seien zugelassen, doch seien die Geschwindigkeiten der ein- 
zelnen Fligelpunkte in Richtung der imaginaren Achse klein gegen u,, , ferner mégen sich die 
Fligeltiefen nur unwesentlich andern. Unter den Vernachlassigungen der Theorie diinner 
Profile * wird dann das y-te Profil ersetzt durch eine Wirbelbelegung [',(z,,s) auf der Profil- 
_,sehne‘‘. Diese ,,Sehnen“ seien alle parallel zur reellen Achse und mégen die Profile nach Lange 
und gegenseitiger Lage ,,méglichst gut**+ 
ersetzen (vgl. Abb. 3). Weiter werde an- 
genommen, da die sich von den Profil- 
hinterkanten z, = 6, ablésenden freien 
Wirbel ¢,(z,,s) in der Verlangerung dieser 


punkten haften bleiben. Diese Wirbel- 
schleppen sind also alle ebenfalls parallel 
zur reellen Achse und haben alle die gleiche Linge s. Die eingefiithrten Wirbel induzieren im 
Punkte z, eine Storungsgeschwindigkeit mit den Komponenten u,(z,, s) und v,(z,,s). Langs 


Abb. 3. Ersatzsystem fiir den y-ten Fliigel. 


der w-ten Sehne bestimmen sie sich aus der Gleichung 


b,+s 


by 
. l ES Vo L v Cy, 
ile 9) —¥ (00) = 32 Ge tan eee: (1) 
v a, v b, 


- Hier und im folgenden seien die Summen iber alle Fliigel zu erstrecken; unter dem Symbol f 


moge der Cauchysche Hauptwert des Integrals verstanden sein. 
Gegeniiber der gro®en Anstrémgeschwindigkeit u, werden, wie iblich, die Stérkompo- 
-nenten wu, vernachlassigt. Soll nun das Profilskelett nicht durchstrémt werden, so mu 


U (24. 8) = — W,(Z,, 8) (2) 


_gelten, wenn mit w,(z,, s) diejenige Geschwindigkeit bezeichnet sei, mit der die Strémung auf 


: 
| 


7 


1 Siehe FuBnote 6 von Seite 179. 

2 Wegen der Linearitat der aufzustellenden Gleichungen kann dazu noch ein stationares Glied additiv 
iiberlagert werden, z. B. eine konstante Profilwolbung. Dieses Glied berechnet sich nach bekannten 
Gleichungen und kann daher hier weggelassen werden. 

3 Siehe FufBnote 1 von Seite 179. : 

4 Dieser Ersatz ist in gewissem Umfange willkiirlich und ]a$t sich im Rahmen der Birnbaumschen 
Theorie nicht genauer bestimmen, 
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dem Profilpunkt bei z, zur Zeit t(s) auftreffen wiirde, wenn keine Wirbel I’, und ¢, vorhanden 
wiren. Diese Anstrémgeschwindigkeit w,, setzt sich zusammen aus Beitragen von den Eigen- 
bewegungen der Fliigel, den Deformationen der Profile und den duSeren Stérstrémungen 
(z. B. Béen). Nach (1) und (2) geniigen also die Funktionen J’,(z,, s) und ¢,(z,, s) den Integral- 


bedingungen 
b, +s 


10, (29 8) = so dere, s) Rel =} &, — 2 J e(Cy 8) Ret =z} &, (3) 


fiir..¢, S\4, 0b, und alle tie 1 Toei cs 
Die Forderung nach Erhaltung der Gesamtzirkulation um jeden einzelnen Fligel samt 
zugehériger Wirbelschleppe bedeutet 


by +s 


bn 
if Ty Coes) - ei(Gars)dg pea 0 (je = 2.ee). (4) 


Soll die Strémung an den Profilhinterkanten glatt abflieBen, so miissen dort die Funktionen 
I’ (2,, 8) und ¢,(z,, 8) stetig ineinander tibergehen, es muB also gelten 
LP (Bus 8) = &u(Bys $) (FP See) (5) 
Das vorausgesetzte Haften der freien Wirbel im Ablésepunkt driickt sich aus durch die For- 
derung 
EnlZy2 8) = E,(By, S—2Z, + O,) (Tipe PR Ee) 5 (6) 
Es wird im folgenden unterstellt, da das durch die Gleichungen (3), (4), (5) und (6) ge- 
kennzeichnete Problem eindeutig lésbar ist, d. h. daB bei bekannten Normalgeschwindigkeiten 
w,(u =1,2,...) stets genau ein ,,Satz* von Wirbelverteilungen J’, und ¢, (u = 1, 2,..-) 
existiert, der die Gleichungen (3) bis (6) befriedigt. In Sonderfallen (Ziff. 4 und Ziff. 5) konnte 
dieser Nachweis im Anschlu8B an Séhngen erbracht werden, doch soll hier darauf verzichtet 
werden. Der allgemeine Fall dagegen steht noch offen. 
Zur Bestimmung der Auftriebsverteilung tiber die Fliigeltiefe werde die Dichte y,(z,, s) der 
gebundenen Wirbel durch die Gleichungen 


Pues S)u ah Pile pe S)o Css Y lB yo s) (u = 1, 2, ae -) (7) 
eingefiihrt, worin p,(z,,s) den Druck au der Stelle z, des w-ten Profils je nach dem Index u 


oder o auf der Unter- oder Oberseite bedeutet. @ sei die Dichte des strémenden Mediums. Die 
Funktionen y,, und J", hangen dann ' durch die Relation 


Vlg 8) = Pyle 8) + i PCs) dee tee ee ee (8) 


zusammen, 


3. Umformungen der Grundgleichungen. Das Ziel der nachfolgenden Umformungen ist die 
Gewinnung von Gleichungen, in denen die J", gema8 (8) durch die y,, ersetzt sind. Gleichung (3) 
integriert und durch partielle Integration umgeformt gibt, wenn man (4) beachtet, 


z 


u b, cy 
1 1 
| wor s) de, = Fe pa | L(y» s) dn, Re 2 BEY A, | dé, a 


a 


u 
b, +s by 


iat | | e,( 8) dn, Re | ab Z| dt, — 2 j TAG s) Re {log (a, —€,)} al, — 


au 
» 
b, b,+s 
by +s 


< HD | ers) Re flog (0, —£)} dy. 


* Man vel. etwa H. Séhngen (FuBnote 6 von Seite 179), dochist zu beachten, daB Séhngen entgegen- 
gesetzt zur hier gebrachten Darstellung eine Anstrémung aus Richtung der positivenreellen Achse 
annimmt, 
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(Um eine eindeutige Funktion Logarithmus zu erhalten, denkt man sich fiir jeden Fliigel die 

z-Ebene von b, nach oo so aufgeschnitten, daB dieser Schnitt keine andere Profilsehne schneidet. 
Fir die praktisch interessierenden Fligelsysteme ist das immer miéglich. Damit ist zwar der 
 Logarithmus fiir jeden Fligel anders erklart, doch beeintrachtigt das die weiteren Uber- 
legungen nicht.) 


Durch Differentation nach s findet man daraus mit (6) 


*u by fy 
0 1 o 1 
Be Oe (ie s) dt, = > ie | LAN s) dn, Re peal dl, aI 
ay, ; is a, a, 
by +s 


Face | lbs) Rel ela, — auf) (9) 
aM 
mit 
b, 
BS) = 25> ee | TH w 9) Re flog (a, —&)} aly + 


b,+s 


ERD | elles Re flog — 09} 
b, 


Addiert man (3) und (9), so lassen sich schlieBlich nach (8) daraus die Funktionen I’, eliminieren 
und man kommt auf das Fredholmsche Integralgleichungssystem erster Art 


b, 
=> Grits s) Re ie 5| Gee a, 5) es) = © eld (10) 
7 mit 


*u 
0 
H (2,4 8) = W (Zs s) + me | OS (Ge, s) dl, 
ay 
Dies ist das wohlbekannte Integralgleichungssystem der diinnen Profile bei stationarer 
Strémung, wenn man fiir den Augenblick s als einen Parameter ansieht. Die Lésung von (10) 


ist wie dort zu bestimmen unter der Forderung glatten Abflusses an den Profilhinterkanten, 
nach (5), (4) und (8) muB also gelten 

| : VAD nas) 0 Cig 1 2s a) (11) 
Fir den speziellen Fall der strémungsfesten Bewegung gilt 

a) 7) 

| Cap WylZy> 8) + Z Wyleu J=o 

und also auch 

| H (2,4, 8) = h,{(s) - (12) 
‘Damit geniigt also die Dichte y,,(z,, s) der gebundenen Wirbel bei strémungsfester Bewegung 
dem Integralgleichungssystem 


(fyls) = By(s) + h,(s)) 
und dies ist fiir jedes s gerade das Integralgleichungssystem fiir ein System aus ebenen Platten 
mit den Anstellwinkeln 


,,(s) = = fds) i eee soe (13) 


Damit ist der in der Einleitung (Ziff. 1) formulierte Satz vollstandig bewiesen. 
13 


184 Nickel: Uber Tragfliigelsysteme in ebener Strémung usw. Ingenieur-Archiv 


4, n Profile hintereinander. In diesem Sonderfall sind von selbst alle Strémungskomponenten 
u, = 0, man kann also in Gleichung (3) den Operator Re {---} weglassen. Weiter kann noch 
attne Kinschrankung z,, durch x, ersetzt werden. Damit geht (10) tiber in 


b, 
ui < wr, 
eat ree Q dé, = H,(%ys s) shy &,AS) (uu ee 1 #18 25 n) . (10b) 


Die Auflésung y,(€,,s) (v =1,...,) von (10b) unter der Joukowskibedingung (11) lautet 
(vgl. K. Nickel*) 


tal (ei ca Uae et Ey 3¢ Ay Xp, Jia £u(S) IT = Xu (v=1, ..., n)y 
Au ax Be Xu fy Gate Xu | 
(1 


Der gesamte Auftrieb aller n Fligel yeertee sich zu | 


b, by, 
sone > | mlb 9) d= 20m >) | i w 5) + g,(s) wo) — a [Ta>* dx,, (15) 
car ay B= on | 


und das Moment um den Nullpunkt lautet 


Ms) =o. > | Eu 8) dB, = 


= 294, Sf [H,(%, 8) +8,(s otf, - FS eg] ]/— Te ay, (16) 


In (14), (15) und oo sind einstweilen noch die Funktionen iF (s) unbekannt. Zu ihrer Be- 
stimmung multipliziert man Gleichung (9) mit den Funktionen 


=I 
(aes 


72(%,) => = (Mts oo Lars : n) ; (17} 

en (ba — xu) (a2 — xu) 
integriert iber < a, Py, und summiert noch iiber alle Gel ck sie, Wegen der Eigenschaft ® 

b 
po (xu) “ 

x oii =9 fir 68,25, (0 =1,.-n) ee 

verschwinden dann die Glieder mit den J’, und man findet 

2 bn 5 by, xy 
oe mn 2 3 
2 Di 6.) | Vl %,) dx, ey DAG 1) | Vi(%,) ee | WiAlE us s) dé, dx., ae 
of dies ay ay 

byt $ 
a De | E(Sy, s) A,(é;) dé, (19) 
v=1 by 


0 falls &, ) auf einer der 
Strecken 


by 
LS gf ral%) — & IIsign (bx —& 
mit A,(é,) = ea ae —— dx, = | eee ae ee) falls & Cat” (20) 


nicht | (A=. ten) 
sail liegt. 


1 Siehe FuBnote 1 von Seite 180. 
® Vel. K. Nickel, Math. Z. 54 (1951), S. 91, Formel (20). 
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In (19) hat man ein lineares, unhomogenes Gleichungssystem zur Berechnung der Funktionen 
£,As) (u eas vey n) gefunden. Da die zugehérige Determinante nicht verschwindet (vgl. 
. oe > oe ee ist dieses System stets eindeutig aufléshar. Damit kénnen die g(s) 
als bekannte Funktionen angesehen werden, die I i FT A 
ee ernie. nur noch von den Wirbelstarken ¢, der An- 

Zur Bestimmung der ¢,(&,, s) (vy = 1 n) geh i i i 

_ Le Ae, =1,...,m) geht man auf die Gleich 3 : : 

tipliziert man sie mit den Wirbelverteilungen a Oges he 


=i Ox Xp Fea RCA) == 


b 
a se SS 
=o = 2) ay — Xu il = b — & 
Vil Xx) 1 IVE daa -H ‘ ~ dé, , (lu, @ Sal oie 5 1) 25 (21) 
“ 


integriert tiber ¢a,,b, > und summiert noch von yu =1 bis n, so ergibt sich nach einigen 
Umformungen wegen 


Ba 
I~ ¢ il Xu) { 1 fur O=% 
ne Ste Reset ee om 
Pa Xu — Sp b % | 0 fiir QO S= 9) (22) 
iD 


unter Beachtung der Gleichungen (4) und (6) das Gleichungssystem 


s 


J 


0 


Bal iy Ota S 


PGn alee | Seat o) |) Th" By 5 4-0) a 
0 


a 


x=1 


nm ie nm 
= oS) | W(X p> 8) /- ae B(x,)dx, (9 =1,.-.4n). (23) 


Dabei wurde zur Abkiirzung gesetzt 


b 
Q , 
sie. 1 “rb, — Eo 
B,(z) =— — SO Use eg (24) 
“9 


In (23) hat man ein Soninesches Integralgleichungssystem zur Bestimmung der ¢, gefunden. 
Bekanntlich besitzt nun die Soninesche Integralgleichung 


J f(0) K(o — s) do = g(s) (25) 
0 
die Lésung ? 
f(o) = g(0) L(c) + [ Lo — s) dg(s) , (26) 
0 
wobei der lésende Kern L(c) die spezielle Lésung von (25) fiir g(s) =1 ist *. Damit wird 


folgende Methode zur Auflésung von (23) nahegelegt: Man betrachtet von den n gesuchten 
Funktionen ¢, die ersten n — 1 fiir den Augenblick als bekannt und lst (23) nach (24) und (25) 


fiir 90 =n auf. Die Lésung ¢,, die natiirlich noch von den e,(9 = 1,..-,m — 1) abhangen 


wird, setzt man in die n — 1 restlichen Gleichungen (23) ein. Wie man sich leicht tiberzeugt, 
ist dieses neue Integralgleichungssystem immer noch vom Sonineschen Typ, kann also ebenso 
weiterbehandelt werden. Nach n — 1 Schritten ist man auf einer einzigen Sonineschen Inte- 


 gralgleichung, etwa fiir ¢,, angelangt, die man nach (26) lst. Die anderen gesuchten Funk- 


tionen ¢€ bis ¢, bestimmt man genau so. Zwar kommt man dabei auf sehr unhandliche Aus- 


_ driicke, doch treten keine prinzipiellen Schwierigkeiten auf. Ob allerdings (23) stets und ein- 


deutig lésbar ist, ist damit noch nicht gezeigt. 
Dem Gleichungssystem (23) entspricht im Falle n = 1 die Wagnersche Integralgleichung. 
Dieser Zusammenhang 1aBt sich noch etwas deutlicher zeigen. Beachtet man namlich, daB die 


1 Siehe FuBnote 1 von Seite 180. 
2 Das Integral ist als Stieltjes-Integral, aufzufassen. ' ; 
3 Vel. etwa W. Schmeidler, Integralgleichungen mit Anwendungen in Physik und Technik I. S. 216 ff. 
Leipzig 1950. 
LB 
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B, die Eigenschaft 


0 falls =) auf einer der Strecken 
PS B,(z) = 1 — fe by —F falls z « a), b,> (27) 
e=1 a a, —z nicht} (4 —1,...,n) liegt 


besitzen !, so findet man aus (23) durch Aufsummieren iiber alle 9 = 1,...,n die eine Glei- 
chung 


Omi 


die vollkommen der Wagnerschen Integralgleichung entspricht. Nur sind fir n > 1 die Funk- 
tionen ¢, durch diese eine Gleichung (23a) sicher nicht mehr eindeutig bestimmt, dazu ist viel- 
mehr das allgemeine Gleichungssystem (23) erforderlich. 

Im allgemeinen Falle, wenn also die Fligelprofile nicht mehr alle hintereinander angeordnet 
sind und ihre Anzahl u. U. nicht mehr endlich ist, kann man zur Bestimmung der Funktionen 
g,(s) und ¢,(C,,s) genau so verfahren. Man multipliziert (10) mit Funktionen Vil2y) der 
Kigenschaft 


x P78, Re| Ese == 0 (Op sls ceed (18a) 


integriert tiber <a, b, > und summiert iitber alle ~. Damit kommt man auf ein lineares Glei- | 


chungssystem fir die g, in Abhangigkeit von den é,. Zu deren Berechnung sucht man weitere 
Zirkulationen y/(z,), die diesmal dem Gleichungssystem 


b 
p 
1 = 1 
Peay f F21e,) Re |: akan =, (g.%=1,2,---) (22a) 
u a, 


geniigen sollen. Die obige Operation auf (3) angewandt liefert dann ein Soninesches Integral- 
gleichungssystem zur Berechnung der Wirbeldichten En(C y> $)- 


} 


Abb. 4. Senkrechtes Fligelgitter. 


5. Senkrechtes Fliigelgitter. Als weiteres Beispiel werde das senkrechte Fligelgitter gemaB 
Abb. 4 behandelt. Mit der Bezeichnung k =a/t (t = Gitterteilung nach Abb. 4) geht (3) 
tiber in 

b+s 


b 
w(x, 8) = Gree s) tg k(x — £) dé + = e(E, s) Ctg k(x —£) dé. (3a) 
a b 


Der Kirze halber sollen die Umformungen entsprechend zu den Ziffern 3 und 4 ibergangen 
und sofort das Ergebnis angeschrieben werden. 


* K. Nickel, a.a, O. (FuBnote 2 von Seite 184), Seite 90. 


s by - | 
: + by bo es %—%y 7. 9 (93 
i . Eu(Oy, 0) |i pee Soiree cn > frat |/ I rmepa iat 
x= L= a, “n= 
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Die aaa y(§, s) itber die Profiltiefe lautet 


‘Gin k(b — €) 
B=.) Sinkeoay Tt 


a x 
w(x, s) + Sf w(n, s) dy 


2k SS 
+ *yemEe —H Game ay a 
= ( )) Gin k(x —&)/ Gin k(b — x) Gin k(x — a) get (28) 


mit 
b +b 
Coj k (« DR ) 
Qk ; b—a z 
ee ewes" [wie dx — 
Sin k(b — a) ssieah es ( \Smie = a cies a 
o(8) + 190) ky(s) + f ky(s — 0) ary) (29) 
0 
wobei der letzte Summand als Stieltjes-Integral aufzufassen ist. 
Dabei ist 
b 
= Gin k(x — a) 
ier rte s) Ginkb =)" (30) 
und : 
. Cojk(s—o +2— 2) 
k,(s) = &,(b, 0) : do , (31) 
; \ Gin k(s — 0) Gin k(s — 0 + b —a) 
wobei ¢,(, o) die spezielle Lésung der Wagnerschen Integralgleichung 
Gin k(s —o + b —a) i 
fut 0) SERIA a + Ts) = 0 32) 


fir [°,(s) = 1 ist. 
Der Auftrieb eines Gitterprofils berechnet sich zu 
b 


Als) =u. | 116.5) a5 sg ry eee <— Gin k 


3 


a 


a 


° i Sit Ages tad | 
etee. | |e s) a ce s) | Best Prater) sa (33) 


/ Gin k(b — x) Sin k(x — a) 


a 


und das Langsmoment um den Nullpunkt wird 
b 


M(s) =o, [ 7169) £48 =o, lei L Go gees 
tee |laina +2 feeder re pbain 8) 
eee ee ACI : = 3 d MESS ee 

i ae w(x “Del Bs J w(7, s) dy ee eee 

b 
Cop k(x — 2 
ALG C08 dx , (34) 
Cojk ; 


wobei mit Arc cos der Hauptwert bezeichnet werde. 
Diese Formeln entsprechen vollkommen denen fiir den Einzelfligel und gehen fir k — 0, 


also unendlich werdende Gitterteilung, auch in diese iitber?. Da fiir strémungsfeste Bewe- 
- gungen der zweite Summand in (28) verschwindet, erkennt man hier ohne Rechnung die Giltig- 


keit des in der Einleitung formulierten Satzes. Uberdies greift dann noch, wie beim Kinzel- 


-fliigel, nach (33) und (34) der Auftrieb eines Profils aus dem Gitterverband stets im vorderen 


Neutralpunkt an. 
1 Vel. H. Séhngen (FuBnote 6 von Seite 179), doch beachte man FuBnote 1 von Seite 182. 
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Fiir den Anfahrvorgang des senkrechten Plattengitters, wenn also das Gitter aus der Ruhe 
heraus unter dem konstanten Anstellwinkel « mit einer gleichférmigen Geschwindigkeit u,, 
bewegt wird, vereinfachen sich mit 


as ull 0 fur s-=40 
CNS Tu Saini tur sis 0 
die obigen Formeln zu 
fo(s) 2. a a sin o = Gin k —! 
a(s) = 2m u, sind ook 5 
Coj k 5 4 
: : ; ky(s) Gin k(b — &) 
Auftrieb teil me = 
uftriebsverteilung —y(&, s) 2 wu, sin & a ae ene: 
2 
Auftrieb A(s) = 209 ut, sin a hy(s) — Tg k- >" 
und die Lage des Auftriebsmittelpunktes (gleichzeitig die Abszisse des Neutralpunkts) 
M a+b 1 b—a b—a 
er (s) oe 5} — =, ta k 5} log Coj k 5} A 


Auftriebsverteilung y(&,s) und Auftrieb A(s) einer Gitterschaufel gehen also beim Anfahr- 
vorgang durch Multiplikation mit dem ,,Abminderungsfaktor“ k,(s) aus den stationaren 
Werten hervor. 


(Eingegangen am 7. Juli 1954.) 
Anschrift des Verfassers: Dr. Karl Nickel, Instituto Aerotecnico, Cordoba, Argentinien. 
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Drehschwingungen in Kreuzgelenkwellen. 


Von F. Weidenhammer. 


1. Aufgabenstellung. Mit Hilfe eines Kreuzgelenkes (Kardangelenkes) kann man Dreh- 
-momente iiber abgewinkelte Wellen weiterleiten, wovon z. B. der Kraftfahrzeugbau bei seinen 
Konstruktionen haufigen Gebrauch macht. Allerdings muf man bei Verwendung nur eines 
Kreuzgelenkes in Kauf nehmen, da® in der abgewinkelten Abtriebswelle eine periodisch 
schwankende Drehgeschwindigkeit (Winkelgeschwindigkeit) und damit ein ebenfalls schwan- 
‘kendes Moment entsteht, selbst wenn der Antrieb mit konstanter Drehzahl bei konstantem 
Moment erfolgt. Um diese Schwankungen in tragbaren Grenzen zu halten, beschrankt man 
sich in der technischen Praxis auf kleine Ablenkwinkel y < 20°. Trotzdem besteht wegen dieser 
periodischen Schwankungen die Gefahr der Anregung von Drehschwingungen in besonderem 
Ma8e, worauf bereits H. Dietz! hingewiesen hat. Im Folgenden soll daher das Drehschwin- 
gungsproblem fiir eine Welle mit einem Kreuzgelenk allgemein behandelt werden. Wie 
H. Dietz+ ebenfalls gezeigt hat, sind die Wellen auch periodischen Beanspruchungen durch 
Biegemomente unterworfen. Es werde jedoch eine solche Lagerung der Welle vorausgesetzt, 
daB Biegeschwingungen merklicher Amplitude nicht erregt werden und im weiteren auBer An- 
satz bleiben kénnen. 


| 2. Die Schwingungsgleichungen. Wie bei Drehschwingungsrechnungen bewahrt, werden 
' die Wellenstiicke als torsionsfahig, aber masselos angesehen und nur die Drehmassen J, und J, 
an den Wellenenden beriicksichtigt. Ferner sei auch die Masse des Kreuzgelenkes selber ver- 
nachlassigbar klein gegen die Drehmassen 


J, und J, (Abb.1). Durch eine kinematische 1? 


| Abb. 1. Die Kreuzgelenkwelle, gezeichnet fiir ein Gelenk 
mit Querverschieblichkeit und statisch bestimmter Lagerung Abb. 2. Die elastische Verdrehung der Welle und die Momentbelastung 


der Welle. der Wellenstiicke. Mg; steht fir Mg; bzw. Mgi (Gir 05L) 5 


Betrachtung der Ubertragungsverhaltnisse am Kreuzgelenk findet man? fiir die Drehgeschwin- 

digkeit ,(t) auf der Abtriebsseite abhangig von der Drebgeschwindigkeit «,(t) und dem 
Drehwinkel g(t) auf der Antriebsseite (gerechnet von der Schnittgeraden der beiden Ebenen 
aus, in denen sich die Gabelzapfen drehen): 


Oy = —* mm, (0S 74 < 90. (1) 


Td erie in2 
1—sin? x sin? My 


Zur Betrachtung des Drehgleichgewichtes bei stationarem Betrieb der Welle ohne Schwin- 
gungen werden auBer den durchgeleiteten Arbeitsmomenten M, und M, noch die beiderseits 
des Kreuzgelenkes angreifenden Momente My und M,, eingefiihrt (Abb. 2). Dabei kenn- 
zeichnet ein Querstrich bei den Kreuzgelenkmomenten den stationaren Zustand ohne Schwin- 
gungen. Da die Lastmomente M, und M, nur fiir den stationaren Betrieb berechnet werden 
und spater als nicht wesentlich durch Schwingungen geandert mit den gleichen Werten in die 
Schwingungsgleichungen ttbernommen werden, eriibrigt sich hier eine solche Kennzeichnung. 


1H. Dietz, Z. VDI. 82 (1938), S. 825. 
20. Kraemer, Getriebelehre, S. 199, Karlsruhe 1950. 
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/ 
Aus Abb. 2 liest man als Bedingungen fiir Drehgleichgewicht ohne Schwingungen ab: / 


M, — Jo %) — My = 9, | (2) 
M, —J,4,—M, =0. J 
Auf Grund des Arbeitssatzes und der kinematischen Grundgleichung (1) gilt fiir die Gelenk- 


momente M,, und M,): 
Mgo On cos X 


a oa ene ? 
Mey Wo ] —sin? x sin? @ 


(3) 


so daB also ein eingeleitetes konstantes Moment M,, als ein periodisch schwankendes Moment 
My in die Abtriebswelle weitergeleitet wird. Aus den Gleichgewichtsbedingungen (2) und der 
Ubertragungsgleichung (3) folgt dann fiir die durchgeleiteten Arbeitsmomente M, und M, die 
Beziehung 

M,—Jo ®o __ cos X% 

M,+J,0,  1-—sin? ysin? gy ° 


Daraus erkennt man, daB die Welle bei stationarem Betrieb mit gleichférmiger Drehung auf 
der Antriebsseite 

® =0, wW =konst., % = Wt 
und periodisch schwankender Drehgeschwindigkeit w,(t) auf der Abtriebsseite lauft, wenn 
zwischen den Arbeitsmomenten die Beziehung 


— J, % (4) 
besteht. In der fir alle weiteren Betrachtungen als geltend vorausgesetzten Gleichung (4) 


wurde zur Abkirzung 


cos X 
SD) pa acne 
1 — sin? x sin? @,t 


g(t) = 


fiir die in (1) und (3) vorkommende periodische Funktion der Zeit gesetzt. 

Zum stationiren Lauf der Welle mit gegebenen Drehgeschwindigkeiten w) und @,(w ) 
werden nun zusatzliche Drehschwingungen betrachtet, welche als differentiell klein voraus- 
gesetzte elastische Verdrehwinkel #) und #, bei den Drehmassen und als im gleichen Sinne 
klein angesehene Verdrehwinkel },, und 3,, beiderseits des Gelenkes zur Folge haben (Abb. 2). 
Wahrend % und #, unabhangige Verdrehwinkel zur Beschreibung der Drehschwingungen sind, 
besteht zwischen den Winkeln 3, und @,, beiderseitig des Gelenkes zunachst eine durch (1) 
gegebene Bindung. Denn wegen w, = dq/dt, w, = dg,/dt und der vorausgesetzten Kleinheit 
von #,, und #,, verlangt (1) fiir diese Verdrehwinkel die Gleichung 


Boy = E(t) Oy - (6) 
Eine weitere Bindung zwischen 0. und 0, ergibt sich aus der Ubertragungsgleichung im 
Kreuzgelenk fiir die schwingende Welle. Setzt man namlich die bei den Drehschwingungen 


auftretenden elastischen Riickstellmomente den jeweiligen Verdrehwinkeln proportional an, 
so gilt 


mit W) = konst. (5) 


Meo = % (8o— 9—0)> Mg = 6s (On — 14), (7) 
wenn ¢, und c, die Verdrehsteifigkeiten der Wellenstiicke und eine aufgesetzte Tilde den 
Schwingungszustand bezeichnen. An Stelle von (3) gilt jetzt 


Meo =o (9 ee Po) a= Oi + bg 


Me, cy (Og, — 94) @ + Seo ; 


und die Beziehung liefert nach Linearisierung und mit der Abkirzung (5) eine weitere Glei- 
chung zwischen #) und 1: 


Cy (I — Boo) — g(t) cy (Oni Gh: (8). 
Mit Hilfe von (6) und (8) lassen sich #,) und 3, durch # und #, ausdriicken, so da® die Ver- 
drehwinkel 9) und 0, als Schwingungskoordinaten verwendet werden kénnen: 
Lys aon __ Bey Oo + B74 OD 
D0 i, ein ee eats rs 2 = (9), (10) 
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Damit lassen sich die Schwingungsgleichungen fiir Drehschwingungen beim stationdren 
Betrieb der Welle angeben. Zunachst erhalt man an Stelle der Gleichgewichtshedingungen (2) 
nun fiir die schwingende Welle mit ® = 0 


My — Jo 9; — My = 0, Mo Jay Os) =, =—(. 


Mit den Riickstellmomenten (7) und der Gleichgewichtshedingung (4) ergeben sich daraus 
dann die Schwingungsgleichungen 

Jo Fy + & (I — Fe) = My, J, OF +6 (6, — Ig) = — ea (11) 
worin die rechten Seiten bei gegebenen Lastmoment M, bekannte Funktionen der Zeit sind 
und die J, und #,, gemaB (10) durch # und #, ausgedriickt werden kénnen. Da der Ablenk- 
winkel x klein ist, bietet sich hierbei eine Entwicklung nach Potenzen von sin y zur Verein- 
fachung an. Setzt man 

b= Bin? y", 
so ist gemaB (1) 
cos ¥ 


2 
g(t) = =1—* cos 2 wt + 0(k4), 


Sa Geel 
1 —sin’ y sin? wy t 


worin in sehr guter Naherung vierte Potenzen von sin y fir 7 < 20° vernachlassigbar sind. 
Mit dieser Entwicklung und den Abkiirzungen 
al 


arr a cae 


fir die Verhaltnisse der Torsionssteifigkeiten vereinfachen sich die Gleichungen (9), (10) zu 
Bao = Yo (1 + k* y, cos 27) By + n(l +h B af zi cos 2 1) 0, 4 0(K4), 


Ba = Yo(1 + Bly, — F| 008 22) By + 74 (1 +H [yy — U] cos 2 2) By + 0184), 


worin noch durch 

T=, t 
eine dimensionslose Zeitzahlung eingefiihrt ist. Da die Antriebsdrehzahl konstant voraus- 
gesetzt wurde, lauten die Schwingungsgleichungen (11) ausfiihrlich geschrieben, 


# il 
§ Jo By + [e0 (1 — 70) — keg Yo 71 Cos 27] Dy — k NAR ey (1 oa 4 cos 2 (| = val 
(12) 


ws J; ay = [ce (1 — 7;) — k? Cy V1 ("4 — 1) cos 2 T] 0, 
AVE +a ro(n = 54 cos 2] o = —M— M, * cos 21, 


wenn Punkte nun Ableitungen nach t = w,t bezeichnen und weiterhin Potenzen k* = sin‘ 
und héhere unberiicksichtigt bleiben. Das Gleichungssystem (12) beschreibt erzwungene 
rheolineare Drehschwingungen mit zwei Freiheitsgraden, obgleich nur zwei Drehmassen be- 
ricksichtigt sind. Dieser zunichst unerwartete Sachverhalt fiihrt in der Schwingungsrechnung 
zu Besonderheiten, und es ist daher zweckmabig, die Schwingungsgleichungen (12) durch Ein- 
fihrung von Hauptkoordinaten zu vereinfachen. 


3. Einfiihrung von Hauptschwingungskoordinaten. Setzt man in den Schwingungsglei- 
chungen (12) y = 0, so miissen sich die bekannten Beziehungen eines Zweimassen-Dreh- 
schwingers mit nicht abgewinkelter Welle ergeben. In diesem Grenzfall (k = sin y = 0) liefert 
dann der Drallsatz fiir das vorliegende Schwingungsproblem Jy %) + J, 0, = 0 und gestattet 
so, die Zahl der Freiheitsgrade auf einen zu verringern, indem man die Differenz der Verdreh- 
winkel 9, — % als neue Schwingungskoordinate einfihrt. Die hierdurch auch fiir den allge- 
meinen Fall (k = sin y + 0) bedingten Vereinfachungeri nutzt die nachfolgende Transfor- 
mation auf Hauptkoordinaten aus. Mit der Abkiirzung 

Pe Ceti (13) 


Cy + Cy 
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und den Matrizen 


1 
ip 0 ¢ = € Cc Yo Cc (7: =< (14), 
Mt od r YA, — my 1 = i] (15), 
Dae: mee ne ‘ (v ve | pages (16) 


und den Vektoren 
M. 
ry) (r) 0 A 
d=|( ° : = k 1718 
ea) t (e M,(1 +* cos 2«)] Ce 


lassen sich die Schwingungsgleichungen (12) zu 
wr Md°+ Wd — kcos27- Wd =x (19) 


zusammenfassen. Das zugehdrige kinetische Potential ist als Differenz der kinetischen Energie 
und der potentiellen Energie mit (14) bis (18) durch 


L=—%omd +49 %d — = cos2e-d’ Md —d't (20) 


gegeben, wenn ein Akzent den Ubergang zur gespiegelten Matrix bezeichnet. Die Einfihrung 
von Hauptkoordinaten verlangt nun die Einfiithrung solcher Koordinaten, daf die ersten 
beiden quadratischen Formen in (20) rein quadratisch in den Verdrehgeschwindigkeiten und 
den Verdrehwinkeln werden. Da { eine Diagonalmatrix ist, fihrt dies auf eine Matrizen- 
eigenwertaufgabe der Zwischenstufe 1. Man fihrt in solchen Fallen die Hauptkoordinaten 
X(t), x(t) durch die Transformation 


b=MA2 Uy (21) 


hk X(T) 22 
: ee (22) 


we — Fiat 3) 
0 Uys, 


sowie eine Transformationsmatrix |] vorkommen. Zur Bestimmung der Matrix U hat man 


eine spezielle Eigenwertaufgabe fiir die charakteristischen Zahlen x und die Eigenvektoren } 
zu lésen: 


ein, worin der Vektor 


und eine Diagonalmatrix 


(23) 


ME EN te Me zea 0 (24) 
Man findet zunachst die Eigenwerte aus der zu (24) gehérenden Determinante 
Sy — 7? poet 
Jo” elie J 
PE AY (25) 
SA OAL eee iy 2 
VJo Ji Jt 
mit den Werten 
= anwar! 
me Vc ae yh (26) 


und kann dann die orthogonalen (t); = 0) und normierten (0500 = 1, 9) 9, = 1) Lésungs- 
vektoren {jy und t), durch Auflésung von (24) leicht berechnen. Damit lat sich die in (21) 
vorkommende orthogonale Transformationsmatrix 11 zusammenstellen: 


eA 1 Aas 
ei ae BV ree 

und es gilt UW’ U = € und Ut = UW. Mit den Werten von (22), (23) und (26) wird das kine- 

tische Potential (20) durch die Transformation (21) in die einfachere Form 

Oon 5 Uo ke , , 

Dag EE a ak RD Or ag CORD Ee eee Ih (28) 


* Vel. L. Collatz, Kigenwertaufgaben, S, 253, Leipzig 1949. 
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ibergefiihrt. In (28) kommen die symmetrischen Matrizen 


x2 0 
wy =(° 
f ) (29) 
und 
Sau a me n= (* %) 
C10 C1 
fer (e es J V5oJs [Y Lyd \ 
eg Wd ; Daler d ee (> a )( ee 2 v1 
—7 ( Joe a J; nds Jo tJs Jo mate 2 Jo J as | (30) 
{ 
Tas i 1 1 VA Jo Yo Bes, V1 
Mal So ae ora 
Gee |. 2 wa J J; aka Ne Mp Seg 
vor. Ferner bezeichnet ${ in (28) den Vektor 
Ry M, ( es Var 3 cos 2 7) 
R, M, *— cos 2T 
eee 2 
Das kinetische Potential (28) liefert nun iiber die Lagrangeschen Gleichungen 
anol OL : 
ea ie i ase) 


die Schwingungsgleichungen (12) in den Hauptkoordinaten: 


A Jona : 
Mex, + x2 x — k* cos 2 t (eg X% + Co 1) = M A Joh Va a a Ecos 2), és 
2 te =) fe2 —_ ke 
Oo? %, k? cos 2 (¢yp % + Cy, x1) = —M, ae cos 27. 


Die freien Schwingungen derartiger rheolinearer Schwinger sind in anderem Zusammenhang 
von E. Mettler 1 untersucht worden, und auch gewisse Resonanzeigenschaften der erzwungenen 
Schwingungen sind bekannt ®. Allerdings liegt hier insofern ein Sonderfall vor, als eine ,,Kigen- 
frequenz‘* verschwindet. 
Um die mechanische Bedeutung dieser Entartung hervortreten zu lassen, werde nochmals 
der Grenzfall der nichtabgewinkelten Welle betrachtet: y = 0 d.h. k = 0. Dann werden die 
Schwingungsgleichungen (32) sehr einfach: 


or (0) Ss ae a (0). 
ws Xo ae Hoe x= M, 2 . ws “= 0, (33) 
(0) 


wobei der Index oben Null den Grenzfall k = 0 bezeichnet. Um die Hauptkoordinaten x), 
(0) (9) (0) 
x, in den urspriinglichen Koordinaten, den Verdrehwinkeln 0, 0,, auszudriicken, wird die 


Umkehrung der Transformation (21) bendtigt: 


p= (%) — wm = —1 _ (I Jo Si % — Vo 1), 
4 VJo+ Ji Jeera vd 1 
Damit erkennt man die Gleichungen (33) als eine Drehschwingungsgleichung mit der Eigen- 


kreisfrequenz x, fiir den einzigen bei nicht abgewinkelter Welle noch verbleibenden Freiheits- 


grad: 


: (0) (0) : (0) (0) Jot St 
(3, — B,) + 5 (Fo — A1) = My Jo Ji 
und eine Aussage 
rae 1 (0) (Q. 


Tag, Jo % +591) = 9, 


1 E. Mettler, Ing.-Arch. 17 (1949) S. 418. 
2 F,. Weidenhammer, Z. angew. Math. Mech. 32 (1952) S. 154. 
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0 0 
die auf den Drallsatz Jo i, + J, , = 0 hin ausgewertet werden kann. Fir die abgewinkelte 
Kreuzgelenkwelle (k + 0) gilt offenbar die letzte Beziehung nicht mehr. Die vektoriellen Aus- 
sagen dieses Satzes verlangen dann die Einbeziehung auch der Momente der periodisch schwan- 
kenden Lagerbelastungen und fihren also nicht zu hier allein interessierenden Aussagen tiber 
das Drehgleichgewicht der reibungslos rotierend gedachten Kreuzgelenkwelle. Auf Grund 


der eigentiimlichen Ubertragungsverhdltnisse im Kreuzgelenk (M,o/M,, + 1) erfolgen die 
Drehschwingungen mit zwei Freiheitsgraden. Tatsachlich wiirde die Mitbericksichtigung 
vierter Potenzen von k = sin y in den Schwingungsgleichungen (12) auch eine zweite Higen- 
kreisfrequenz x? = 0(k*) liefern. Da jedoch bei den vorliegenden technischen Aufgabenstel- 
lungen ¥ < 20° ist, sollen Potenzen k* = sin 7! als vernachlassigbar klein auBer Ansatz bleiben, 
so daB x7 = 0 gesetzt werden muf. 

In den Schwingungsgleichungen (32) miissen die Zwangsglieder in bezug auf ihre Periodi- 
zitat betrachtet werden, wenn man das Resonanzverhalten verstehen will. Im allgemeinen 
(z. B. bei Antrieb durch Kolbenmotore) wird das durchgeleitete Arbeitsmoment M, nicht 
konstant sein, sondern mit dem Drehwinkel periodisch schwanken. Gilt also die Fourier- 
zerlegung 


M,(t) = & (a, cosrtv + b, sin rT) (34) 
r=0 
mit bekannten Entwicklungskoeffizienten a,, b,, so sind die Erregerglieder Rj und R, durch 
(31) ebenfalls leicht als Fouriersche Reihen vom gleichen Typ anzugeben: 


Ryasd) (on, cos r 0 85 sin Tea (t==0, 1) (35) 
=0 


worin die Koeffizienten «;,, und §;,, im weiteren jedoch zahlenmafig nicht bekannt zu sein 
brauchen. In den Hauptschwingungsgleichungen (32) kommen demnach als Erregerglieder 
im allgemeinen trigonometrische Reihen mit der Periode 2 vor. Wenn jedoch M, spezieller 
zusammengesetzt ist, kann dies bei rheolinearen Schwingern wesentlichen EinfluB auf die Lage 
der Resonanzen haben. Bei den vorliegenden Aufgabenstellungen diirfte jedoch nur der Fall 
eines konstanten Drehmomentes Mo, wie es etwa bei elektrischem Antrieb gegeben ist, Be- 
deutung haben. Fir M, = konst. sind die Erregerglieder (35) unmittelbar durch (31) selbst 
gegeben. Ry) und R, sind demnach in diesem Falle gerade Funktionen der Zeit mit der kiirzeren 
Periode z, was ein Absinken der Resonanzdrehzahlen auf kleinere Werte zur Folge hat, wie 
sich zeigen wird. 


A. Die Resonanzfrequenzen. Zur Auffindung der resonanzkritischen Drehgeschwindig- 
keiten (Resonanzfrequenzen) Wy = a) der Antriebswelle sollen nun die durch (32) beschriebenen 
erzwungenen Schwingungen berechnet werden. Zu dem Zweck werden die Fourieransatze 


x; = DS, (Az, cos rr + B,; sind). (60,1) (36) 
r=0 


in (32) eingesetzt gedacht und der Ablenkwinkel y durch sin y = k als Kleiner Parameter k in 
die Rechnung aufgenommen. Das dann entstehende unendliche Gleichungssystem fiir die 
Fourierkoeffizienten A; und B;,, braucht im einzelnen nicht niedergeschrieben werden}. Es 
geniigt lediglich zu erkennen, daB dieses System in vier voneinander unabhangige Teilsysteme 
zerfallt. Denn auf Grund der in (32) allein auftretenden Produkte 


cos 21+ A; cos rt = 55 (cos (r + 2) 7 + cos (*—2)1), | 
Bi,r (37) 
cos 2+ B,, sinr y= —~ (sin (r + 2) t + sin (r — 2) t) | 


kann erstens keine Mischung der cos- und sin-Ansatze von (36) eintreten, so daB also die 
Koeffizientengleichungen entweder die A,,, oder die B;,, je fir sich allein enthalten. Zweitens 
bedingen die Identitaten (37) auch noch, daf in einer Koeffizientengleichung jeweils nur A,,, 
oder B;, mit geradzahligem oder ungeradzahligem Zweitindex vorkommen kénnen. Demnach 


verbleiben vier voneinander unabhingige inhomogene Gleichungssysteme zur Berechnung 


1 Vel. Formel (9) und (10) der in FuBnote 2 auf Seite 193 genannten Arbeit. 
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der Fourierkoeffizienten aus den Ansatzen (36). Die jeweils unendlich viele Unbekannte ent- 
haltenden Systeme seien durch 


Are pte eas 9s) — Oe a 1 Bay 4 1) > (38a) 
(Ao,2, Ai. oy) == (00,27 ~ &1,27)5 (38b) 
(Soran es Bieri) = (fo 97-- 1 Pdr) (38c) 
(Boar (Bia =(Po2+ Bier) (38d) 


bezeichnet und enthalten als inhomogene Seiten die Fourierkoeffizienten der Zwangsglieder (35). 
Die Auflésung der Gleichungen ist méglich, sofern keine der (konvergenten) unendlichen 
Koeffizientendeterminanten fir A; und B;,, verschwindet. Ist dies jedoch fiir gewisse Para- 
meterwerte @, und k bei auch nur einer Determinante der Fall, so liegt ein Resonanzfall vor, 
denn erzwungene Schwingungen endlicher Amplituden lassen sich nicht mehr durch den 
Fourieransatz (36) berechnen. Allerdings tritt dieser Resonanzfall nur ein, wenn mindestens 
ein Fourierkoeffizient der Erregung auf der inhomogenen Gleichungsseite dieses Teilsystems 
von Null verschieden ist. Im Falle des konstanten Drehmomentes M, hat diese Einschrankung 
eine Verminderung der Resonanzgefahr zur Folge, wie sich ergeben wird. 

Zur Berechnung der resonanzkritischen Werte ao) = w)(k) kénnte man nun die Deter- 
minantengleichungen zu den Gleichungen (38) heranziehen und oy aus diesen berechnen: 


Big ere Ais pat) = 0% (39a) 
A(Ac2, “Ai2,) =; (39b) 
A(Bo2r +1 Bior41) =9, (39c) 
A(Bo,2, —Bi,21) = 0. (39d) 


Doch kann man die Werte a, tiber die freien Schwingungen wesentlich einfacher berechnen. 
Offenbar besagt das Verschwinden einer der Koeffizientendeterminanten (39), daB die zuge- 
hérigen homogenen Gleichungen eine bestimmte periodische Lisung besitzen, da ja ihre 
Fourierkoeffizienten aus den homogenen Gleichungen von Null verschieden zu berechnen sind. 
Uber diese freien Schwingungen kann man auf Grund der Determinantengleichungen noch 
weitergehende Aussagen machen. So beschreiben (39a) und (39c) Schwingungen mit der 
Periode 2 2, wobei (39a) nur cos-Glieder und (39c) nur sin-Glieder enthalt. Dagegen fihren 
(39b) und (39d) auf freie Schwingungen mit der kleineren Periode 2, wobei (39b) wiederum nur 
cos-Summanden und (39d) nur sin-Glieder enthalt. 

Nachdem auf diese Weise die Zusammensetzung der freien Schwingungen weitgehend be- 
kannt ist, kann man diese nun mit einer kurzen Stérungsrechnung sehr einfach berechnen und 
die zugehérigen Frequenzparameter o,(k) angeben. Diese Eigenfrequenzwerte a) sind dann 
die gesuchten resonanzkritischen Drehgeschwindigkeiten. Zur Auffindung dieser Werte 
Wy = % werden sowohl die Frequenzparameter wie auch die periodischen Lésungen der 
Schwingungsgleichung nach dem kleineu Parameter k entwickelt: 


(0) (2) 
OK = Oy + kh? aw, + 0(k*), (40) 
(0) (2) ; 
ee he at Oks ae (e081). (41) 
Fiihrt man diese Ansatze in die homogen gemachten Gleichungen (32) (M) = 0) ein, so erhalt 
man zur Berechnung der freien Schwingungen nach einem Vergleich gleichhoher Potenzen von 
k die rekursibel auflésbaren Differentialgleichungen 


(0), (0 (0) (0). (0) 


) oe 
AG X + 5% = 0, a9 x, = 0; (42a) 
(0) (2), (2) (0) (0) (0) (2) (0), 
OG X) + % % = cos2Tt (Coo Xo + £01 x;) — 2 Oy % %X, > (42b) 
(0) (2), (0) (0) (0) (2) (0), | 
a8 w= COB 2 v (C49 % + C11 1) — 2 My oq X13 °° *- 


Durch Auflésung der inhomogenen Differentialgleichungen sollen nacheinander die periodi- 
schen Lisungen berechnet werden, deren Vorhandensein aus den Gleichungen (39) abzulesen 


waren. 


196 Weidenhammer: Drehschwingungen in Kreuzgelenkwellen. Ingenieur-Archiv 


(0) 
a) cos (2 r + 1) t-Typ entsprechend (39a). Als Ausgangslisung ist hier x) = cos (2r + 1)t 
(r =0,1,...) zunehmen. Dann folgt aus den Gleichungen (42a) zunichst 


(0) Ko 


( : 
und weiter muf Hs = 0 sein, da eine periodische Partikularlésung aus den Gleichungen (42) 
berechnet werden soll. Die erste Gleichung (42b) lautet mit den gefundenen Werten dann 
2 (2). , (2) i (0) (2) 
(a7 bap L 42 Xe =" (cos (2 r + 3) 7 4- cos(2 r — 1) 7) + 2 M% 0% (2 r + 1)? cos (27 + 1)T, 
(2) (2) 
und man erkennt, da® x) nur dann frei von nichtperiodischen Gliedern bleibt, wenn man ag 


speziell wahlt: 


(2) Cee (2) , 
ica a Leet r= 0» nnd. = 08 tin fe an 
sa) 


Damit sind folgende Resonanzfrequenzen gefunden: 


oy = Hy — HLH + Of) (r = 0), | (43) 


aie 4 = CP Ls . 
Faeroe seme 7 aa) (44) 
Auf Grund von (38a) sind diese Werte dann resonanzgefahrlich, wenn das durchgeleitete 
Arbeitsmoment Fourieranteile vom cos (2 r + 1) t-Typ (Periode 2 z) enthalt. 
(0) 
b) sin(2r-+1)z-Typ entsprechend (39c). Die Ausgangslésung x) = sin (2r-+1)t 
(r = 0,1,2,...) fiihrt nach einer entsprechenden Rechnung wie in a) auf die Resonanz- 
frequenzen 


Op = 


X = Xo ai Ke om a O(k*) (r ar 0) z ) 
oy = et OR) H1,2,.. 4, (46) 


die gema8B (38c) resonanzgefahrlich sind, wenn in der Erregung Anteile sin (2 r + 1) t (Periode 
27) vorkommen. 


(0) 
c) cos2rt-Typ entsprechend (39b). Als Ausgangslésung muf hier « = cos2rrt (r = 1, 


2,...) genommen werden. Dann liefert (42a) zunachst 
es 
OX as or 


(0) 
und wieder x, = 0, wahrend aus (42b) 
0) (2) 


%o\2 (), 9 (2) Cee ( 1 
i Ky + %o Xo = —" (cos (2r + 2)t + cos (2r — 2) 7) + 2 9 Mp Gg cos Art 


(2) 
folgt. Also mu o) = 0 sein, damit auch x, rein periodisch wird. Es gilt also hier 
My = 52 + O(h4) (r = 1,2,...), (47) 


und diese Frequenzen sind resonanzgefahrlich, wenn cos 2 r 7-Anteile (Periode z) in den Er- 
regergliedern Ry und R, enthalten sind. Dies ist nicht nur der Fall, wenn das Arbeitsmoment 
M, derartige Anteile enthalt, sondern auch schon dann, wenn M, = konst. ist und nur durch 
den. Kreuzgelenkmechanismus eine Schwingungserregung gegeben ist. 


d) sin 2 rt-Typ entsprechend (39d). Von der Ausgangslésung x ='sin 27? (P=-1, 250) 
ausgehend erhalt man hier auf gleichem Wege wie in c) die Frequenzbeziehung 
x 
% = me + 0(k4) (48) 
als resonanzgefahrlich, wenn sin 2 r t-Anteile in My (und damit auch in R, und R,) vorkommen. 


In der vorstehenden Stérungsrechnung sind genau wie in den mechanischen Ausgangs- 
gleichungen Potenzen k* = sin* y und hohere als klein vernachlassigt. Nur auf Grund dieser 
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Vernachlassigungen fallen die Frequenzgleichungen (44) mit (46) und (47) mit (48) zusammen. 
Kine ausfiihrliche Rechnung wiirde ergeben, daB diese Frequenzkurven o,(k) lediglich fiir 
k = 0 eine Berithrung héherer Ordnung eingehen. 

Die resonanzkritischen Drehgeschwindigkeiten der Antriebsseite j) = (Kk) sind in ihrer 
Lage zueinander aus Abb. 3 fiir cy) > 0, k4 = 0 zuentnehmen. Die Resonanzen bei den Werten 
w/t (r = 2,3, 4,...) nennt man Subresonanzen; sie sind bis r = 5 in das Diagramm aufge- 
nommen. Ferner bezeichnen die Abkirzungen s, die sinrt-Typen von (39c) und (39d), 
wahrend c, fiir die cos rt-Typen von (39a) und (39b) steht. 
Man kann so die in a) bis d) diskutierten Resonanzmiglich- 
keiten ablesen. So sind fir konstantes M, alle c, mit ge- 
radem r kritisch, also alle Frequenzen %,/2 r(r = 1, 2,...), 
nicht jedoch z. B. x) selbst. In diesem Falle haben also allein 
die Subresonanzen Bedeutung. Dies Verhalten ist durch die 
periodischen Koeffizienten in den Schwingungsgleichungen 
(32) bedingt und steht in bemerkenswertem Gegensatz zu HH) Hy Ky Xp 
den Drehschwingungen, die sich durch lineare Differential- AS Ae ¥: f 
gleichungen mit konstanten Koeffizienten beschreiben lassen, Pe a a 
bei welchen es keine Subresonanzen gibt. 

Die bei rheolinearen Schwingern méglichen instabilen freien Schwingungen haben bei den 
Kreuzgelenkwellen wohl keine Bedeutung, da man bei diesen Wellen auf Grund der starken 
Lagerbelastungen vermutlich stets eine Dimpfung haben diirfte, welche freie Schwingungen 
schnell abklingen laft. Im Gegensatz hierzu behalten bekanntlich die erzwungenen Schwin- 
gungen auch bei Anwesenheit dimpfender Krafte ihre Bedeutung, weswegen sie ausfihrlich 
untersucht wurden. 


ke 


5. Zusammenfassung. Die Schwingungsgleichungen zur Berechnung der Drehschwin- 
gungen in einer abgewinkelten Kreuzgelenkwelle (Ablenkwinkel 7) wurden bei Beriicksich- 
tigung je einer Drehmasse auf der Antriebs- und Abtriebsseite aufgestellt. Sie beschreiben auf 
Grund der Ubertragungsverhaltnisse im Gelenk erzwungene rheolineare Schwingungen mit 
zwei Freiheitsgraden. Dabei entsteht eine Schwingungserregung auch dann schon, wenn das 
durchgeleitete Arbeitsmoment M, keinen Schwankungen unterworfen, sondern streng konstant 
ist. Die resonanzgefahrlichen Drehgeschwindigkeiten wy =a, der Antriebswelle wurde bei 
Vernachlassigung von sin* y fiir kleine Ablenkwinkel (etwa < 20°) berechnet. Dabei erweisen 
sich im allgemeinen Falle (M) + konst.) die Werte 


% = % + sin? ¥ oy O(sin* y) 
A Xo 
und die Subresonanzen 


te == ir --O(cint-y) (eee se ae 


als resonanzgefahrlich, wenn x, die Eigenkreisfrequenz (26) und cy) den durch (30) gegebenen 
Wert bezeichnen. Fiir den Sonderfall eines streng gleichférmigen Arbeitsmomentes (M) = 
konst.) bleiben nur die Subresonanzen 


% = 52 + 0(sint x) (ea ee Sica 


ibrig, da die Resonanzmdglichkeiten der rheolinearen Schwinger von der Zusammensetzung 
der Erregung abhangen. 


Aus dem Institut fiir Mechanische Schwingungstechnik der Technischen Hochschule Karlsruhe. 
(Eingegangen am 27. Juli 1954.) 
Anschrift des Verfassers: Priv.-Doz. Dr. F. Weidenhammer, Karlsruhe i. B., Hertzstr. 16 Bau 40. 
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Asymptotische Absaugegrenzschichten an langsangestromten 
zylindrischen Korpern*. 


Von W. Wuest. 


1. Einfiihrung. Ein besonders einfacher Sonderfall einer Absaugegrenzschicht tritt an einer 
unendlich ausgedehnten Wand mit gleichférmig verteilter Absaugegeschwindigkeit auf. Das 
Geschwindigkeitsprofil u der Grenzschicht ist an jeder Stelle gleich und nur von der Koordi- 
nate y senkrecht zur Wand abhangig entsprechend der Beziehung 

_ %y 
a =l]—e ’ 
(U bezeichnet dabei die konstante Geschwindigkeit auBerhalb der Grenzschicht und vw) die 
konstante Absaugegeschwindigkeit.) Eine unendlich ausgedehnte ebene Wand ist freilich nur 
ein mathematischer Idealfall. Bei einer Platte endlicher Ausdehnung mit gleichmabig verteilter 
Absaugegeschwindigkeit tritt das genannte Geschwindigkeitsprofil als asymptotische Lésung 
fiir geniigend groBen Abstand von der Vorderkante auf. 

Es liegt nahe, allgemeinere Falle solcher Absaugegrenzschichten zu untersuchen. Eine 
mégliche Verallgemeinerung besteht darin, langsangestrémte zylindrische Flachen in Betracht 
zu ziehen. Wenn man sich auch hier 
wieder auf die asymptotischen Lé- 
sungen beschrankt, gelangt man zum 
Grenzfall unendlich langer zylindri- 
scher Kérper, die in Achsrichtung an- 
gestrémt werden und bei denen die 
Verteilung der Absaugegeschwindig- 
keit nicht von der x- Koordinate inAchs- 
richtung abhiangt. Offenbar miissen 
in diesem Fall alle Geschwindigkeits- 
komponenten unabhangig von x sein, 
so daB in den Navier-Stokesschen Diffe- 
Abb. 1. x, 7, ¢-Koordinatensystem, dargestellt am Beispiel des elliptischen rentialgleichungen alle Ableitungen 

Se nach x verschwinden. 

Es ist fiir das Folgende zweckmafig, von den rechtwinkligen x, y, z-Koordinaten zu x, 7, 
¢-Koordinaten iiberzugehen, wobei 7, ¢ krummlinige, isometrische Koordinaten sind. Man 
kann diese durch konforme Abbildung eines Quadratnetzes gewinnen, wobei wir den ortlich 
veradnderlichen MafSstabsfaktor mit H bezeichnen wollen. Das Koordinatensystem soll dabei 
so gewahlt werden, daB die Begrenzung des Kérpers durch 7 = y) = konst. gegeben ist. Abb. 1 
zeigt diese Festlegungen am Beispiel des elliptischen Zylinders. Die Kontinuitatsgleichung 
und die drei Bewegungsgleichungen lauten dann, wenn u, v, w die Geschwindigkeitskom- 
ponenten in den x, 7, ¢-Richtungen bezeichnen, 


ea 


cae 4 am ao (1) 
oH y+ 0H we tot oe Og ao le | gee pee meee ope 


* Auszugsweise vorgetragen auf der Tagung der Gesellschaft fiir angewandte Mathematik und Me- 
chanik in Miinchen am 21. 4. 1954. 
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Die Geschwindigkeitskomponente u kommt nur in (2) vor; v, w kénnen also unabhingig von u 

berechnet werden. Doch kann aus einer fiir v, w gewonnenen Lisung noch nicht gefolgert 

werden, da8 man auch fiir u eine Lésung erhalt, welche die Randbedingungen befriedigt. 
Fihrt man noch die x-Komponente der Rotation des Geschwindigkeitsvektors 


1 Lee zw pe 


rot =20= 
2D Pl By ae 


(5) 


ein, so kénnen (3) und (4) nach Elimination des Druckes zusammengefaBt werden in 


ow fala) Welle) Fok) 
oH wi — e 
on le oc 4 on? 7 | (6) 
Wir kénnen offenbar die Falle unterscheiden: 
a)w=0 Querstrémung ist wirbelfrei, also eine Potentialstrémung, 


b) w =konst. im Unendlichen dreht sich die Fliissigkeit als starrer Kérper, 

) = @(n). Voder oy = w(C), 

Pro aotnst) fate 0 Hamelsche spiralférmige Strémungen als Querbewegung, 
w + 0 allgemeiner Fall. 


2, Absaugegrenzschichten mit wirbelfreier Querbewegung. Da die Querbewegung (v, w-Kom- 
ponenten der dreidimensionalen Strémung) eine Potentialbewegung! ist, kénnen wir das 
Koordinatensystem so wahlen, daf w = 0 ist. In (1) bis (6) fallen dann die unterstrichenen 
GréBen weg. Aus der Kontinuitatsgleichung (1) und der Beziehung (5) folgt mit w = 0: 


He Kons (7) 


Im Fall rechtwinkliger Koordinaten ist H = 1 und demnach v = konst. Dies entspricht aber 
der homogenen Absaugegrenzschicht an einer unendlich ausgedehnten ebenen Wand und wir 
kénnen daher die Bedingung (7) als eine Verallgemeinerung der homogenen Absaugegrenz- 
schicht auf beliebige zylindrische Kérper auffassen. 

Aus v H = F = konst. folgt 


» oe _ st oH 
an Tei) 
so da (3) und (4) geschrieben werden kénnen: 
eS ot oP HE SL) ee ee oP 
DW on = Q On ‘ oy rela Pai o a 


mit der Lésung 
Q F 1 Q 
I Po 2 ; ie Po Z _ 


Fir die Geschwindigkeitskomponente u erhalt man die Differentialgleichung: 


Ou Cu Cw 
a én “(Fe + A : (8) 
die man durch den Ansatz 
) u = C, 7 +62 + C, (9) 


| lésen kann. Durch Einsetzen in (5) folgt 
Fo =v (x? + B*). (10) 


| Wenn die Wandkurve durch a7 + fC = konst. gegeben ist, kann mit diesem Ansatz also die 
- Randbedingung u = 0 an der Wand und u = U auferhalb der Grenzschicht befriedigt werden. 


3. Beispiele. a) Ruhender Kreiszylinder. ZweckmaBig wahlen wir hier 7 = In'r, Hz, 
F = — a; dabei ist a der Halbmesser des Kreiszylinders und v) die konstante Absaugege- 
schwindigkeit an seiner Wandung. Da die Wandung durch 7, = In a = konst. dargestellt ist, 
- setzen wir in (6) 8 = 0 und erhaltena = — va/y. Die Geschwindigkeitsverteilung der u-Kom- 


1 Die Potentialstrémung ist bekanntlich eine Lésung der Navier-Stokesschen Differentialgleichun- 
gen. Doch erfiillt die Potentialstrémung im allgemeinen nicht die zu stellenden Randbedingungen. 


| 14 
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ponente ist also bei Beachtung der Randbedingungen 


é ey ous A 
et ee ay, 
| 


a 


In Abb. 2 ist die Geschwindigkeitsverteilung fiir verschiedene Werte von vw a/y aufgetragen. 
Das je Einheit der Zylinderlange | 

verdrangte Flissigkeitsvolumen ist | 

gegeben durch 


(12) 


Wenn wir uns das verdrangte Flissig- | 


he ‘ ; EAerte keitsvolumen auf einen Ringraum mit 

Abb, 2. Geschwindigkeitsprofile der u-Komponente an einem Kreiszylinder ; 
bei gleichmaBiger Absaugung. Im schraffierten Bereich haben die Gee dem AuSenhalbmesser a, konzentriert | 
Bob uindis Reltepeotue made Uiee ete nes tan heeee core denken, kénnen wir eine Verdrangungs- | 
| 


dicke der Grenzschicht 


6=a,—a 


definieren. Nun ist aber 


2 71a" 


AQ =x (a2 — a?) 3 
o so dak i 
i 
pana (13} 


ak yes fiir > 2e 
ren In Abb. 2 sind diejenigen Geschwindig- 
keitsprofile, die zu unendlich groBer 


Verdrangungsdicke fiihren, schraffiert. 
b) Rotierender Kreis-| 
zylinder. Wir setzen hier 
n til =(1 +iy) n(x +7y) 
= (1 -+721y)Inre?, 
Abb. 3. Koordinatennetz bei einem rotierenden Kreiszylinder mi eich- 
: EEBleoe Absadpang ee ns) ee iar £0 dab | 
y=Inr—yO, €=ylInr+ @ (14)) 
und He 2S 1 
(tty S| 
wird. Das Stromliniennetz, das logarithmischen Spiralen entspricht, ist in Abb. 3 aufgetragen 
firy = 2/4. Fir F = v Hist an der Zylinderwandung offenbar (vgl. Abb. 3) 


Pee pa iee 
und demnach F = v, a. Die Zylinderwand ist durch r = a gegeben, so daB die Lésung (6) 
u = konst, e(* + #) nr + (6 —a © 4 konst. 
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nur dann mit den Randbedingungen vereinbar ist, wenn 


poy =0 
ist. Dann ist aber wegen (7) 
Fy (1k?) 0 5 


und die u-Komponente ist gegeben durch 


U = (16) 
Die u-Verteilung wird also durch die Rotation gar nicht beeinfluBt. 
ec) Elliptischer Zylinder. Wir wahien hier 
x =co\j7cosl, y =c Ginn sine, (17) 
ff: Coj 2 —cos 2¢ 
Hic / 5 : (18) 


Wenn 7 = 7 die elliptische Wandung darstellt, mu also wegen v H = F = konst. die Ab- 
saugegeschwindigkeit verteilt sein wie 


F (See 
vy =— 3 
c 2 


(19) 


Nun ist aber der Kriimmungshalbmesser ¢ einer Ellipsenkurve 7 = 7) im Punkt ¢, wenn a und b 
die groBe und die kleine Halbachse bezeichnen, 


0 = a? b2 ie ae aM 


a* b4 


g ; ses cos? ¢ _sin?® ¢\3/? 


Lees 1 Coj 27 — cos 20\3/2 HE 
Gin No Coj 7% 2 tab 


Demnach muB8 die Absaugegschwindigkeit verteilt sein wie: 


(or Poe > (19a) 


Ve ab 
also umgekehrt proportional der dritten Wurzel des Kriimmungshalbmessers und nicht mehr 
gleichmaBig wie in den vorhergegangenen Fallen, Die Gesamtabsaugemenge ist gegeben durch 


Ona fod =20F. 


Die u-Verteilung ist offenbar bestimmt durch 


= (in — ny) 
eee a (20) 
und das gesamte verdrangte Flissigkeitsvolumen je Langeneinheit des Zylinders ist 
AQ =| fe sv wr) Het a. 
No 5% 


Die Ausfiihrung der Integration ergibt 


AQ =z 


Co 2 7 (F/v) + 2 Gin 2 mo 
(F/v)? — 4 ; 


Da die Linien konstanter u-Geschwindigkeit konfokale Ellipsen sind, mufs auch die Ver- 


dringungsdicke eine zur Wandellipse konfokale Ellipse bilden, deren groBe und kleine Halb- 


_achse wir mit a,, b, bezeichnen (Wandellipse a, b). Es muf dann offenbar gelten 


AQ =x (a,b, — ab). 
14* 
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Bei Beachtung von b, = Va == OP ze Va? — @ +b? und a = ¢ Wp) %, b =e Sin yj erhalt man | 
schlieBlich 


ul (een ey ee (Bh aa en (21) 
a ( ‘ a! (Fj)? — 4 : 


a a 


Im Grenzfall b/a = 1 erhalten wir wieder das gleiche Ergebnis wie beim Kreiszylinder. Wie | 
beim Kreiszylinder ist auch hier F/y = 2 eine Grenze, an und unterhalb der die Grenzschicht- | 
dicke keinen endlichen reellen Wert hat. 

d) Strémung lings zweier sich schneidender Ebenen. Wir be- | 
viitzen hier die Abbildungsfunktion 


‘ i Z\nl7 L nly. 
Vs Viel C=const y9+1¢ = — tr =— le CORY (22) | 
mit | 
y= a ae sin oe (22a) | 
LTV og 82 | 
= ( 7) (22h) 
und dem Mafstabsfaktor 
1—al/y 
H =H) (23) | 
Hierbei ist y der Winkel zwischen beiden Ebenen 
und / ein noch unbestimmtes LangenmaB. In Abb. 4 
ist der Fall y = 7/2 dargestellt. Hierbeiist 
= ace sin2 3 (24) 


und 


Abb. 4. Koordinatensystem fiir die Grenzschicht- 
strémung langs zweier sich rechtwinklig schneidender (E 
Ebenen bei linear zunehmender Absaugegeschwindigkeit. 15 . (25) 


Die Absaugegeschwindigkeit an der Wand ist verteilt wie v» = — r F/I?, so daB dvg/dr = F/I? ist. 
- Damit wird die u-Geschwindigkeitsverteilung 
BL — 5 FS (FZ) sito 
u y 2 dr v \l 
Gale’ =l—e : (orm 


Die Linien konstanter u-Geschwindigkeit sind hier gleichseitige Hyperbeln und wir kénnen 
eine Verdrangungshyperbel 7, definieren, indem wir zunachst die Verdrangungsdicke in einem 
endlichen ¢ Bereich berechnen und dann die Integrationsgrenze C)—> oo gehen lassen. Mit 
F/y = — a ist 7, also definiert durch 


Une fo co fy 
lim | H® dy dt = lim [{# e—* dy dl , 
SP fe pois 
. (27) 
eal See aay ates oe | 
lim {= |, Ue Gin 2° —Ve2 + y? +. Co|t = lim — | e447 Wr Sin=2 dy. 
(a) 2 1 s Co > 2 U7] 
0 
Der Grenziibergang fiihrt auf 
1 
Ne eee (28) 


Der Abstand r, des Scheitelpunktes der Verdrangungshyperbel vom Eckpunkt ist daher 


i oae SST eee | 
Ager 12 1 =I 2% =Vaen- (28a) 


e) Gekerbte oder wellige Wand. Wir gehen aus von der Abbildungsfunktion 


. «. Je 
y+io =Insin -, (29) | 
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die einer unendlichen Folge von Senkenlinien entspricht, wobei 21 der Abstand der einzelnen 
Senkenlinien voneinander ist. Durch Zerlegen in Real- und Imaginarteil findet man 


eT or 7 é a sin 8 O0j> ¥ 
Te ah tN (Co} y — cos x), € = arc cos - (30) 


5 Gof y —cos x) 
wenny =ya/l und *% = xz/list. Der Ma®stabsfaktor ist 


__ of ¥y — cos ¥ 

~ Cof y —cos x* (31) 
Abb. 5 zeigt das so gewonnene Koordinatennetz. Wahlt man als Wandkontur 7 = 7) = konst., 
so entspricht 


No < 9,346 einer unendlichen Folge von langsangestrémten ovalen Zylindern mit Ab- 
saugung, 

% = 9,346 einer gekerbten Wand, 

My > 9,346 einer welligen Wand. 


Die Absaugegeschwindigkeit ist jeweils wie vy = — F/H verteilt. Die u-Geschwindigkeitsver- 
teilung lautet 
= (ean) es Fin YX Coty — 008 
ies eye eo Me a ee ae i 
U (32) 


C=const 


Abb. 5. Senkenreihe zur Darstellung einer gekerbten oder welligen Wand. 


4, Absaugegrenzschichten konstanter Wirbelstarke. Auch bei Annahme von wm = w,)=konst. 
ist die Gleichung (6) ohne weiteres erfiillt und man erkennt, daB wirbelfreie Lésungen und solche 
mit konstanter Wirbelstarke fiir die v, w- Komponenten linear tiberlagert werden kénnen. Wir 
erhalten somit nach (1) und (5) das Gleichungssystem 


O(v A) Ow H) _ 

LH) Ae!) <0, (33) 
O(wH) av H) 

pa 2 w, H?. (34) 


Wir beschranken uns im Folgenden auf die Falle des Kreiszylinders und des elliptischen 
Zylinders. 

a) Kreiszylinder. Im Fall einer zylindrischen Begrenzung wahlen wir wieder 
y =Inr, H =r = e", und Lisungen der Gleichungen (33) und (34) lauten 
OE, == leanne, =e 1m, 
w HH = w,e?" + konst. , | 
(35) 


also 


1 
WwW = WoT + — konst. 
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Die Gleichung (2) fiir die u- Komponente 14Bt sich lésen durch den Ansatz 


Fn 


ule, 


der auch die Randbedingungen erfillt. Ein Vergleich mit (11) zeigt, daB die u-Verteilung | 


genau gleich wie im wirbelfreien Fall ist und nur von der Absaugemenge abhangt. 


b) Elliptischer Zylinder. Wahlen wir wieder die durch (17a) und (17b) ge- 


gebenen elliptischen Koordinaten, so ist 
Js hy = £ (Co 24 —cos2€). 
Lésungen der Gleichungen (33) und (34) sind in diesem Fall gegeben durch 
vH=C, + asin 2, wH = Cy +o Sin 24. (37) 
Fir groBe Werte von 7 geht dies wieder in eine starre Drehung iiber wie im vorhergehenden 
Fall. Im Gegensatz zum Kreiszylinder ist aber die Verteilung der Absaugegeschwindigkeit 


von der Drehung abhangig. Dementsprechend ist auch die u-Verteilung von der Drehung 
abhangig und als Lésung der Differentialgleichung 


2 ; 0 2 : Ou. eu Ou 
(G. +5 eysin 22) F i (Ce }- Fm Gin 2m | Fe = (Fat aa) (38) | 


gegeben. Ks ist zweckmafig, den Lésungsansatz 
j=l + e—%1 Ai(C) + e— (4 +2)0 A,(Z) fe. (39) 


zu wahlen, durch den man Differentialgleichungen 1. Ordnung fiir die A;(¢) erhalt. 


5. Absaugegrenzschichten mit nur von 7 abhingiger Wirbelstirke. Im Fall w = w(y) | 


geht (6) tiber in 
CH Wy = 9.0), (40) 


woraus v H = F(y) folgt. Wenn wir uns auf den Fall einer kreisférmigen Begrenzung be- 
schranken und H =r = e” setzen, muB offenbar auch wH eine Funktion von 7 allein sein, 
und die Kontinuitatsgleichung wird nur erfillt, wenn v H = F = konst. ist. Dann findet man 
aber als Lésung von (40) 


Ba. 
Hi Dye (41) 

und man erhalt 

w HH =—2%0_ (Fj + 2) + konst (42) — 

F/y + 2 6 
was auch geschrieben werden kann 
Tete 
Fe 24 Do r v ait konst. 
es . 
v 


Wenn nur — F/y >1 ist, so klingt w fiir groBe r immer ab und man kann mit der Lésung (42) 
sowohl w = 0 als auch w = konst. fiir r = a befriedigen. Da wir vorher bereits unter gleichen 
Bedingungen an der Wand und im Unendlichen die Lésung w = 0 im ganzen Striémungsfeld 
erhalten hatten, bestimmen offenbar in einem ringférmigen Strémungsfeld die Randbe- 
dingungen nicht eindeutig die Lésung. Vielmehr kénnen Wirbelstrémungen iiberlagert werden. 
Fir die u-Komponente erhalt man wieder die durch (11) gegebene Verteilung. 


6. Absaugegrenzschichten mit nur yon ¢ abhiangiger Wirbelstirke, Im Fall @ = w(C) geht 
(6) tiber in 


wHao' =v", (43) 
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wobei die Striche jetzt Ableitungen nach ¢ bezeichnen. Offenbar mu8 jetzt w H = G(¢) sein, 
so daf aus der Kontinuitatsgleichung folgt 


vH=n7G(0)+ Fo). 
Durch Einsetzen in (5) ergibt sich 


4G’ 4+ F’ = — 2 H%a(l). (44) 
Offenbar ist nur im Fall rechtwinkliger Koordinaten, also fir H =1 ein sinnvolles Ergebnis zu 
erwarten. Um (44) zu erfiillen, mu8 G’’= 0 sein. Da w nur eine Funktion von allein ist, an 
der Wand (7 = konst.) aber verschwindet, muS w im ganzen Strémungsfeld gleich Null sein. 
Wegen (43) ist dann w’’='0, und aus (44) folgt F’’= 0, also F=a+bC + cl2. Dies ent- 
spricht aber einer Poiseuille-Strémung. Durch Koordinatenverschiebung kann das in z lineare 
Glied zum Verschwinden gebracht werden. Wir kénnen also b = 0 voraussetzen. Fiir das 
Druckfeld erhalt man nach (3) und (4) 

1 op Op aan 


3) en Se Ph > 8 ==) > 
also 
p =konst. + 2ucy. 
Zur Lésung der Gleichung (2) fiir die u- Komponente kann man den Ansatz 


u 


Foooae L—e—*" [((Fo(m) 4-0? Fe(n) C4 Fy (ny) +s +] (45) 


machen. Doch ist die Berechnung der Funktionen F; und die Befriedigung der Randbedin- 
gungen an der Wand mihsam. 


7. Spezielle Losungen, bei denen die Querbewegung die Strombahnen einer Potentialstromung 
besitzt. Wir suchen nach Lésungen, fiir die zwar w = 0, aber v H = F(C) + konst. ist. Die 
Querstrémung soll also zwar den Bahnen einer Potentialstrémung folgen, ohne aber selbst 
eine solche zu sein. Die Méglichkeit derartiger Strémungen ist bereits von Hamel} mit dem 
Ergebnis untersucht worden, daB nur Strémungen auf logarithmischen Spiralen (mit den 

_ Grenzfallen der reinen Radial- und der reinen Kreisstrémung) diese Bedingung erfiillen. 
Die umgeformten Gleichungen (3) und (4) 


Ft (re) arene: & (im): (46) 


oo hae eae oa at \ H? 

alm 
1 im) o) ee Pee (47) 
‘oa goes . hie. én \ H? 


kénnen zusammengefabt werden in 


i 
Fe OP EN, OF F ne 
oy =" al) + ae (e) . 
mit H = r/¥1 + y? und Inr = (7 + yC)/(1 + y?) vgl. (15) und (14) folgt dann aber 

—2FF’ =7[4F'+F" (14+?) —4yF’]. (49) 


Wenn wir uns weiterhin auf die Radialstrémung (y = 0) beschranken, ist die Lésung nach 
Hamel durch elliptische Funktionen gegeben. 
Die Differentialgleichung (2) fiir u 14Bt sich umschreiben in 


2 
Behr ges eget ps) (50) 


Mit dem Lésungsansatz 


un, = rei g(C) 
1G, Hamel, Jahresber. Dtsch. Math, Vereinig. 25 (1916) S. 34. 


206 Wuest: Asymptotische Absaugegrenzschichten usw. Ingenieur-Archiv 


erhalt man fir f,, die Differentialgleichung 


aay sree) (51) 


Da F/y eine in€ periodische Funktion ist, handelt es sich offenbar um eine Hillsche Differential- 
gleichung. Fir n = — 2 lautet eine spezielle Lésung f, = F’. Wir betrachten folgende beiden | 


Aufgaben. 


a) Kreiszylinder mit ungleichfirmiger (aber nicht willkir- 


liehe'r) Verterlhua'g der Absaugegeschwindigkeit. Wir kénnen hier fir | 


die Lésung den Reihenansatz 


eo hee sy pg iat 


U n=—l 


machen. Da die f, selbst wieder periodische Funktionen sind, kénnen sie in eine Fourierreihe 
entwickelt werden, und die c, kénnen sicherlich so gewahlt werden, daB fir r = a (Zylinder- | 
oberflache) alle periodischen Glieder sich gegenseitig aufheben und die konstanten Glieder die | 


Summe 1 ergeben. Das Problem ist also offenbar lésbar. 


b) Strémung im Winkelraum zweier Ebenen mit Absaugung| 


in der Schnittgeraden. Fir die radiale Senkenstrémung zwischen ebenen Wanden 
sind von Hamel Lésungen entwickelt worden, welche die Randbedingungen fiir die Querstré- 


mung erfiillen. Man erkennt jedoch leicht, daB es nicht moglich ist, hierzu eine u-Strémung zu | 
finden, die fiir geniigend groBen Abstand von den Wanden in den konstanten Wert u = U | 


ibergeht und die Haftbedingung an der Wand erfillt. 


8. Nachbarlésungen zur wirbelfreien Querbewegung. Im Folgenden sollen Nachbarlésungen 
zuvH = F, = konst.; wH = 0 untersucht werden. Dabei interessieren vor allem geschlossene 
Kérper, bei denen die Abweichungen A(v H) von der wirbelfreien Verteilung in eine Fourier- 
reihe zerlegt werden kénnen. Zunachst liefert schon die Annahme w = 0 periodische Lésungen 
der Form 


A(v H), = X Cre" —*9), A(w H), = —i S C,e—"— #9), (52) 


mit denen aber die Randbedingung w = 0 an der Wand nicht erfillt werden kann. Unter der | 


Annahme, da die Abweichungen von der Grundlésung v H = F, nur gering sind, gewinnen 
wir naherungsweise weitere Lésungen dadurch, daB wir in der Differentialgleichung (6) quadra- 
tische Glieder der Abweichungen vernachlassigen und somit die linearisierte Differential- 
gleichung 


Fy dw ao) ow 
> Oy Ot OF (64 


gewinnen, deren Liésungen wir in der Form 


Des Sree 


n 


F Foe ys 
Ae ee ay a ee ne 
Aus (1) und (5) kann man dann vH und wH berechnen. Diese Lésung soll auf den Fall des 
Kreiszylinders angewandt werden. 


Wegen H? = e?” findet man aus (54), (1) und (5) 


AGH), = re ee | 


ansetzen mit 


+ 2)? — n2 


(—An + 2)? — n? 


n 


55 
sintt= SBP ease reine, | a | 


(54) 


| 
| 
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Aus der Zusammenfassung von (52) und (55) bei Beachtung der Bedingung w = 0 fiir 7 = 1% 
findet man 


aes poe are eee re ein’, | 
| (56) 


A(wH) = > —iC, leo ght + (—In + Din oi sek 


n 


Die C, bestimmen sich aus der vorgegebenen Verteilung von v an der Wand!. Fir/, = 2 +n 
oder F,/y= (4+ 4n)/(2 — n) verschwindet iibrigens das betreffende Teilglied von w H fir alle 7. 
Eine nur durch dieses Teilglied hervorgerufene Stérung klingt nach aufen ab wie 


A(v H) = Cye—" im) tine (57) 


oder mitr =iny 
Tr nm 
Alery = ¢; (") cosn€. 
rt. 
Wir zerlegen auch die u-Verteilung in eine Grundverteilung u, und eine Stérungsverteilung 


Au. Mit dem gleichen Grad der Anniaherung gilt dann 


Au _ Au Au 
oan on” F oc? 


— A(vH) Ss. (58) 


Die homogene Gleichung (d. h. ohne das letzte Glied) hat die gleiche Form wie (53) und also 
auch entsprechende Lisungen. Fiir die inhomogene Gleichung machen wir den Ansatz 


F, 
eee —H n a Reng He cali ear rene) ai 
Au; = De mah pea gear a7? q ( yin Bliee ° (59) 
Die Konstanten D,, und E,, bestimmen wir durch Einsetzen dieses Ansatzes in die Differential- 
gleichung, wobei wir A(v H) aus (56) entnehmen und 0u,/07 = — (Fy/) e» —" beriicksichtigen. 
Wir finden 
Cn G 
pee ‘Zs 2 
Fyn LRG ES Fie tee on el 


Insgesamt erhalt man also als Lésung fiir die Stérgeschwindigkeiten Au 


Cn 


n2 


ret a we 
ton (—dn + 2) {(— dn + 2) + ZL (— dn + 28 — vf] 


0 


et te + (—In + 2) (1 — my) 


Aig 


Fo —) tin’ 
x £0” 7 lo) Be Reg wie ne 
wobei die Konstanten K,, aus dem Verschwinden fiir 7 = 7) bestimmt werden, wahrend die C, 
wie oben durch die Verteilung der Absaugegeschwindigkeit an der Wand (vgl. Formel (39)) 
gegeben sind. 


Bei einer speziellen Stérung, welche der Bedingung (40) geniigt, ist 
Atrc=Wv. 


9. Zusammenfassung. Die Untersuchung bezieht sich auf die Strémung langs einer unend- 
lich langen zylindrischen Flache, an deren Oberflache die Grenzschichten laufend abgesaugt 
werden. Die Geschwindigkeitskomponenten sollen nicht von der Koordinate in Achsrichtung 
abhangen und in groBem Abstand von der Zylinderflache in eine gleichmafhige Parallelstrémung 
in Richtung der Zylinderachse iibergehen. Aus den allgemeinen Differentialgleichungen dieses 
Problems lassen sich verschiedene Sonderfalle heraussondern. Bei Verschwinden der Wirbel- 
komponente in Hauptstrémungsrichtung ergeben sich einfache Lésungen, die als Verall- 
gemeinerung der homogenen Absaugegrenzschicht an einer ebenen Wand aufgefaBt werden 
kénnen. Als Anwendungen werden u. a. die Strémung lings eines ruhenden und rotierenden 
Kreiszylinders, eines elliptischen Zylinders sowie im Winkelraum zweier sich schneidender 


1 Dieser Sonderfallist auf eine etwas andere Weise bereits von G. Hamel a. a, O. untersucht worden, 
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Ebenen behandelt. Nachbarlésungen hierzu werden durch Linearisierung der Differential- 
gleichungen gewonnen. Weitere weniger einfache Sonderfalle ergeben sich dadurch, da man 
die Wirbelkomponente in Hauptstrémungsrichtung als konstant oder nur von einer Verander- 
lichen abhangig voraussetzt. SchlieBlich sind noch spezielle Lésungen méglich, bei denen die 
Querstrémung eine Hamelsche spiralférmige Strémung ist. 

Die untersuchten Faille lassen folgendes vermuten: 

1. Bei einem im Endlichen geschlossenen zylindrischen Kérper (z. B. elliptischer Zylinder) 
und endlicher Absaugemenge (je Langeneinheit des Zylinders) existieren asymptotische Lésungen 
fiir die u-Komponente der Strémungsgeschwindigkeit in Richtung der Zylinderachse, welche 
die Hafthedingung an der Wand erfiillen und in geniigenden Abstand von der Wand in eine 
konstante Parallelstrémung (u = U) iibergehen. Die mit diesen Randbedingungen gewonnenen 
Lésungen sind jedoch nicht eindeutig, weil noch Wirbelstrémungen iiberlagert werden kénnen. 

2. Bei einem offenen zylindrischen Kérper, dessen Flachen sich ins Unendliche erstrecken 
(z. B. zwei sich schneidende Ebenen), existieren asymptotische Lésungen mit den vorher- 
genannten Eigenschaften nur dann, wenn nicht unendlich groBe Teilabschnitte der Oberflache 
je Langeneinheit des Zylinders frei von Absaugung sind. Dabei wird vorausgesetzt, da die 
Grenzschicht nicht in anderer Weise Energie zugefithrt erhalt (z. B. durch Druckgefalle in 
Richtung der Zylinderachse oder durch eine bewegte Wand). 


(Eingegangen am 2. August 1954.) 
Anschrift des Verfassers: Dr. W. Wuest, Gottingen, Bunsenstr, 10. 


, 
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Uber den Auftrieb eines Blasfliigels. 
Von H. B. Helmbold. 


1. Zirkulationsauftrieb und Impulsauftrieb. Wird ein Strahl aus einem Blasschlitz tangen- 
tial zur Fliigeloberflache nach hinten geblasen, so entsteht iiber den ohne Strahl bereits vor- 
handenen Auftrieb hinaus ein Zusatzauftrieb. Der Gesamtauftrieb besteht dabei im all- 
gemeinen Fall aus zwei verschiedenartigen Anteilen, dem Zirkulationsauftrieb A,, der fir das 
Auftreten von freien Lingswirbeln in der Strémung hinter dem Fligel und einen entsprechen- 
den induzierten Widerstand verantwortlich ist, und dem Impulsauftrieb 4,, der die Reaktion 
fir die sehr weit hinter dem Fliigel tibrigbleibende Abwartskomponente des Impulsstroms im 
Strahl selbst darstellt und dessen Erzeugung mit einem negativen Widerstand, dem Blas- 
schub, verbunden ist. Eine Abschatzung dieser Anteile auf Grund von Impulsbetrachtungen 1 
zeigt, daB der Gesamtauftrieb reiner Zirkulationsauftrieb ist und der Impulsauftrieb verschwin- 
det, wenn die Spannweite gegen Unendlich geht, womit man zur ebenen Strémung gelangt; 
dabei wird der Blasstrahl, der den Fliigel unter einem Winkel 6y gegen die Anblasrichtung 
verlassen hat, von der umgebenden Strémung schlieBlich vollstandig in die Richtung der un- 
gestérten Strémung umgelenkt. 


2. Der Blasfliigelauftrieb in ebener Stroémung. Nun lehrt die Erfahrung, da8 der Zusatz- 
auftrieb 4A eines Blasfliigels in ebener Strémung ein Vielfaches N der Abwartskomponente 


AJ, =obh[(V + v)? — V?] sin dg (1) 
des Impulsiiberschusses im Strahl an der Hinterkante betragt (9 = Luftdichte, b = Spann- 
weite, h = Dicke des Blasstrahls, V = Geschwindigkeit dicht auBerhalb des Strahls, 
v = Ubergeschwindigkeit im Strahl). Beispielsweise ergibt die Auswertung einer von 
P. Rebuffet 2 veréffentlichten Versuchsreihe * fiir den Vergleichsfaktor N einen Wert nahe 
an 5. Durch die Umlenkung des Strahls in die Anblasrichtung wird also bedeutend mehr 
Zusatzauftrieb erzeugt als dem urspriinglich erzeugten Abwartsimpuls 4J, unmittelbar ent- 
spricht. Um das erklaéren zu kénnen, fragt man am besten zuerst einmal nach den KrAaften, 
die den Strahl umlenken und damit den urspriinglichen Abwartsimpuls 4J, abfangen. Das 


kénnen nur zusatzliche Druckkrafte ASS 


sein, die die umgebende Strémung 


(,,Grundstrémung**) auf die Ober- Ko Blastutt SS : 
flache des Blasstrahls ausibt, Geblise lll) a 


Druckkrafte von der gleichen Art 
wie der Auftrieb am Tragfliigel. Es 
muB also, wenn wir den Blasstrahl 
als Fremdkérper betrachten, in der Grundstrémung iiber den Blasstrahl hinweg noch ein 
Drucksprung vorhanden sein, der wegen der Konstanz des Gesamtdrucks in der Grund- 
strémung von einem entsprechenden Geschwindigkeitssprung begleitet ist: 


V4 (702) 


Abb. 1, Ebene Platte mit Blasstrahl aus der Hinterkante (schematisch). 


Fn ae I 
TE We Wea = (V7 Va) (2) 
oder, wenn wir mit V,, als ungestérter Geschwindigkeit und mit AV als Ubergeschwindigkeit 
At =VL+AV umd V,—V,=6 (3) 
schreiben, Peo—Pi =0(V, +AV)e. (4) 


Die Bezeichnungen der Gleichung (2) sind durch Abb. | erklart, die den einfachsten Fall, die 


_ ebene Platte mit Blasstrahl aus der Hinterkante, illustriert. Den Geschwindigkeitssprung ¢ 


1 H, B. Helmbold, On the Lift of a Blowing, University of Wichita, Engineering Study No. 135, 

April 1954. 
es P. Rebuffet, Aérodynamique expérimentale, S. 461, 2. Aufl., Paris et Liége 1950. 

3 Bei diesen Versuchen wurde die Blasluft einem gegeniber dem Fliigel ruhenden Luftraum ent- 
nommen. Bei den folgenden Uberlegungen ist angenommen, daf die Blasluft der Strémung entnommen 
wird und mit der Geschwindigkeit V in den Fliigel eintritt, um ihn an der Hinterkante mit der gleich- 
gerichteten Geschwindigkeit V + v zu verlassen (Abb. 1). 
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darf man als Dichte einer Belegung des Strahls mit tragenden Wirbeln betrachten. Der Im- 
pulssatz, angewandt auf eine den Strahl dicht umschlieBende und an der Hinterkante sowie 
im Unendlichen senkrecht schneidende Kontrollflache, besagt, daB wegen Fehlens auferer 
Krafte die Aufwartskomponente der Druckkrafte mit der Anderung der Abwartskomponente 
des Impulsiiberschusses im Strahl im Gleichgewicht steht: 


0= | e(Va + AVe)eds— 4. (5) 
L 


Hier ist s der Abstand von der Fligelvorderkante lings der Wirbelflache gemessen. L die 
Fliigeltiefe und AV, die x-Komponente der Stérgeschwindigkeit [(V, + 4V)dx =(V, 
-- AV,) ds]. Dabei ist schon beriicksichtigt, daB ohne Strahl (v = 0) mit dem ImpulsiiberschuB 
AJ =obh{[(V + v)? — V2] auch die Wirbeldichte ¢ verschwinden muf [siehe auch Gl. (14)]. 
Die Abwartskomponente ist AJ, = AJ sin dy. Der Auftrieb des Tragfligels mit Blasstrahl 
ergibt sich aus einer entsprechenden Anwendung des Impulssatzes auf eine den Fligel dicht 
umschlieBende und den Strahl an der Hinterkante senkrecht schneidende Kontrollflache: 


B 
A ANAS 
A = | (Ve + AV«) (vy +47) ds +. (6) 
i 
Hier haben wir die Dichte y der Belegung des Fliigels mit tragenden Wirbeln in den ohne 
Strahl bereits vorhandenen Grundanteil yy und den vom Strahl verursachten Zusatzanteil Ay, 
dessen Zustandekommen wir noch zu diskutieren haben, zerlegt. Durch Addition der beiden 


letzten Gleichungen ergibt sich 


L 


A ie} 
gen Ue: AV:) (Yo + Ay) ds + [e(V.+ AV,) ¢ ds. (7) 
L 


fi 
Der Auftrieb des Tragfliigels mit Blasstrahl kann aber auch durch Anwendung des Impuls- 
satzes auf eine Fliigel und Strahl einschlieBende Kontroliflache festgestellt werden, die iberall 
in sehr groBer Entfernung vom Fliigel verlauft und den Strahl im Unendlichen senkrecht 
schneidet (Satz von Kutta-Joukewski), 


L (oe) 
as ea yf 
: eY.|f 6 + Ay)d eS (8) 


Da die Ausdriicke (7) und (8) iitbereinstimmen miissen, ergibt sich die Beziehung 
L 2 
[ AV. (v9 + dy) ds + [ AVe-eds =0, (9) 
0 L 


die sich unabhangig von ihrer Herleitung durch Anwendung des Induktionsgesetzes verifi- 
zieren laBt (Gegenseitigkeit der Wirbelinduktion). 
Subtrahieren wir nun vom Auftrieb des Fliigels mit Blasstrahl nach (8) den Auftrieb des 
Fliigels ohne Blasstrahl 
L 


A 
B= 0Vm [ rods, (10) 
0 


so erhalten wir den vom Strahl verursachten Zusatzauftrieb 


L ce) 
AA 
<4 = Vo([ dy ds + | eds). (11) 
L 


0 


Dieser Zusatzauftrich ist nach (5) um folgenden Betrag gréBer als die Abwartskomponente 
des Impulsiiberschusses im Strahl an der Hinterkante: 


BA AS, r re 
AS =0v, | dy-ds —0[ AV. eds. (12) 
0 L 
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Das laBt sich mit Hilfe von (9) umformen in 


Ze ; 

50 | V+ AV.) dy + AVe: has. (13) 

0 

Der Zusatzauftrieb ist also mindestens dann gréfer als die Abwirtskomponente des Impuls- 
iberschusses an der Hinterkante, wenn die vom Blasstrahl verursachten Zusatzgréfen A V,, 
und Ay langs der gesamten Fliigeltiefe positiv sind. Das ist aber bei positivem Ausblaswinke] 
Ox des Strahls immer der Fall. Denn da nach (15) die Wirbeldichte ¢ hinter dem Fliigel itberall 
positiv ist und diese Wirbelbelegung mit Bezug auf die Anblasrichtung tiefer liegt als die Hinter- 
kante, induziert sie lings der gesamten Fligeltiefe eine Stérgeschwindigkeit mit positiver 
x-Komponente AV, und mit positiver y-Komponente (Aufwind); der induzierte Aufwind 
aber bedingt eine itberall positive Zusatzbelegung Ay. 

Anschaulich gesprochen steht der Fliigelin einer durch den Strahl gestérten Strémung. Die 
Stérung besteht darin, daB der Uberschu® des schnellen gekriimmten Strahls an Zentrifugal- 
kraften gegeniiber der Grundstrémung den Strahl nach unten und damit die Grundstrémung 
aus ihrem urspriinglichen Verlauf heraus drangt (Abb. 2). Der Blasstrahl, den man zu einem 
diinnen Blech erstarren lassen darf, wenn man die beim Erstarren fortfallenden Zentrifugal- 
krafte durch auBere Krafte ersetzt, wirkt dabei wie eine mit dem Fliigel nicht verbundene 
Fortsetzung desselben, so daB sich die 
tragende Wirbelflache iiber die Hinter- AA 
kante hinaus als nunmehr ebenfalls 2=== 
tragende Wirbelflache fortsetzt und —. 
die wirksame Fligeltiefe vergréBert mit Strahl 
wird, wodurch dann auch die Um- Abb.2. Spaltungs- und ZusammenfluB-Stromlinien (schematisch). 


_strémung der Vorderkante verstarkt 
und die Spaltungsstromlinie in Fligelnahe nach hinten unten verlagert wird. Das bedeutet er- 


’ 


: 
; 


héhte Zirkulation um den Fliigel selbst, wo die Stérung sich in einem induzierten Aufwind, einer 
induzierten Ubergeschwindigkeit und einer endlichen Wirbeldichte an der Hinterkante ausdrickt. 


3. Die Integralgleichungen der Wirbelverteilung. Um dem Uberschuf des schnellen ge- 
kriimmten Strahls an Zentrifugalkraften das Gleichgewicht zu halten, bedarf es eines Druck- 
sprungs von der GriéBe 

h h 

Po — Pi =e ((V ++)? — V*]— = 2 Ap, —, (14) 
worin V = V, + AV die értliche Geschwindigkeit der Grundstrémung, v die Ubergeschwin- 
digkeit des Strahls gegeniiber V, r den 6rtlichen Kriimmungsradius und Ap, den UberschuB 
des Strahls an Gesamtdruck bezeichnet (vgl. die eingangs zitierte Arbeit des Verfassers). 

Durch Vergleich mit (4) erhalt man die Wirbeldichte hinter dem Fligel 

5 Ape By As, as 

fan pe me re (15) 
Sind die vorkemmenden Neigungen und Kriimmungen hinreichend klein, so daB AV < V,, 
und h = const angenommen werden darf, dann ergibt das Gesetz der Wirbelinduktion zwei 
Integralgleichungen zur Bestimmung der Funktionen y(x) und ¢(x); ihre linken Seiten lauten 


aS Es _ im Bereich 0X x<L, (16) 

x x 

= lisa jee im Bereich x >L, a0) 
x 


wenn yp und yz die Ordinaten der tragenden Wirbelfliche am Fliigel und am Strahl bedeuten. 


| ‘Die rechten Seiten lauten in beiden Bereichen formal iubereinstimmend 


L oe) 
dy(«x) Bey 1 y( x’) deat 1 Apg h d*ys(x’) dx’ (18) 
dy Fave a x 


Dak Ge dx’* x’ —x 


mit q,. = @ V2,/2. Die Funktion yp(x) ist durch die Gestalt des Fligelprofils vorgegeben; 
die unbekannten Funktionen ys(x) und ¢(x) sind durch Gl. (15) miteinander verknipft. 
Lésungen dieses Gleichungssystems stehen zur Zeit noch nicht zur Verfiigung. 


(Eingegangen am 2. August 1954.) 
Anschrift des Verfassers: Professor H. B. Helmbold, 646N. Bluff Avenue, Wichita 6, Kansas USA. 
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Wachstum der laminaren Grenzschicht an schrag angestrémten 
Zylindern bei Anfahrt aus der Ruhe*. 


Von H. Wundt. 


1. Einleitung. Die dreidim ensionalen Grenzschichtstrémungen inkompressibler Medien sind 
der rechnerischen Behandlung wesentlich schwieriger zuganglich als die ebenen, weil die Me- 
thode der Stromfunktion zur Integration der Kontinuitatsgleichung fiir drei Ortskoordinaten 
im allgemeinen nicht anwendbar ist. Es ist daher von Bedeutung, daB in Sonderfallen das 
System der Grenzschicht-Differentialgleichungen ,,zerfallt“‘, wie zuerst Prandtl+ bemerkt hat. 
Unter diesem als ,,Unabhangigkeitsprinzip“ bezeichneten Sachverhalt versteht man, da® nicht 
in jeder der Gleichungen zugleich samtliche Geschwindigkeitskomponenten der Grenzschicht- 
Strémung auftreten, weshalb man erst ein Teilsystem fiir zwei Komponenten lésen kann. Die 
dritte Komponente wird nachtraglich aus der letzten Gleichung bestimmt und beeinfluBt die 
anderen Komponenten nicht mehr. 

Ein solcher Fall liegt beim schrag angestrémten Zylinder vor, auf den sich die folgenden Unter- 
suchungen beziehen. 


2. Aufstellung der Grenzschichtgleichungen. Wir verstehen, wie neuerdings allgemein tiblich, 
unter (x, y, z) ein krummliniges Koordinatensystem derart, daB x die krummlinige Wandbogen- 
lange langs der Kontur eines ebenen Schnitts senkrecht zu den Erzeugenden des umstrémten 
Zylinders ist, y der senkrechte Abstand von einem Bezugs-Normalschnitt, und z der nach auben 
positive Abstand von der Zylinderoberflache. Die hierauf bezogenen Geschwindigkeitskomponen- 
ten der Grenzschichtstrémung heiBen u, v, w. p sei der Druck, t die Zeit, 9 die konstante Dichte 
und y die konstante kinematische Zahigkeit des Mediums. 

Mit diesen Bezeichnungen lauten die Gleichungen fiir eine allgemeine dreidimensionale, in- 
kompressible, instationire Grenzschichtstr6mung (vgl. etwa Schlichting?) : 


Uy + Wily + vy + wu, =—— Py +9 a (1a) 
1 (Bewegungsgleichungen) 

Sp crest Ca Oy ht Ua tee tate else (1b) 

uz +o, +w,=0. (Kontinuitatsgleichung). (1c) 


Die Indizes bedeuten partielle Ableitungen nach den betreffenden Variablen. 
Hierzu treten die Randbedingungen 


furz=- 0: i (X. 920) — 0.) ; 

ah oe (Haftbedingungen) 

w(x, ¥,t)=0, (Undurchlassigkeit der Wand) (2) 
fiir z /Re —>o: u(x, y, t)h—> U(x, y, t), 


v(x, y, th» V(x, y, t) . 


U° und V° sind die Werte der mit grofen Buchstaben U, V, W bezeichneten AuBenstrémung an 
der AuBenflache der Grenzschicht; Re ist die mit einer charakteristischen Lange und einer charak- 
teristischen Geschwindigkeit gebildete Reynoldssche Zahl. 

Die Druckglieder sind aus den Bernoullischen Gleichungen 


U? + UU + VU; =—— px, (3a) 
Vo + UW VE+ Ve = —— py (3b) 


* Die Arbeit entstand am Mathematischen Institut der Universitat Freiburg i. Br., Abteilung Angewandte 
Mathematik. Sie wurde durch das ,,Air Research and Development Command, United States Air Force“ 
finanziell gefordert. 

 L. Prandtl, Uber Reibungsschichten bei dreidimensionalen Strémungen, Betz-Festschrift (1945), S. 134. 

® H. Schlichting, Grenzschichttheorie, Karlsruhe 1951 (vgl. bes. Kap. X und XI). 
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zu bestimmen. Bei gegebener Kérperform kann man die aus den Eulerschen Bewegungsglei- 
chungen der idealen Fliissigkeit zu berechnenden Geschwindigkeitskomponenten U® und V° und 
damit den Druckgradienten langs der Kérperoberflache als bekannt ansehen. 

Von dem System (1) mit (2) und (3) ist, abgesehen von einer einzigen Arbeit!, noch keine 
Lésung bekannt, in der alle Terme miteinander in Wechselwirkung stehen. Wir suchen daher 
solche Falle, in denen u und w aus (1a) und (1c) unabhangig von v bestimmbar sind, wahrend v 
nachtraglich aus (1b) berechnet werden soll. Dazu miissen die v enthaltenden Terme in (1a) und 
(1c) verschwinden, d. h. v, = 0 sein und wegen v + 0 auch uy = 0. Die Tangentialkomponenten 
der Geschwindigkeit in der Grenzschicht sollen von der einen Wandkoordinate y nicht ab- 
hangen. Infolgedessen miissen wir die Randbedingungen und die Anstrémverhiltnisse y-unab- 
hangig wahlen, was eine von y unabhingige AuSenstr6mung U, V, W (als Funktionen von x, z, t) 
erfordert. 

Die genannten Bedingungen sind erfiillt fiir einen schrag angestrimten, unendlich langen 
Zylinder in rotationsfreier AuBenstrémung, wobei die nachtraglich zu berechnende Geschwin- 
digkeitskomponente diejenige parallel zu den Erzeugenden ist. Wird der Zylinder in groBem 
Wandabstand unter einem festen Winkel a gegen einen Normalschnitt angestrémt, der durch 
tga = V. /U,, bestimmt ist”, so ist die Zeitabhingigkeit von U und V dieselbe, und V ist nicht 
von x abhangig. Insbesondere gilt an der Wand bzw. in grenzschichttheoretischer Naherung am 
AuBenrand der Grenzschicht V° = V°(t). 

Auf den damit bezeichneten Fall beziehen sich die folgenden Untersuchungen. Die Ber- 
noullischen Gleichungen (3) reduzieren sich auf 


Up Ue Ut =—— px, (4a) 
1 
Ve ee aoe, (4b) 
Wahrend der Druckgradient y-unabhangig ist, gilt dies fiir den Druck selbst nicht; er ist linear 
in y. 

Die Grenzschichtgleichungen (1) fiir den schrag angestrémten Zylinder lauten damit 

a, +uu, + wu, = U; + UOUS + yuy,, (5a) 
v,tuv, +wr, =V; + Y Viz 9 

UU, -w, =0 (Sc) 


mit den Randbedingungen 
fiir z — 0: u(x,t) — 0 (x, 1) = w(x, 1) = 0, | (6) 
fiir z /Re—> 00: u(x, t) > U(x, t); v(x, t) > V(t). \ 

Aus der Giiltigkeit des Unabhangigkeitsprinzips folgt, daB die vor allem interessierenden 
Ablésungserscheinungen der Grenzschicht weder zeitlich noch értlich durch das Hinzutreten einer 
Querkomponente V_, beeinfluBt werden. Alle diesbeziiglichen Daten sind vielmehr bereits der 
Lésung fiir die Senkrecht-Anstrémung allein zu entnehmen?. 

3. Integration der Grenzschichtgleichungen fiir verschiedene Bewegungsgesetze mittels der 
Blasiusschen Zeitreihe. Die Integration des Systems (5) kann nun in zwei getrennten Schritten 
vorgenommen werden: 

(a) Ermittlung von u und w aus (5a) und (5c), 

(b) Ermittlung von v aus (5b) unter Beriicksichtigung von (a). 

Fir zwei Falle von Bewegungsgesetzen hat bereits Blasius 1908 in seiner Dissertation’ auf 
Anregung Prandtls Losungsansatze gegeben, die im folgenden verwendet werden, namlich fiir 

(x) ruckartige Anfahrt des Zylinders. aus der Ruhe zu gleichférmiger Geschwindigkeit, 

(B) gleichférmige Beschleunigung aus der Ruhe heraus. 

Gértler zeigte, daB sich die Blasiussche Methode auf Bewegungsgesetze der Form 


U = konst. « (n = 0, ganz) (7) 
verallgemeinern laBt* °, Die Falle (a) und (f) sind darin mit n = 0 und n = 1 enthalten. 


1 K. H. Thiriot, Z. angew. Math. Mech. 20 (1940), S. 1. ith ¢ 

2 V., und U,, sind die Komponenten der ungestorten Anstromgeschwindigkeit fern vom Zylinder. 
3 H. Gortler, Arch. Math. 3 (1952), S. 216. 4 H, Blasius, Z. Math. Phys. 56 (1908), S. 1. 

5 H. Gértler, Ing.-Arch. 14 (1944), S. 286. 6 H. Gortler, Arch. Math. 1 (1948/49), S. 138. 
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Obgleich das Problem (a) der Senkrecht-Anstrémung in den genannten Arbeiten von Gértler 
bereits behandelt ist, soll hier der Gedankengang der. Lésung doch zusammenfassend wiedergege- 
ben werden, da einerseits, wie am Beispiel des Kreiszylinders (Abschnitt 5) ersichtlich, die Glei- 
chungen erneut mit héherer Genauigkeit numerisch gelést werden muBten, zum anderen damit die 
Methodik zur Ermittlung der Querkomponente v bereits gegeben ist und daher dann kurz gefaBt 
werden kann. 

Wir nehmen eine vorlaufig noch beliebige natiirliche Zahl n fiir das Anstroémgesetz (7) an, so 
daB die Potentialstrémung am Aufenrand der Grenzschicht 


(0 ina (Do 
l q(x) fire > 0 


Ua nt) a= (8) 


ist. Durch q(x) werden die Daten des speziellen Problems (Kreiszylinder, elliptischer Zylinder 
o. dgl.) ausgedriickt, berechenbar aus den Eulerschen Gleichungen der idealen Fliissigkeit. 

Blasius integriert die Kontinuitatsgleichung (5c) durch Einfiihrung einer Stromfunktion 
w(x, 2, t) mit u = y,, w = —y,. Die weitere Rechnung beruht nun auf einer Iteration. Bei Be- 
wegungen aus der Ruhe heraus sind fiir ganz kleine Zeiten die quadratischen Tragheitsglieder 
unwesentlich; aus (5a) entsteht in nullter Naherung eine lineare inhomogene Gleichung vom 
Warmeleitungstyp. Das entstehende Defizit gegeniiber der Lésung der vollen Gleichung wird als 
inhomogenes Glied in die Gleichung fiir die erste Naherung eingegeben, und so fort. Praktisch 
lauft dies auf eine Reihenentwicklung der Stromfunktion nach Potenzen der Zeit hinaus, die, 
eingesetzt in (5a), eine rekursiv zu lésende Folge linearer Differentialgleichungen fiir die K oeffi- 
zientenfunktionen liefert. 


Zuvor fihren wir wie in der Warmeleitungsgleichung mittels der Ahnlichkeitstransforma- 


tion 1 = 2/2 yore die dimensionslose Variable 7 ein, wodurch die Anzahl] der unabhangigen Ver- 
anderlichen um eins vermindert wird. Die Grenzschichtgleichung (5a) lautet in den Koordi- 
naten x, 1, t, wobei das aus (8) berechnete Druckglied eingesetzt ist, 


n w 


1 


Ue Fd P UUs -| 2\ve u,=nge—!+qg?™+ 4t¢m: (9) 
Der Strich bei q bedeutet die Ableitung nach seiner Variablen x. 
Fir die Stromfunktion kann man nun den Ansatz machen 
(x41) = 2p vem BA ves x, i(x0). (10) 


Wir entwickeln also nach Potenzen von t"+!, d.h. nach Potenzen des vom Kérper zurick- 
gelegten Weges. Nach Einsetzen in (9) entsteht jetzt durch Vergleich gleicher t-Potenzen die 
Folge rekursiv zu lésender linearer Gleichungen fiir die 7,, ,(x, 7): 


Oy, 0 Ox, 0 Ox, 0 
One +27 Or? aan a ae 4ngq(x) , (11a) 
eyo eee OXny Xn 0 Xn 0 O%n0 Xn 0 
| se n a n n n n vA 

ars! 24) On? 4 (2n + 1) Cn 4 On Ox On Say On? od q() q (2)| 2 (11b) 
oy Oe Ox, OX, 9 OX te) 0? 0 0? to) 0? 

nd 2a? 4 (3n a) fa 4 |Fno- 2 Xnl And Xn0 Xnl Xn1 7% Xn0 , (11e) 

” nN Ox On On Ox On Ox On? Ox On? 

ay i 0 0 0? 0 0? } 2 2 

ans J x 3__4(4n + 3) Xn3 =4| Xn0 Xn2 Xnl Xn1 Xn2 tno Png Png 

On On On On OxOn On Ox On On Ox On Ox On? 


OXn1 Ont O%n 2 Ox, 0 (11d) 
Ox On Ox On? 


ates | 


Die Iteration kommt dadurch zum Ausdruck, daB die rechten, inhomogenen Seiten jeweils nur 
Terme enthalten, die aus den Lésungen schon vorher behandelter Gleichungen berechenbar sind. 

Nun wird gefordert, daB die zur nullten Naherung gehérende Beiwertfunktion xno bereits 
die wesentlichen Randbedingungen erfiillt, wahrend die folgenden 4%n,a (A > 0) am Innen- und 
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am AuBenrand der Grenzschicht keine Beitrage mehr liefern sollen, m. a. W.: 


0 
Pasa = ne a(a) » 
fiir ep == Oe On, 4 fiir "—> 00: i (12) 


Die x-Abhangigkeit der y,,, kommt nur in den rechten Seiten der ersten beiden Gleichungen 
(11) hereint. Aus diesem Grund sind die An, (x, 1) in Produkte g,(x)C,,,(7) aufspaltbar, 
wobei die g,(x) eindeutig durch die q(x) ausdriickbar und die Cn, (4) universelle, d. hh. vom spe- 
ziellen Zylinder unabhangige Funktionen sind. Man berechnet sie auf Vorrat und kann dann Bei- 
spiele durch einfaches Einsetzen der betreffenden Potentialgeschwindigkeit q(x) in die Aus- 
driicke fiir u, v und w behandeln. 

Da die rechte Seite von (11a) nur einfach q(x) enthalt, setzen wir an 


Hn o(%> 1) = 4(*) Cno(m) - (13a) 
Beim Einsetzen dieses Ausdrucks in die rechte Seite von (11b) bemerkt man, daB alle Terme 
den Faktor qq’ besitzen. Daher setzt man 


; Anr(xm) = 9(x) (x) Cnai(n) - (13) 
Die rechte Seite von (llc) erhalt damit Terme mit zwei verschiedenen Formen der x-Ab- 
hangigkeit, namlich mit Faktoren qq’? und q?q’’. Wegen der Linearitat der Gleichungen kann 
man superponieren: 

Hn a(x) = 9(x) 9g? (x) Sn,2a(77) + 4°(x) 9'"(%) Sn,20(7) « (13 ¢) 


Entsprechend ist aufzuspalten: 


Hn s(%o) = 4(%)q9(X) Sn, 3a) + 9°(*) q(x) a" (%) Sn, 30(9) + CH) (4) nae)» (13. d) 

Die Anzahl der je Stufe bendtigten Funktionen € wichst weiterhin rasch an. Es erscheint 
nicht sinnvoll, noch héhere Naherungen zu berechnen, da die Grenzschichtvoraussetzungen 
ohnehin nur fiir kleine Zeiten gelten, jedenfalls nicht mehr zu lang nach dem Beginn der Grenz- 
schichtablosung, weil dann das starke Anwachsen der Grenzschicht die Potentialstrémung ver- 
drangt und somit die 4uBeren Randbedingungen verandert werden. 

Uber die Konvergenz des Blasiusschen Reihenansatzes (10) ist nichts bekannt. Fiir prak- 
tische Zwecke ist er bei kleinen Zeiten jedenfalls brauchbar, da der Maximalbetrag der Beiwert- 
funktionen €,,, ,(7j) mit wachsendem A, soweit berechnet, immer kleiner wird. 

Die Ansatze (13) ergeben wiederum eine Folge rekursiv zu losender Differentialgleichungen 
fiir die ¢,,, (Striche bedeuten hier Ableitungen nach dem Argument 7): 


mo +2 luo —4nlio =—A4n, (14) 
Ont +2 Cr. —4(2n+1) Sir =4 Cro —Cnolno — 1), (15) 
Cn. 2a fe 2 4 Cn2a— 4 (3 Aa 2) Cote = 4 (2 ape CA — bnobni= Cho ei) \ 


(16) 


Cn, 26 at 2% nob — 4 (3 UA 2) eae =4 Cn0Cn1 — bro Gai): 


Osa + 2 Cr3a—4 (40 4+ 3) Crise = 4 (2 Goo Cages te Or = Cao On,20. 57 Sn 0 Onda 
"— ni Cni)s | 
Orn +2 6% 35 —4 (404 3) Case = 4 (2 Co oCnte + 3 bn 0 lm2b + Sni— CnoSn,2b Co AD) 
— 2 Oho bn,20 — 2 Sno Sn,26— Cn Sn)» 
O50 +2 Cr se—4 (40 4+ 3) lise = 4 (Croln,2o — Or Cn,28) 


ere e reese eee eeee 


(bar: . _{[ono1, 
EE eT ERNE 


Die Methode zur Gewinnung einer Folge gewohnlicher Differentialgleichungen fiir die Bei- 
wertfunktionen der Querkomponente v ist mit geringen Abweichungen die gleiche. 

Damit in groBer Entfernung vom Zylinder eine von x und y unabhangige stationare Potential- 
strémung besteht, muB das Bewegungsgesetz fiir U? und fiir V° die gleiche Zeitabhangigkeit 


(18) 


1 Ist q = konst., so gibt bereits die nullte Naherung, wie wir sehen werden, die volle Lésung. Dies ist der 


- Fall der Anfahrt einer ebenen Platte in ihrer eigenen Ebene. 


15 
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aufweisen: 
{ 0 fur +< 0, 
V(t) =| Ge fiir 10 (@ = konst.). (19) 
Das daraus berechnete Druckglied wird in die Grenzschichtgleichung (5b) eingesetzt. 
Diese lautet in den auch fiir (9) eingefiihrten (x, 7, t)-Koordinaten 
n w Biya 1 
eR Maree gaia st mick igs Yan (20) 


Die Gleichung fiir v ist also im Gegensatz zu denjenigen fiir u und w von vorneherein linear. 
Wir machen den Zeitreihenansatz 


w(t Ast) == t" D> Eee Do aixe i) (21) 
Rez 


Da wir in (10) den entsprechenden Ansatz fiir die Stromfunktion gemacht hatten, miissen wir 
erwarten, daB die Beiwertfunktionen q,,, vergleichsweise eine Differentiationsordnung tiefer 
stehen als die y,,,, d.h. daB die w,,, den 07,,,/0n entsprechen. 

Mit dem Ansatz (21) fiir v sowie den sich aus (10) durch Differentiation ergebenden Aus- 
driicken fiir u und w gehen wir in (20) ein und erhalten die Differentialgleichungen 


te, ee Bae oe (22) 
a 2q Pt _ 42m +1) Peete (22b) 
oe 42n BB os (3n+2)O2=+4 ogee (Gieee cc ee oe » (22¢) 
a ap — (4n +3) Qm3= +4 fag'Cnn 2c8 Tee: 
— [q+ 2909) on20 + (24g q" + ea!" \bn,20) (22d) 
AG eect a er Cre ae : 
- a eee werden entsprechend gefordert 
fiir 7 = 0s\0p,9.— 0. fiir pspon BELG Geen, (23) 


Dann zerfallen wir die @,,, in gleichfalls zwangslaufiger Weise gemah 
Qn o(%s7)) = 4 Cno(M) » 
Pna(#o7) = 44 Salt) is (24) 
n2(%> 1) =F 1g? Cn,2a(7) +99 Sn,20(1) 15 
Pn (X71) = 714? Ln, 30m) + 99 4Cn,s0(0) + 24” En, s(t). 


Daraus folgen die gewéhnlichen Differentialgleichungen fiir die universellen Beiwertfunk- 


tionen ¢,,, der Querkomponente v 
Cro +2 Cro —Anl,» =—4n, (25) 
Car +29 Cir —4(2m +1) Snr =—4lno eo, (26) 
xe +20 Ca2e—4 (80 +2)bn20=—4 (Cla tar + ono loads | 
Coos +2 Crp —4 Bm 42) ane = £4 Cobo — Ceo ads | 


(27) 
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C3 + 2n Capes (4 n + 3) Cis — = 4 (Ciro Gata +on1 ae 1 Cro Ges 
Crab +27 Cn35—4 (4n + 3) C35 = +4 ee tage Cee wes Caine 


7 wr es = (28) 
4 —2fi0 Cn,2a— 25no Cn,2b — Cn Cni— no Cn,2b)» 
n3e + 21 Cn3e—4(4n + 38) Cre = +4 (plo Cn, 2b) 
mit den Randbedingungen 
a => Cat >1 
fir 7 = 0% Cz,=—0; fiir 17> oo: lz - 


Durch Vergleich von (25) mit (14) ersieht man unmittelbar: £9 = C/o. Diese Identitit ist 
in den rechten Seiten von (26), (27) und (28) bereits verwendet. 


Mit den Beziehungen (13) und (24) lautet nunmehr die Darstellung der gesuchten Strom- 
funktion und der Geschwindigkeitskomponenten der Grenzschichtstr6mung am Zylinder 


(x, 1, 1) = 2) vt q(x) {Lno(m) + tt aq’ Lar(y) + EO [42 Cy 20(n) + 94 Cn,20(N)] 


+ +) [4/3 f,34(n) + 99’ 4" Cn,39(N) oy) 
an qq’ fn,3-(n) |] -- o- ie 


1 , I gel 
(U(x, 7,1) = y= Siei” q(x) t” {2n0(m) + +1 q’Cni(m) + 2+ [9/2 2, 20(N) 
hey 3 

99 Cn,20(%) | 4PC*) [g*C,-20(7) ie 
+ 9q'q’'Cn,30(y) + 979" Cn,3e(4)] + °°}, 

o(x.9.t) = Ge" {Croln) + +8 q'Larln) + OCT) [G7Cn,20(n) +99 Cn, 20(7)] (a2) 

+ +1 [93 ln,30(n) +999 Sn30(1) + 49" Sn, se(m)] ++ * “hs 
w(x, 7, t) = — ys = —2V vt fq’ Cno(n) +t 8(G? + gq')Enrln) 
+ P+ DT(q's + 2 44/9") Suzan) + (2 49/9" + 474") Cn, 207) ] (33) 


) + BO+D[(q'* + 3.99/79) Cn, 309) + (249'2Q"’ +. 47g? +479 9") Cn, 50(n) 
) + (3 qq” + Pq) Snel) + °°} 


) Die u- und die v-Komponente stimmen in der Bauart véllig iiberein; sie unterscheiden sich 
‘nur in den Beiwertfunktionen. Wie die numerische Auswertung im nichsten Abschnitt zeigt, 
sind die Beitrage der ¢,., durchschnittlich geringer als die der entsprechenden Cray so daB 
demnach fiir gleiche Zeiten sich die u-Grenzschicht etwas weiter nach aufen erstreckt als die 
v-Grenzschicht. Fiir konstantes q(x) reduzieren sich die Reihen fiir u und v auf das erste Glied, 
wahrend w identisch verschwindet. Es wird zudem, wie zu erwarten, u~ v, und man erhalt 
eine der bekannten exakten Lésungen der Navier-Stokesschen Gleichungen, namlich den Fall 
der in ihrer Ebene angestrémten diinnen Platte, bei der x- und y-Richtung gleichwertig auftreten. 


4, Berechnung der universellen Beiwertfunktionen. Die universellen Funktionen ¢,,, bzw. 


€,,, wurden teilweise schon friiher berechnet und zwar Coo, Coy, C4) und C1, jeweils mit erster und 
-zweiter Ableitung durch Blasius}, 6,9, Cro. Cro fiir n =0,1,2,3,4 durch Gértler®, Co.2q und C25 


durch Goldstein und Rosenhead?. Fiir die C,,, lagen seit 1946 einige eigene unveréffentlichte 
Rechnungen geringerer Genauigkeit vor’ (graphische Integrationsmethoden). Die Arbeiten 
wurden damals nicht weiter verfolgt. ~ 

Wahrend Blasius und spater Goldstein-Rosenhead explizite Ausdriicke fiir ihre Funktionen 
-angaben, rechnete Gértler mit Hilfe numerischer Integrationsmethoden. Die Genauigkeit betrug 
_durchweg 3 bis 4 Dezimalen, teilweise jedoch bei sehr grofer Schrittweite in 7. 


1 H. Blasius, Z. Math. Phys. 56 (1908), S. 1. 
2 H. Gortler, Ing.-Arch. 14 (1944), S. 286. 
3 S. Goldstein und L. Rosenhead, Proc. Cambr. Phil. Soc. 32 (1936), S. 392. 
4 Vel. H. Gértler, Bericht tiber ,,Zahe Fliissigkeiten“ in ,,FIAT Review of German Science“ (1939—1946), 
Band 5, Teil III, S. 68. 
15% 
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Eine Auswertung der Ausdriicke (31) bis (33) am Beispiel des Kreiszylinders (siehe im fol- 
genden Abschnitt) zeigte, da ein Abbrechen der Reihen schon hinter den Gliedern mit ¢,,; die 
bei inkompressiblen Medien energetisch nicht sinnvolle Erscheinung einer ,,Uberhéhung“ der | 
Geschwindigkeitsprofile ergab. Es stand zu offen, daB® der Effekt durch héhere Reihenglieder 
wieder abgebaut werde. Daher schien es geboten, wenigstens ein weiteres Glied der zeitlichen 
Entwicklung zu ermitteln. Eine genaue Nachpriifung der Rechnungen von Goldstein- Rosenhead 
ergab leider verschiedene Fehler in den exakten Ausdriicken und infolgedessen auch in den 
numerischen Werten fiir Co.2_ und C62, (dort mit F, und f, bezeichnet) ?. 

Die prinzipielle Lésbarkeit der Differentialgleichungen (14), (15) und (16) durch exakte Aus- 
driicke war sichergestellt, und die Vermutung lag nahe, daf dies auch fiir die ibrigen in Frage 
kommenden Differentialgleichungen gelte. Da die Genauigkeit einer durch numerische Inte- 
gration erhaltenen Lésung infolge des Iterationsprozesses mit wachsendem / rasch abnimmt?, 
erwiesen sich die vorhandenen Tafeln der €,,, (A=0;1) als zu ungenau. 

Daher wurde der Weg der expliziten Lésung eingeschlagen, der jeden Verfahrensfehler ver- 
meidet. Die umfangreichen Rechnungen wurden zunachst nur fiir den Fall n=0 (ruckartige 
Anfahrt zu gleichférmiger Geschwindigkeit) durchgefiihrt. Es ergab sich dabei das jetzt zu 
schildernde verhaltnismaBig allgemein anwendbare Verfahren: 

Man geht aus von Gleichung (14) fiir den Fall n= 0 

Coo + 27 Coo = 0 (14") 
unter Beriicksichtigung der Randbedingungen (18) mit der Lésung 
b0(n) = aa et 
\x 


Durch einmalige Integration folgt 


Soul) = = [ ene ds. 


) 
Jetzt fihren wir die Bezeichnungsweise 
mot 
Te | e~* ds —1=@,(y) (komplementares Fehlerintegral) (34a) 
i 
und 
2 2 
——e—" = @,(n) (Fehlerfunktion) (34b) 


|x 
ein. ®)(7) und ®, (7) haben beide die Eigenschaft, im Unendlichen exponentiell zu verschwinden, 
so da auch nach Vorsetzen beliebiger Polynomkoeffizienten die Randbedingung fiir 7—>oo mit 
Ausdriicken in ®) und @, leicht zu erfiillen ist. 
Nunmehr erhalten wir durch nochmalige Integration unter Beriicksichtigung von (18) 
1 


Soult) = dln) +3 Ail) +n 


1 Ks sei hier erlaubt, die genannten Irrtiimer der Goldstein-Rosenheadschen Arbeit richti 
» die ¢ : t tellen. 
der Fehler auf S. 398 scheint indessen ein echter Rechenfehler zu sein. S. 395, Gl. (7), feicter Toad. ee 


ia 2 =-n2 a 2 ° _ -72 a 2 
2-1 (e T 1 @ Hig) lies 20° (e T 1 ¢@ 27) « S. 395, Ausdriicke (11), letzte Zeile: anstatt 
i eee 2 
ny => a? (20? + 1) erty + 6 


: 1 ij =n? 
lies n, = re ? (Qn? + 1) erfy + ne; S. 397, Ausdruck (25), letzter Term: anstatt a \/? (4 n> —10n) e 7 


: -1/2 -72 
lies a7 1/ (47? + 107)e 7; S.398, Ausdriicke (39), Term JT, erster Summand: anstatt ae [10 73 
0 


=F 1 

re aa) nates -1/2 -7? 

+ (20 — 150)7 — 8B] e” lies 30% / [10 0 4° + (15% — 20) 7 — 8B] e” ; S.399 Ausdriicke (44), Term, | 
j cia ; 8 

Soke Se anstatt A — Wm lies d = — eo Die Tafelwerte fiir f,(j) und F,(y) 

(einschl. graphische Darstellung) sind wegen des Fehlers auf S. 398 falsch. 


2 Es handelt sich um inhomogene Gleichun Bei lati i 
le ingen. Bei der Variation der Konstanten werden d ft 
Integranden fiir groBe 7 sehr groB, so das die absolute Genauigkeit der Loésung um jeweils yon noe ae 
potenzen geringer ist als die der rechten Seiten der Differentialgleichungen. 
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Man bemerkt die Eigenschaft von @, und @,, sich bei Integration und Differentiation mit 
gewissen Koeffizienten zu reproduzieren. 


Rechnen wir nun das inhomogene Glied der nachsten Gleichung (15) aus, so sehen wir, daB 
es aus Termen mit P?, D, D,, B® besteht; die wir als homogen vom Grade zwei bezeichnen wollen, 
ferner solchen mit ®, und ,, die entsprechend homogen vom ersten Grade sind, versehen jeweils 
mit Polynomkoeffizienten in 7. Man rechnet leicht nach, daB der Homogenitatsgrad der Terme 
bei Differentiation invariant ist. Der Differentialoperator in (15) ist linear; aus diesem Grund 
machen wir fiir die Lésung den Ansatz 


Cou(y) = (427? + 4p + Ay) Dy + (Bin + By) By ®, + C, BF 
+ (Dz 7? + Dyn + Dy) By + (E,9 + Ep) ®,. (35) 


Wegen €9:(0co) = 0 muB ein Absolutglied fehlen. 


Beim Einsetzen von (35) in (15) liefert der Koeffizientenvergleich zusammen mit einer Glei- 
chung zur Erfillung der Randbedingung (),(0)=0 gerade so viele lineare Gleichungen wie 
unbekannte Koeffizienten A,, B,,... vorhanden sind. Die Koeffizientenmatrix zerfallt jedoch 
in fast ganz getrennte Blécke auf der Hauptdiagonalen, fiir jeden Homogenitiatsgrad einen, 
so da die Behandlung verhaltnismaBig einfach ist. Das Resultat fiir €4;(7) ist unter (37) an- 
gegeben. 


Im nachsten Schritt ist £5; (7) zu integrieren. Auch hier empfiehlt sich ein Ansatz der Form (35), 
nicht etwa Teilintegration. Dabei stellt sich jedoch heraus, daf das System der Terme vom Homo- 
genitaétsgrad 2 offenbar noch nicht vollstandig ist. Man mu in den Ansatz fiir €),(7) noch ein 


Glied! mit @, (j2 7) aufnehmen, damit das entstehende lineare Gleichungssystem lésbar wird. 


Derselbe Vorgang, nur betrachtlich umfangreicher, wiederholt sich zur Bestimmung von 
Co2a(y) und fo 2,(7). Die Berechnung der rechten Seiten von (16) zeigt, daB ein Ansatz mit 
Termen bis zum Homogenitatsgrad drei notwendig ist. Diese sind die folgenden: 


Dn), Dr(y)®,(y) - Dy (4) B2(q), DN), Dy(3 7)» Bo(y2 y) Pn), Do(V 2 n) Py(n) - 


Ein solcher Ansatz fihrt durch Einsetzen in (16) auf ein System von 32 linearen Gleichungen 
fiir die verfiigbaren Koeffizienten A,, B,, ..., dessen Lésung die gesuchte Funktion €0,2.(7) 
bzw. Co,25(7) liefert. 


Es gelingt allerdings nicht mehr, mit dieser Methode auch Co2,(7) und Co,2,(7) zu gewinnen?, 
Sie wurden daher mittels sorgfaltiger numerischer Quadratur ermittelt. 


Genau das Analoge gilt fiir die Berechnung der £o,4(n) . Die linearen Gleichungssysteme unter- 
scheiden sich von denen fiir die £,,(7) nur durch geinderte inhomogene Glieder. 


Zur Wiedergabe der sehr langen expliziten Darstellungen fiir die Co,,(7) usw. wird nun die 
folgende unmi®verstandliche Symbolik verwendet, wobei die zur gleichen Naherung gehérenden 
Funktionen zunachst in gemeinsamer Form mit allgemeinen Koeffizienten angegeben sind; die 
zugehérige Tabelle gibt dann die Zahlwerte fir jeden einzelnen Fall: 


Cool) =n Be(n) + Pil) +1 aa ; 
(36) 


5 d ee geet 
1 @, (27) ist namlich vom Homogenitatsgrad zwei, wie a ®,(l 27) = 2 @?(n) zeigt. Wir rechnen 
somit alle Ausdriicke, deren erster Differentialquotient auf die Form P(y) ®j(n) ®{() mit» -- « = nzu bringen 
ist, zum Homogenitatsgrad n. : 


2 Das liegt daran, daB (7) Din) offenbar in Termen ®, nicht geschlossen integrierbar ist. 
Infolgedessen ist eine explizite Darstellung der Funktionen €),, C),, . . - gleichfalls nicht mehr moglich. 
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Colt) = (5 0? — 5m) @8) + (37? — 1a) Pol) Pale) + 3 Pin) 


12 gay +-|(1-+ shar +(— 3 +s9)9— gia] Oo 


le ! a | ( a i) Pao) : ( Poe ! st?) 


| 
I 


Cou(1)) | _ (An? + Ao) ®f(q) + Byrn Do(n) Pan) + Cy Pi(n) + (Dan? + Do) Po(n) 
Coun) J + (En + Ey) ®,(y) - 
Couln) £..(n) Con) foi(n) 
A int eed A aig er Bete 
: ue 1 ee andes ps be > 8c 
Ay a ey D BU il n 2 1 Ree 4, 
f: 3 1 0 2 Je 2 3a 
1 EX Wat pte e 
il 0 Ey cae ane 3 32 
Co Se 9 E, pale oe Siu 
3/2 3) x 
Soul) | _— Ar ®5(n) + (Ban? + Bo) Po(7) Pin) + Cry Pi(m) + Din Po(n) 
C01(7) J + (E, 7? + E,n + Ey) ®,(m)- 
Coin) AO) Coa(n) f1(n) 
A ie Ene 8 16 
1 nla eee ig a DE i Ora a eae 
B at see Ee Wa at 
M 2 2 A A 
1 1 E, a aieey a aT 
Ci —— es 3a 3) x 
2 2 fe 3, 4 Loge 
M ee Noe ae De 


(37) 


(38) 


Co,2a(7) und €o,,(7) sind nicht geschlossen angebbar (siehe vorhergehende Fufnote); 


(Ay 7* + A, 9? + Ap) O3(n) + (Bs 7? + By) O(n) Py(n) 
+ (Cy 4° +> Co) Po(n) Pi(n) + (D3 4? + Dy) y(n) 


£o,2a(1) + (Ey 7? + E, 7? + Ep) By (V3 n) + Fy D,( y2 n) ®,(n) 
F0,20(1) | Bs + (Gq7* + G7 + Gn + Go) Be(n) 
£o,2a(7) | & + (H, 7? + Hy 7? + Hy + Hy) ®o(7) ®,(n) 
Co,20(n) + (Ip? + En + Ip) ®i(n) 
+ (Ky* + K, 7? + Ky + Ky) ®,(7) 


+ (Lz 7% + Ly? + L1H + Ly) ®,(n). 


(39) 


Die hier vorliegenden Ausdriicke wurden numerisch ausgewertet! und zwar im Hinblick auf 
Verwendbarkeit zu eventuellen weiteren iterativen Verbesserungen mit méglichst hoher Stellen- 
zahl®?. Es wurde durchgehend mit 10 Dezimalen gerechnet; nach Abstrich der im Differenzen- 
schema angezeigten Rauhigkeit durch Rundungsfehler sind die nachfolgend angegebenen Werte 
auf eine halbe Einheit der letzten Dezimalen genau. Dasselbe gilt auch fiir die durch numerische 


Integration erhaltenen Funktionen €o,(7) und €o,22(1)- 


1 Die numerischen Rechnungen wurden von Frl. Herlinde Kompe ausgefihrt. 


2 Als Grundlage diente das Tafelwerk Probability Functions Vol. I, des National Bureau of Standards, 


Washington D. C. 


| 
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Tabelle 1. Die universelle Beiwertfunktion €,)(n) mit erster und zweiter Ableitung. 


Ui Coon) coo(n) Cooln) 
0 0 0 1,128 379 167 1 
0,1 0,005 632 511 4 0,112 462 916 0 1,117 151 606 8 
0,2 0,022 418 327 8 0,222 702 589 2 1,084 134 787 1 
0.3 0,050 028 899 1 0,328 626 759 5 1,031 260 909 6 
0,4 0,087 938 007 9 0,428 392 355 0 0,961 541 298 8 
0,5 0,135 451 644 8 0,520 499 877 8 0,878 782 578 9 
0,6 0,191 745 786 8 0,603 856 090 8 0,787 243 431 7 
0,7 0,255 908 682 4 0,677 801 193 8 0,691 274 860 4 
0,8 0,326 991 022 2 0,742 100 964 7 0,594 985 786 3 
0,9 0,404 012 094 7 0,796 908 212 4 0,501 968 574 2 
1,0 0,486 064 958 1 0,842 700 792 9 0,415 107 497 4 
1,1 0,572 275 791 9 0,880 205 069 6 0,336 479 597 8 
12, 0,661 859 363 8 0,910 313 978 2 0,267 344 347 0 
Hes 0,754 124 738 3 0,934 007 944 9 0,208 207 986 8 
1,4 0,848 480 438 0 0,952 285 119 8 0,158 941 707 7 
1,5 0,944 433 280 7 0,966 105 146 5 0,118 930 289 2 
FG 1,041 582 359 1 0,976 348 383 3 0,087 229 058 6 
ila 1,139 609 7163 0,983 790 458 6 0,062 711 040 5 
1,8 1,238 269 181 1 0,989 090 501 6 0,044 191 723 3 
1,9 1,337 374 602 1 0,992 790 429 2 0,030 524 740 4 
2,0 1,436 788 439 2 0,995 322 205 0 0,020 666 985 4 
etl 1,536 411 361 4 0,997 020 533 3 0,013 715 650 0 
oe, 1,636 173 232 1 0,998 137 153 7 0,008 922 155 1 
2,3 1,736 025 618 8 0,998 856 823 4 0,005 689 017 2 
2,4 1,835 935 807 4 0,999 311 486 1 0,003 555 648 7 
Deo) 1,935 882 178 6 0,999 593 048 0 0,002 178 284 2 
2,6 2,035 850 752 1 0,999 763 965 6 0,001 308 050 0 
25 2,135 832 680 6 0,999 865 667 3 0,000 769 924 8 
2,8 2,235 822 483 5 0,999 924 986 8 0,000 444 207 9 
2,9 2,335 816 838 0 0,999 958 902 1 0,000 251 210 9 
3,0 2,435 813 771 5 0,999 977 909 5 0,000 139 253 1 
ail Papssyatsy (3114 BYE I 0,999 988 351 3 0,000 075 663 3 
oe, 2,635 811 282 7 0,999 993 974 2 0,000 040 297 6 
Sh 2,735 810 844 7 0,999 996 942 3 0,000 021 037 2 
3,4 2,835 810 624 2 - 0,999 998 478 0 0,000 010 764 9 
sy 2,935 810 515 3 0,999 999 256 9 0,000 005 399 4 
3.6 3,035 810 462 6 0,999 999 644 1 0,000 002 654 6 
a7 3,135 810 437 5 0,999 999 832 8 0,000 001 279 3 
3,8 3.235 810 426 0 0,999 999 923 0 0,000 000 604 3 
3,9 3,335 810 420 7 0,999 999 965 2 0,000 000 279 8 
4,0 3,435 810 418 4 0,999 999 984 6 0,000 000 127 0 
4,1 3,535 810 417 2 0,999 999 993 3 0,000 000 056 5 
4,2 3,635 810 416 6 0,999 999 997 1 0,000 000 024 6 
4,3 3,735 810 416 6 0,999 999 998 8 0,000 000 010 5 
4,4 3,835 810 416 5 0,999 999 999 5 0,000 000 004 4 
4,5 3,935 810 416 5 0,999 999 999 8 0,000 000 001 8 
4,6 4,035 810 416 5 0,999 999 999 9 0,000 000 000 7 
4,7 4,135 810 416 5 1,000 000 000 0 0,000 000 000 3 
4,8 4,235 810 416 5 1,000 000 000 0 0,000 000 000 1 
4,9 4,335 810 416 5 1,000 000 000 0 0,000 000 000 0 
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Tabelle 2. Die universellen Beiwertfunktionen Cox(m)) und £y,() mit Ableitungen. 


Cox() 


0 
+-0,007 383 53 
+0,027 039 47 
+-0,055 507 79 
+ 0,089 733 23 


+0,127 091 32 
-+-0,165 401 25 
+ 0,202 924 41 
+0,238 348 37 
=+-0;270 757 27 


-++0,299 590 17 
+0,324 590 31 
+0,345 748 55 
+0,363 244 79 
+0,377 391 22 


-+-0,388 580 45 
+0,397 241 23 
+0,403 803 05 
+0,408 670 23 
-+0,412 205 11 


0,414 719 06 
410,416 469 87 
40,417 664 07 
10,418 461 80 
+0,418 983 73 


0,419 318 19 
0,419 528 12 
0,419 657 19 
+.0,419 734 91 
40,419 780 76 


+0,419 807 26 
+.0,419 822 25 
40,419 830 57 
0,419 835 09 
0,419 837 49 


+0,419 838 74 
+0,419 839 38 
+.0,419 839 70 
0,419 839 86 
0,419 839 93 


0,419 839 97 
0,419 839 99 
40,419 839 99 
+.0,419 840 00 
0,419 840 00 


of (7) 


0 
+-0,141 283 66 
0,246 060 94 
+-0,318 208 01 
+0,361 923 47 


+0,381 602 94 
-+-0,381 699 89 
+0,366 580 11 
+0,340 379 04 
-+0,306 871 99 


+0,269 367 16 
+0,230 629 95 
+0,192 843 86 
+0,157 609 73 
+0,125 980 39 


+0,098 524 61 
+0,075 411 01 
0,056 501 90 
0,041 447 09 
0,029 769 72 


0,020 938 23 
0,014 421 68 
+0,009 727 99 
0,006 426 55 
+0,004 158 09 


+0,002 635 02 
+0,001 635 53 
+-0,000 994 32 
+0,000 592 10 
+0,000 345 37 


0,000 197 33 
0,000 110 44 
40,000 060 55 
+-0,000 032 52 
-++0,000 017 11 


0,000 008 82 
0,000 004 45 
--0,000 002 20 
0,000 001 07 
£0,000 000 51 


+-0,000 000 24 
0,000 000 11 
+0,000 000 05 
0,000 000 02 
+0,000 000 01 


+0,000 000 00 


017) 


-+1,607 278 16 
+-1,224 170 49 
+0,877 853 47 
+0,572 106 00 
-+-0,309 552 01 


-++0,091 483 49 
—0,082 245 18 
—0,213 237 11 
—0,304 474 91 
—0,360 140 57 


—0,385 341 34 
—0,385 764 50 
—0,367 295 80 
—0,335 643 63 
—0,296 009 92 


—0,252 840 88 
—0,209 676 71 
—0,169 103 93 
—0,132 799 27 
—0,101 643 42 


—0,075 877 76 
—0,055 277 15 
—0,039 315 99 
—0,027 311 57 
—0,018 535 86 


—0,012 293 67 
—0,007 969 89 
—0,005 051 38 
—0,003 130 62 
—0 001 897 50 


—0,001 124 92 
—0,000 652 39 
—0,000 370 16 
—0,000 205 49 
—0,000 111 63 


—0,000 059 35 
—0,000 030 88 
—0,000 015 73 
—0,000 007 84 
—0,000 003 83 


—0,000 001 83 
—0,000 000 85 
—0,000 000 39 
—0,000 000 18 
—0,000 000 08 


—0,000 000 03 
—0,000 000 01 
—0,000 000 00 
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Coi(m) Coun) 

0 +0,170 581 18 
+0,017 093 83 +0,171 442 93 
+0,034 237 88 +0,170 808 11 
+0,051 066 69 +0,164 560 69 
+ 0,066 871 36 +0,150 048 47 
+0,080 771 39 +0,126 448 58 
+0,091 890 91 +0,094 763 95 
+0,099 536 55 +0,057 480 84 
+0,103 313 23 +0,017 989 91 
+0,103 182 14 —0,020 086 19 


+0,099 451 85 
+0,092 714 76 
+0,083 747 99 
-++0,073 402 30 
40,062 498 76 


+0,051 748 45 
-+-0,041 702 58 
+0,032 733 52 
-+-0,025 042 04 
+0,018 682 78 


+0,013 599 57 
0,009 663 00 
+0,006 704 58 
10,004 544 15 
0,003 009 66 


+0,001 948 04 
+ 0,001 232 77 
++0,000 762 85 
0,000 461 69 
+0,000 273 33 


-+-0,000 158 31 
+0,000 089 73 
+ 0,000 049 76 
--+0,000 027 01 
+0,000 014 35 


-+0,000 007 47 
+-0,000 003 80 
+0,000 001 89 
-++0,000 000 92 
-+-0,000 000 44 


+0,000 000 21 
+0,000 000 09 
--0,000 000 04 
+ 0,000 000 02 
+ 0,000 000 01 


-+-0,000 000 00 


—0,053 514 82 
—0,079 913 16 
—0,097 986 88 
—0,107 544 87 
—0,109 332 43 


—0,104 752 08 
—0,095 546 20 
—0,083 504 75 
—0,070 239 84 
—0,057 045 48 


—0,044 840 18 
—0,034 176 54 
—0,025 295 40 
—0,018 202 18 
—0,012 746 65 


—0,008 693 82 
—0,005 779 15 
—0,003 746 32 
—0,002 369 46 
—0,001 462 80 


—0,000 881 80 
—0,000 519 21 
—0,000 298 70 
—0,000 167 94 
—0,000 092 29 


—0,000 049 59 
—0,000 026 06 
—0,000 013 39 
—0,000 006 73 
—0,000 003 31 


—0,000 001 59 
—0,000 000 75 
—0,000 000 34 
—0,000 000 16 
—0,000 000 07 


—0,000 000 03 
—0,000 000 01 
—0,000 000 01 
—0,000 000 00 


Wundt: Wachstum der laminaren Grenzschicht an schrag angestromten Zylindern usw. 


Tabelle 3a. Die universellen Beiwertfunktionen C0,2a(yj) und Co,20(7). 


n €0,2a(7) 


0 
—0,001 245 5 
—0,005 025 9 
—0,011 412 0 
—0,020 416 0 


—0,031 937 8 
—0,045 742 2 
—0,061 461 5 
—0,078 618 0 
—0,096 660 0 


—0,115 005 5 
—0,133 086 7 
—0,150 390 9 
—0,166 492 7 
—0,181 074 9 


—0,193 937 7 
—0,204 995 5 
—0,214 2645 
—0,221 843 3 
—0,227 889 9 


—0,232 598 5 
—0,236 178 2 
—0,238 835 6 
—0,240 762 4 
—0,242 127 1 
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£0,20(7) 


0 
—0,000 338 2 
—0,001 359 1 
—0,003 057 5 
—0,005 390 7 


—0,008 271 3 
—0,011 574 2 
—0,015 150 6 
—0,018 843 9 
—0,022 505 3 


—0,026 005 3 
—0,029 242 1 
—0,032 1449 
—0,034 674 0 
—0,036 817 8 


—0,038 587 6 
—0,040 011 9 
—0,041 129 9 
—0,041 9865 
—0,042 627 4 


—0,043 095 8 
—0,043 430 3 
—0,043 663 8 
—0,043 823 1 
—0,043 929 4 
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£0,2a(7) 


—0,243 071 4 
—0,243 710 0 
—0,244 132 0 
—0,244 404 6 
—0,244 576 7 


—0,244 683 0 
—0,244 747 1 
—0,244 785 0 
—0,244 806 9 
—0,244 819 2 


—0,244 826 0 
—0,244 829 7 
—0,244 831 6 
—0,244 832 6 
—0,244 833 1 


—0,244 833 4 
—0,244 833 5 
—0,244 833 6 
—0,244 833 6 


Co,20() 


—0,043 998 8 
—0,044 043 1 
—0,044 070 8 
—0,044 087 7 
—0,044 097 8 


—0,044 103 7 
—0,044 107 1 
—0,044 109 0 
—0,044 110 1 
—0,044 1107 


—0,044 1110 
—0,044 111 2 
—0,044 111 3 
—0,044 111 3 


Ingenieur-Archiy 


AuBer den in der Tabelle zu findenden Zahlenwerten mégen noch die folgenden Anfangswerte 
explizit angegeben werden: 


B50) tees ae — + 1,128 379 1671, 
MA 


” 2 8 
0) = ees = + 1,607278 1594, 

ie 1 /1r 89 — 10873 256 
o0,20(0) = Tiga eee A = — 0,248 091 6725, 

i 1 7, —M+27/3 , 128\ | 

£0.24(0) =-( 4 jee s5ae) ~~ 090674203315, (40) 
ATC) ee aE URAL TLD or 

Vu 32)x 

, 1 23.0) 2451 Poo yiser 256 
C0,2a(9) el ae = | ee = + 0,0808175001, 
zr ed 23 , —253+351/3 512 
Co,20(0) = Ta ( ae fo ae =) = + 0,045 230 3896. 


5. Beispiel: Ruckartige Anfahrt eines Kreiszylinders. Die Ausdriicke (31), (32), (33) sollen 
nun zur Darstellung der Grenzschichtentstehung am schrag angestrémten Kreiszylinder ver- 
wendet werden und zwar, da die Beiwertfunktionen nur fiir n —0 vorliegen, fiir den Fall der 
ruckartigen Anfahrt des Korpers zu gleichformiger Geschwindigkeit. 

Die Potentialstrémung lings eines Normalschnitts eines Kreiszylinders ist bekanntlich 

fe x 
O\ X) sin (41) 
wo x die von der vorderen Staulinie aus gemessene Wandbogenlainge, U,, die Komponente der 
Anstromgeschwindigkeit parallel zur Normalebene, und R der Zylinderradius ist. Die V°-Kom- 
ponente ist natiirlich konstant: 


Gian (42) 
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Tabelle 3b. Die erste Ableitung der Funktionen €0,2 a() und £o,2 1(7) sowie die Funktionen f0,2a(n) und Co,20(n). 
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£0,2 a(n) 


0 
—0,025 000 1 
—0,050 739 5 
—0,077 011 2 
—0,102 905 7 


—0,127 138 5 
—0,148 333 6 
—0,165 250 6 
—0,176 952 4 
—0,182 910 3 


—0,183 046 3 
—O2177 718.0 
—0,167 653 6 
—(,153 851 5 
—0,137 460 0 


—0,119 655 0 
—0,101 531 0 
—0,084 017 3 
—0,067 826 1 
—0,053 433 6 


—0,041 089 9 
—0,030 850 2 
—0,022 619 1 
—0,016 198 1 
—0,011 3319 


—0,007 745 6 
—0,005 173 6 
—0,003 3773 
—0,002 155 0 
—0,001 344 2 


—0,000 819 8 
—0,000 488 8 
—0,000 285 0 
—0,000 162 6 
—0,000 090 7 


—0,000 049 5 
—0,000 026 4 
—0,000 013 8 
—0,000 007 0 
—0,000 003 5 


—0,000 001 7 
—0,000 000 8 
—0,600 000 4 
—0,000 000 2 
—0,000 000 1 


—0,000 000 0 


€0,25(n) C0,2a(n) £0,20(7) 
0 0 0 

—0,006 782 2 -+0,008 130 8 -+0,004 545 8 
—0,013 6315 -+0,016 331 0 -+-+0,009 081 6 
—0,020 267 1 +0,024 160 0 +0,013 287 4 
—0,026 249 6 +0,030 853 6 -+-0,016 709 7 
—0,031 1511 +0,035 583 4 +0,018 948 2 
—0,034 656 3 --+0,037 685 4 +0,019 781 0 
—0,036 607 7 +0,036 824 7 +0,019 2115 
—0,037 008 8 -++0,033 073 6 -+0,017 448 5 
—0,036 001 7 +0,026 897 4 +0,014 840 3 
—0,033 8285 -+0,019 061 7 --0,011 791 3 
—0,030 789 4 +0,010 489 2 -+-0,008 684 0 
—0,027 201 0 +0,002 102 1 +0,005 823 7 
—0,023 363 1 —0,005 316 2 +0,003 411 0 
—0,019 533 0 —0,011 2149 --0,001 539 5 
—0,015 912 1 —0,015 315 0 +0,000 213 4 
—0,012 639 6 —0,017 596 8 —0,000 6265 
—0,009 7965 —0,018 248 0 —0,001 074 2 
—0,007 412 3 —0,017 592 0 —0,001 2345 
—0,005 477 3 —0,016 012 2 —0,001 205 6 
—0,003 954 3 —0,013 887 1 —0,001 068 8 
—0,002 789 9 —0,011 545 5 —0,000 884 7 
—0,001 9241 —0,009 239 2 —0,000 694 0 
—0,001 2975 —0,007 138 3 —0,000 520 4 
—0,000 855 7 —0,005 336 8 —0,000 375 3 
—0,000 552 0 —0,003 867 8 —0,000 261 2 
—0,000 348 3 —0,002 721 0 —0,000 176 0 
—0,000 215 0 —0,001 860 4 —0,000 115 1 
—0,000 129 9 —0,001 237 3 —0,000 073 1 
—0,000 076 8 —0,000 801 1 —0,000 045 2 
—0,000 044 4 —0,000 505 3 —0,000 027 2 
—0,000 025 2 —0,000 310 7 —0,000 015 9 
—0,000 014 0 —0,000 186 3 —0,000 009 1 
—0,000 007 6 —0,000 108 9 —0,000 005 1 
—0,000 004 0 —0,000 062 2 —0,000 002 8 
—0,000 002 1 —0,000 034 7 —0,000 001 5 
—0,000 001 1 —0,000 018 9 —0,000 000 7 
—0,000 000 5 —0,000 010 0 —0,000 000 4 
—0,000 000 3 —0,000 005 2 —0,000 000 2 
—0,000 000 1 —0,000 002 7 —0,000 000 1 
—0,000 000 0 —0,000 001 3 —0,000 000 0 

—0,000 000 6 

—0,000 000 4 

—0,000 000 2 

—0,000 000 1 

—0,000 000 0 


Es ist zweckmaBig, noch die Reynoldssche Zahl Re=U,R/vy sowie eine dimensionslose 
Zeit t= U.,t/R einzufiihren. 7 mift den vom Zylinder zuriickgelegten Weg in Hinheiten seines 


Halbmessers. 
Bis zum quadratischen Glied angeschrieben lauten damit die Ausdriicke (31) bis (33) 


ule +> | = 2.U, sin | Coot) +27 cos Fe Foun) +47? [eos* Fe £0.20") | 


= sin" z £o,20() 


ions" 


a sin? = Co,2 Hn) | peg 


+ . 


| 


Cool) +2 r008-% foln) +47] cos? Loan) | 


: 


(43) 


(44) 
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w (i > r) ese vel | cos = Cool) + 27 (coe e = su x) Cox(") 
+473 [eos Ro = — sin? z) Co,2a(%) (45) 


— 3 cos z . sin > Co 20(n)| +: | ‘ 


Zwecks Vertafelung hat man sich nun auf einen bestimmten Anstrémwinkel tg «0 = V.,/U,, 
sowie eine bestimmte Reynoldssche Zahl festzulegen. Es wurde « =409°, d.h. Vy =U. (=e 
gewahlt. Im Hinblick auf die bei der experimentellen Druckmessung von Hiemenz! aufgetretene 
Reynoldssche Zahl Re =9360 wurde hier die runde Zahl Re =10000 zu Grunde gelegt. 


Abb. 1. Die Linien konstanter Steigung in der Potentialstrémung um einen unter ~ = 45° angestrémten Kreiszylinder. 
(Stromlinien gestrichelt.) 


Die Wiedergabe der Strémungsverhiltnisse durch ein anschauliches Bild bereitet wegen 
des Auftretens aller vier Koordinaten x, y, z,t betrachtliche Schwierigkeiten. Die Zeit wird in 
naheliegender Weise so beriicksichtigt, daB man fiir aufeinanderfolgende Zeitpunkte verschiedene 
Bilder zeichnet. 

Die Grenzschichtdicke mu} notwendigerweise in der Zeichnung passend tiberhéht werden. 
Dies hat zur Folge, da8 fiir die Zeichnung von Stromlinien auch die radiale Geschwindigkeits- 
komponente in diesem Mae zu iiberhohen ist, was der vorangegangenen Verkleinerung mit dem 
Faktor Re~!/? gerade entgegenwirkt. Man wirde in diesem Fall praktisch eine Strémung mit 
ziemlich kleinem Re zeichnen, jedoch als ob hierfiir die Grenzschichttheorie giiltig sei. Die dann 
am Rand der Grenzschicht durchaus merkliche w-Komponente kann nicht stetig — und da- 
mit nicht kontinuitatsgerecht — in die Potentialstrémung iibergehen, die ja dort W°—0 
verlangt. 


Diese prinzipielle Schwierigkeit, die fiir das Verstandnis der Senkrecht-Anstrémung eine Rolle 
spielt, wird hier umgangen, weil es sich in dieser Arbeit ja nur um die Darstellung des Effekts 
einer zusdtzlichen Querkomponente v handelt. Unter Beriicksichtigung der Tatsache, da8 man 
bei senkrechter Projektion der Stromlinien auf Normalschnitte des Zylinders wieder das alte 
Stromlinienbild der Senkrecht-Anstrémung erhalt, kann naimlich die Stromliniendarstellung 
iiberhaupt unterbleiben. Wir beschranken uns daher auf die Zeichnung eines Feldes von Linien 
gleicher Steigung (Isoklinen), welche die Stromlinien nunmehr gegeniiber Normalschnitten des 


1 K. Hiemenz, Dinglers Polytechn. Journal 326 (1911), S. 321. 
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Zylinders besitzen. Dies bedeutet, daB in der (Re 1) Bhene die Kurven 
eee ee kk 
are onst. = tg B (46) 
as 


zu zeichnen sind. Der Scharparameter f wird geeignet aquidistant variiert. 

Am Au®enrand der Grenzschicht mii®ten asymptotisch auch die Steigungsverhiltnisse der 
Potentialstrémung an der Kontur erreicht werden. Dies ist wegen der nicht verschwindenden 
w-Komponente, die nun aber nicht tiberhéht zu werden braucht, nicht genau erfiillt. Die 
starksten Abweichungen finden in der Gegend der Symmetrieebene der Strémung statt, weil 
dort u identisch verschwindet und daher 
w sogar uberwiegt. Lediglich aus Ver- 
standnisgriinden wird daher das ent- 
sprechende ,,Steigungsfeld“ fiir die ge- 
samte Potentialstrémung und den unter 
« = 45° angestrémten Kreiszylinder 
wiedergegeben. Hierbei kénnen die 
Stromlinien mit eingezeichnet werden, 
da keine Uberhéhungen notwendig sind! 


(Abb. 1). 


Zur numerischen Auswertung von 
(43) bis (45) wurde ein Gitter mit den 
Maschenweiten 


x ° ° ° 

aie 0° (10°) 180°, » = 0 (0,2) 4,0, 
a= 0.(0,1) 05 

zu Grunde gelegt?. 


Daraus wurde das Steigungsfeld (46) 
berechnet und durch die Gitterpunkte 
hindurch mittels graphischer Interpola- 


tion in «/R- und in y7-Richtung die a 
oki; it P M a oe ee ae ee ean & 
soklinen mit Farameterwerten 0 20° 40° 60° 0° mo’ = tea’ O° 160° = 180° 
B ae 0° (5°) 90° Abb. 2. Einige », Wandstromlinien’‘ der Potentialstrémung in Seitensicht. 

cape Die Pfeile auf der Mantellinie x/R = 90° geben die Steigung « = 45° der 


3 ungestérten Anstrémung an. 
gelegt. 
Die Grenzschichtdicke (aber nicht etwa w) erscheint in den Abbildungen gegeniiber den 
Zylinderdimensionen etwa 20-fach ttberhéht. Entsprechend dem Grenzschicht-Dickenwachstum 


1 Die Linien gleicher Stromliniensteigung V../U2 + W2 = konst. = tg sind hier zugleich auch Linien 
gleichen Geschwindigkeitsbetrags |/ U2-+ W? und damit nach Bernoulli gleichen Drucks; nur der Scharparameter 
hat jeweils verschiedene Bedeutung. Ihre Gleichung in kartesischen Koordinaten (x die Anstrémachse, 
Mittelpunkt des Zylinder-Normalschnitts im Ursprung) lautet fiir den Zylinder vom Halbmesser 1 
ogee 2 
pete ae) 
(x2 + 2)? 
Es handelt sich also um algebraische Kurven vierter Ordnung, wahrend die (gestrichelten) Stromlinien selbst 
mit der Gleichung 


= konst. = ctg”fB . 


1 
¥y ( ae al konst, = 6 
Kurven dritter Ordnung sind. Insbesondere zerfallt die Kurve mit b = 1 (d. h. B = 45°) in eine gleichseitige 
Hyperbel und den unendlich fernen Kreis. Die Steigungsverhaltnisse der durch keinen Korper gestorten An- 
strémung, also der Wert 4, finden sich an der Zylinderkontur auf Mantellinien, die von der (vorderen bzw. hin- 
teren) Staulinie um den Zentriwinkel 30° entfernt sind. 

2 Bei einer vorlaufigen Auswertung von (43) bis (45) nur bis zu Gliedern mit 1" einschlieBlich (vel. Fub- 
note 4 von S. 217) zeigte sich der bei inkompressiblen Medien vermutlich energetisch unzulassige Effekt einer 
Geschwindigkeitsiiberhéhung der u- und w-Profile auf der Zylindervorderseite. Darunter wollen wir ver- 
stehen, da innerhalb der Grenzschicht héhere Geschwindigkeiten auftreten als am Aufenrand asymptotisch 
erreicht werden sollten. Die Uberhéhung wird im Rahmen der vorliegenden Genauigkeit durch das folgende 
Reihenglied mit t? wieder abgebaut. Damit ist die zu Beginn des Abschnitts 3 behauptete Notwendigkeit 
genauerer Berechnung der Geschwindigkeitskomponenten begriindet. 

3 Die numerischen Rechnungen wurden von Frl. Herlinde Kompe ausgefibrt. 
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eo yr wird der Umfang der Bilder fiir wachsende Zeiten gréBer. Fiir tr —0 existiert noch keine 
srenzschicht. Man kann lediglich Wandstromlinien zeichnen, auf denen die Geschwindigkeit 0 
nerrscht und deren Steigung 

™ | 


tg By) = lim > — 


n70 Us Un 


(47) 


n= 0 
st. Da sie die potentialtheoretischen Steigungsverhiltnisse am AuBenrand der Grenzschicht fiir 


lle z angibt, sei ibr Verlauf in Seitensicht auf dem abgewickelten Halbzylinder (i Phen) 


wiedergegeben (Abb. 2). Ihre Gleichung ist Hee: 
{ x Ve Vii \ a1 
ve) = | iD a( 4 5 In tg 7 + konst. (48) 
Wie schon oben (Abschnitt 2) bemerkt, werden die Ablésungsorte und -zeiten der Grenz- 
. %, schicht durch die Querkompo- 
nente v nicht verandert. Es mége 
a y daher geniigen, einige Zahlwerte 


hierzu anzugeben, die aus der 
verbesserten Berechnung der Co,, 
resultieren. 


Abb. 8. Die Strémungsverhiltnisse im Wirbel- 
gebiet (schematisch). Die Stromlinien sind ge- 
< strichelt. a) Ablésepunkt, b) hinterer Staupunkt, 
f ° c) Wirbelpunkt, d) weiterer Punkt mit verschwin- 

Abb. 7. + = 0,5. denden Geschwindigkeitskomponenten u und w. 


In der hier vorliegenden zweiten Naherung ergibt sich der Ablésungsbeginn in der hinteren 
Staulinie zur Zeit 7 =0,319 507. Die Ablésestellen liegen fiir 7 = 0,4 bei = = 143,95°, fiirr =0,5 
bei — = 131,21°). 

Interessant mag auch noch der (natiirliche, zeichnerisch nicht iberhdhte) Ablésewinkel y, 
1.h. der Winkel zwischen Ablésestromlinie und Kontur, sein: 

Se y 
poh ae sas u Run R \n=0 ” oe 
er hingt im Gegensatz zu Ablésezeit und Ablésestelle wesentlich von der Reynoldsschen Zahl und 
yom Zylinderradius ab. Fiir Re =10000 und R=10 cm ergeben sich in den betrachteten Be- 
reichen zeitlich langsam anwachsende Winkel von unter 10’, also durchweg sehr flache Ablésung. 

Die Darstellung der zeitlichen Entwicklung der Steigungsverhaltnisse in der Grenzschicht 
um den schrig angestrémten Kreiszylinder (Abb. 3 bis 7) 1a8t nun folgendes erkennen. 

Die Linien gleichen Stromlinienanstiegs, die bei Reibungsfreiheit innerhalb der diinnen wand- 
zahen Schicht noch praktisch radial verlaufen wiirden, erscheinen schon beit =0,1 an der Kontur 
yegen die vordere Staulinie zu abgebogen; es entstehen Steigungsprofile. Der auBere fast radial 
verlaufende Teil dieser Profile wird, als nicht mehr Grenzschichtcharakter tragend, gestrichelt 
wiedergegeben. 


1 Fiir t -> 00 wird bekanntlich asymptotisch eine Ablésestelle zwischen 109° und 110° erreicht. 


230 Wundt: Wachstum der laminaren Grenzschicht an schrag angestrémten Zylindern usw. Ingenieur-Archiv | 


Fir tr = 0,2 verstarkt sich das Abbiegen der Isoklinen. Vor dem Druckminimum beim Zentri- 
winkel 90° bildet sich eine wandnahe Zone mit Untersteigung aus, d.h. ein Gebiet, in dem ge- 
ringere Steigung herrscht als es die Potentialstrémungsverhaltnisse zulassen wirden. Die 
Isoklinen hinter dem Druckminimum zeigen bereits einen Wendepunkt, der auf die kommende 
Ablésung hinweist. 

Abb. 5 zeigt die Verhaltnisse bei r =0,3 kurz vor der Ablésung. Das Untersteigungsgebiet | 
hat sich ausgedehnt. Immer mehr Profile haben Wendepunkte. Die Zone starker Steigung hat | 
sich von der hinteren Staulinie aus an der Wand weit nach vorne verbreitert. 


Bei t = 0,4 hat die Ablésung schon stattgefunden. Der (in der Ausdrucksweise einer Ge n)- 


Projektion) hintere Staupunkt hat sich in sechs einzelne Punkte mit u=w—=0 aufgespalten. | 
Dies sind Orte mit 90°-Steigung, also einer Strémung parallel zu den Erzeugenden. Man liest | 
aus Abb.8, die die charakteristischen Linien im Wirbelgebiet schematisch wiedergibt, die Lage | 
dieser Punkte ab. 
Bei rt =0,5 hat sich das Gebiet starker Steigung um die hintere Symmetrieachse bereits sehr / 
weit ausgedehnt. Wahrend zwischen Ablésepunkt a und Wirbelpunkt c noch eine Briicke 
ziemlich hoher Steigung bleibt, bildet sich im wandnahen Wirbelgebiet eine Zone geringerer | 
Steigung aus, die sich rasch eintieft. Raumlich gesehen sind die Stromlinien dort Schrauben- 
linien wechselnder Steigung. Im betrachteten Zeitabschnitt bis 7=0,5 ist die Steigung im 
gesamten Feld stets positiv. Fiir viel gréBere 7 lieBe sich, beginnend auf der hinteren Sym- | 
metrieachse, eine Abwartsstr6mung errechnen; doch ist dann wegen der Ausbreitung der Wirbel | 
die Grenzschichtvoraussetzung einer ungestérten Potentialstrémung langst nicht mehr giiltig. 
SchlieBlich mége darauf hingewiesen werden, daB in den Abbildungen auf der Zylindervorder- ; 
seite die Steigungsprofile nach auBen zu (gestrichelt eingezeichnet) wieder etwas absinken. Digs 
beruht jedoch nur auf der endlichen Gréfe von w, ist also ein unerwiinschter Effekt der Grenz- 
schichtvernachlassigungen, die, wie mehrfach betont, den kontinuitatsgerechten Anschlu8B der | 
Grenzschicht an die AuBenstrémung nicht gewahrleisten. 


(Eingegangen am 10. September 1954.) 
Anschrift des Verfassers: Dr. Hermann Wundt, Bensberg bei Kéln, Planstrafe 4 (bei der Ferd.-Stucker-StraBe). 


_ Die Grundlehren der mathematischen Wissen- 
schaften in Einzeldarstellungen mit besonderer Beriicksichtigung der Anwendungsgebiete, 


- Herausgegeben von R. Grammel, E. Hopf, H: Hopf, F. Rellich, F. K. Schmidt, 
B. L. van der Waerden. 
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Von Dr. Hermann Boerner, Professor der Mathematik an der Justus-Liebig-Hochschule, GieBen. 
Mit 15 Abbildungen. XT, 287 Seiten Gr.-8°. 1955. DM CY ane Ganzleinen DM 36,60 
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Ergebnisse der Mathematik und ihrer Grenz- 


pe ebiete. Unter Mitwiskung der Schriftloitung des ,,Zentralblatt fiir Mathematik® 
herausgegeben von L.V. Ahlfors, R. Bacr, R. Courant, J. L. Doob, S. Eilenberg, H. Rademacher, 
F. K. Schmidt, B. Segre, E.Sperner. Neue Folge. 


Die Ergebnisse der Mathematik fanden in der Zeit ihres Erscheinens von 1932 bis 1942 viel Beachtung 
und weite Verbreitung. Sie beginnen'jetzt wieder in neuer Folge in einzeln berechneten Heften und werden 


? 
‘ 


schnell fortgesetzt werden. 


Die Bezieher der ,,Zentralblatt fiir Mathematik“ erhalten die ,,Ergebnisse der Mathematik‘ zu einem 


ad 


gegeniiber dem Ladenpreis um 10% ermaéBigten Vorzugspreis. 
Erstes Heft: Transfinite Zahlen. Von Heinz Bachmann. VII, 204 Seiten Gr.-8°. 1955. DM 29,80 


transfinite Funktionen. — III. Arithmetik der Ordnungszahlen. — IV. Arithmetik der Machtigkeiten und Kardinalzablen 


VI. Probleme des Kontinuums und a zweiten Zahlklasse. — VII, Unerreichbare Zahlen. —Literatu rverze ichnis. | 
— Sachverzeichnis. aie 


Zweites Heft: Equazioni alle derivate parziali di tipo ellittico, Di. Carlo Miranda. 222 Seiten Gr,-8°. 


Indice: Prefazione. — I. Problemi al contorno per le eqyazioni lineari. — II. Funzioni rappresentate da integrali. — III. Tra- 
duzione in equazioni integrali dei problemi al contorno. — IV. Soluzioni generalizzate dei problemi.al contorno. — V. Mag-— 
giorazione a priori delle soluzioni del problema di Dirichlet. — VI. Equazioni non lineari. — VII. Altre ricerche sulle equazioni 
ellittiche. Equazioni di ordine superiore. Sistemi di equazioni. — Bibliografia. — Indice per autori. 


und dabei auch eine gedringte Ubersicht iiber die verschiedenen, beim Studium der partiellen Differentialgleichungen an- ‘ 

gewandten Verfahren zu geben. Er beschiftigt sich hauptsichlich mit den grundlegenden TeiJen der Theorie, die im Gegensata 

zu den Untersuchungen tiber gemischte Probleme, iiber Gleichungen héherer Ordnung und Systeme: von Differentialgleichungen 

schon einen gewissen AbschluB erreicht haben. Die folgenden Verfahren werden besonders herausgestellt: 

1. Die BrrauDsche Methode fiir die Zuriickfihrung der Randwertprobleme auf  Integralgleichungen zweiter Art. 

2. Das Studium der Randwertprobleme mit Hilfe linearer Funktionalanalysis, ‘das in Arbeiten von CACOIOPPOLI, PICONE 
und. WHYL entwickelt worden ist. ; 

3. Die Untersuchungen von BERNSTEIN, HOPF, SCHAUDER, LERAY und CACOIOPPOLI, die auf direkte Majorisierung’ 
der Lésung linearer und nichtlinearer Gleichungen hinauslaufen. — ® 

In zwei vorbereitenden Kapiteln werden verschiedene klassische Grundbegriffe sowie Bigeusshatien verallgemeinerter Po- 

tentiale zusammengestellt. SchlieBlich werden im 7, Kapitel einige spezielle Fragen bei einer Gleichung zweiter Ordnung, bei 

Gleichungen héherer ‘Ordnung, bei Gleichungssystemen sowie die Abhangigkeit von einem Parameter untersucht. Das Heft 

schlieBt mit einer ausgedehnten Bibliographie, die mehr als 600 nach 1924 erschienene Arbeiten umfaBt. > 


’ 


Drittes Heft: Analytische Fortsetzung. Von L. Bieberbach, Berlin, Mit 2 Textabbildungen. IV, 
168 Seiten Gr.-8°. 1955. ; =e ME 24, 80 


“HURWITZ und CRAMER. Die EULERsche Reihentransformation. Ein Test fir singulare Stellen eines Funktionselementes. ‘ 
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Liickensatz von POLYA. Weiteres iiber Koeffizientendichte. Komplementaére Reihen. — Die Haufigkeit der fortsetzbaren 


sitze zum HADAMARDschen Multiplikationssatz. Altere Untersuchungen. Die neuere Entwicklung. — Arithmetische Eigen- 
schaften det Koeffizienten. Potenzreihen mit endlich vielen verschiedenen Koeffizienten. Potenzreihen mit ganzen rationalen 


gemeiner Satz vod LHAU. Der spezielle Satz von Lpav. —Literaturverzeichnis,—Namenverzeichnis. 
— Sachverzeicthnis. 


Als letzter Band der alten Reihe erschien am 30. 11. 1942: ee 
: V. Band, 5. Heft: Nagy, Spektraldarstellung linearer Transformationen des Hilbertschen Raumes. 


SPRINGER-VERLAG / BERLIN. 66 OTTINGEN. HEIDELBERG 


Diesem Heft liegen 2 Prospekte des Springer-Verlages, Berlin | Géttingen / Heidelberg und 1 Prospeht 
des Springer-Verlages, Wien, bei. 


Springer-Verlag, Berlin / Géttingen / Heidelberg. Printed in Germany. 57 275 4022 0,8 K.B./Z. 040. 
. | | 


Inhaltstibersicht: I. Einleitung. Allgemeine mengentheoretische Vorbemerkungen. — II. Ordnungszablen und © 


ohne Auswahlaxiom, — V, Die Konsequenzen des Auswahlaxioms und der Alephypothesein der Kardinalzahlenarithmetik, — ia 
1955, . | DM 28,80 


Dem Verfasser erschien es daher von Interesse, den gegenwirtigen Stand der Forschung in einerMonographie zusammenzufassen 
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Inhaltsiibersicht: Grunidlegende ses ote Die LAPLACE-BORELsche Transformation. Die Koeffizienten als ganze = Hs 


und der nichtfortsetzbaren Reihen. BOREL, STZINHAUS, BOERRNER, POLYA, HAUSDORFF. Banach-Raume. — Zu- 


